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PRELIMINARY REMARKS 



1. Measurable sets. Let E be some set of points of an 
interval 5 = [a, b]. Denote the complement of E with respect 
to 5 by C E ; i.e., by definition C B consists of points which 
do not belong to E. 

There are a variety of ways in which the points of set 
E may be included in a finite or countable system of intervals 

(Xj, a 2 , . . . , a„, . . . 

We denote by 2 a the sum of the lengths of the intervals 
a,, a 2 , . . ., a n , .... For any system of intervals covering E, 

2«>o 

The lower bound of 2 a > which depends solely on the 
set E, is called the exterior measure and is denoted m*E. 
From the definition of an exterior measure it follows that 
for any e > there exists a system of intervals a x , a 21 ... 
. . ., a n , ... which include all points of the set E such that 

The interior measure m M E of the set E is the difference 
between the length of the interval S and the exterior measure 
of the complement of the set; i.e., 

m,£ = b — a — m* C E 

If the exterior and interior measures of E are equal, then 
the set E is called measurable in the sense of Lebesgue 
(Lebesgue measurable, or, simply, measurable), while the 
common value of the measures m* E and m^E is called the 
Lebesgue measure of E (or, simply, the measure of E) and 
is denoted by mE or mes E. 

The measure of the interval (a, b) is its length: mes 
(a, b) — b — a. The set co of points of the interval (a, b) is 
called a set of measure zero if co can be covered by intervals 
the sum of whose lengths is arbitrarily small. 
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2. A function of a real variable f(x) defined on a measu- 
rable set E is called measurable if for any number A the 
set <$> (/ (x) > A), which consists of those points x belonging 
to the set E for which f(x)>A, is Lebesgue measurable. 

Note. The requirement of measurability of the set 
& (/ (*) > A) may be replaced by one of the following three 
conditions: 

(a) the set $ (f (x)^ A) is measurable, 

(b) the set d> (f (x) < A) is measurable, 

(c) the set <§ (f (x) <: A) is measurable. 

3. A function f(x), nonnegative on the interval (a, b), 

b 

is called summable on that interval if J / (x) dx is finite*. 

a 

A function f (x) of arbitrary sign will be summable on 
an interval (a, b) if and only if the function \f(x)\ is 

b 

summable, i.e., if the integral ^ 1 / (jc) | has a finite 

a 

value. 

In the sequel we shall have to do with the basic interval 
/ = (a, b) (or /„ = (0, a) ) and the basic square 

fi{a<jM<6}(or S {0 < jc, / <a} ) 

4. The L % (a, b) space. We say that f(x) is a quadratically 
integrable function on [a, b] if the integral 

b 

$f*(x)dx 

a 

exists (is finite). The class of all quadratically integrable 
functions on [a, b] is denoted by L 2 (a, b) or, simply, L 2 . 

Basic Properties of L 2 Functions 

(a) The product of two quadratically integrable functions 
is an integrable function. 

(b) The sum of two L 2 functions is also an L 2 function. 



* The integral is to be understood in the sense of Lebesgue, but 
if the reader is not acquainted with the Lebesgue integral, the integrals 
may everywhere be understood in the sense of Riemann. 
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(c) If f(x) g L 2 and X is an arbitrary real number, then 

lf(x) <= L % 

(d) If / (x) g L % and g (x) g L 2 , then we have the Bunya- 
kovsky-Schwarz inequality 

( \f(x)g(x)dx) *zlp(x)dxl?(x)dx (1) 

\ a J a a 

The scalar product of two functions f(x) € ^» and g(x) g L 2 
is, by definition, the number 

b 

(/, g)= \f(x)g(x)dx (2) 

a 

The norm of an L a function / (x) is the nonnegative number 



Ul = V(f,f) = y/ r \f*(x)dx (3) 



(e) For /(*) and g(x) taken in Z, 2 we have the triangle 
inequality 

(f) Convergence in the mean. Let the functions f{x) and 
fi (*)> fi ( x )> • • • > fn (•*)> • • • be quadratically summable on 
(a, b). If 

6 

lira Uf n (x)-f(x)]*dx=*0 

tlicn we say that the sequence of functions f^x), f 2 (x) ... 
converges in the mean or, more precisely,, in the mean square, 
to the function f (x). 

If a sequence {/„(#)} of L a functions converges uniformly 
to f(x), then f(x)£L 2 and {/„(*)} converges to f(x) in the 
mean. 

We say that a sequence {/„ (x)} of functions in L 2 converges 
in the mean in itself if for any number e > there exists 
<i number N > such . that 

b 
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for n > N and m> N. Sometimes sequences convergent in 
themselves are called fundamental sequences. For a sequence 
{fn( x )} to converge in the mean to some function, it is 
necessary and sufficient that this sequence be fundamental. 
The space L 2 is complete; i.e., any fundamental sequence 
of functions in L 2 converges to a function which also lies 
in L 2 . 

Two functions f(x) andg(jc) in L 2 (a, b) are called equivalent 
on (a, b) if f(x)=^g(x) only on a set of measure zero. 
In this case we say that f(x) = g(x) almost everywhere on- 



5. The space C U) (a, b). The elements of this space are 
all possible functions defined on the interval [a, b] and 
having, on this interval, continuous derivatives up to the 
/th inclusive. The operations of addition of functions and 
multiplication of functions by a number are defined in the 
usual manner. 

We determine the norm of an element f(x) g C U) (a, b) from 
the formula 



f i0) (x) being equal to f(x). 

Convergence in C il) (a, b) implies uniform convergence 
both of a sequence of the functions themselves and of the 
sequences of their derivatives of order k (fe=l, 2, /). 

The concepts of a measurable set, a measurable function, 
a summable function, etc. are extended to the case of spaces 
of higher dimensionality. For example, the function F (x, t) 
will be called quadratically summable on £2{a^x, t ^.b] 
if the integral 



(a, b). 




b b 




a a 



In this case, the norm of the function F (x, t) is defined 
by the equality 
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6. A function /(z) of a complex variable z, differentiable 
;il every point of a domain G in the plane of the complex 
variable z, is called analytic (regular) in that domain. 

The function f(z) is called entire (integral) if it is analytic 
throughout the plane (with the exception of the point at 
infinity). 

The function f(z) is called meromorphic (or fractional) if 
it can be represented as a quotient of two entire functions: 



In any bounded domain, the meromorphic function /(z)can 
have only a finite number of poles. 

A point z = a is called an isolated singular point of 
the function / (z) if there is a neighbourhood < | z — a | < 8 
of that point in which f (z) is analytic, while analyticity 
of the function breaks down at the point z = a itself. The 
isolated singular point z = a is called a pole of the func- 
tion f(z) if 

lim f(z) = 00 

It is assumed that f(z) is single- valued in the neighbour- 
hood of the point z = a, z=£a. 

For the point z — a to be a pole of the function f (z) 
it is necessary and sufficient that it be a zero of the 

function <p (z) = j^y , i.e., that qp(a) = 0. 

The order of a pole z — a of the function f(z) is the 
order of the zero z = a of the function 

7. The residue of the function /(z) at the isolated sin- 
gular point z = a is the number 



^f/ (z)= 2irj /(2)d2 



where c is a circle |z — a| = p of sufficiently small radius. 
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If the point z = a is a pole of order n of the function 
f(z), then 

res / («) = j^jyy lim ^ {(* - a)» / (,)} 

For a simple pole (n=l) 

res / (z) = lim {(z — a) f (z)} 

If / (z) = -2i£L ( <p (a) 0, and, at the point z = a, (z) has 
a zero of order one, that is, t|> (a) = 0, i|/ (a) =^ 0, then 

res/(z) = ^ 

8. Jordan's lemma. If f(z) is continuous in the domain 
| z | ^ /?„, Imz^a (a is a fixed real number) and lim/(z) = 0, 
then for any %>0 z -** 

lim [ e iXz f (z)dz = 

where c R is an arc of the circle \z\ — R in that domain. 

9. A function /(jc) is called locally summable if it is sum- 
mable on any bounded set. 

Let a complex- valued function <p ( t ) of a real variable 
t be locally summable, equal to zero for t <0, and let it 
satisfy the condition |q> (t)\ < Me**' for all / (M > 0, s > 0). 
Such functions q>(*) will be called original functions. The 
number s is termed the order of growth of the function <p(/). 

T/te Laplace transform of the function <p(t) is the func- 
tion Q>(p) of the complex variable p = s + io defined by 
the equality 

CO 

®(p) = ^e-rtyWdt 
o 

For any original function <p(t), the function <t>(p) is 
defined in the half-plane Re p > s and is an analytic 
function in that half-plane. The fact that the function 
<D(p) is the Laplace transform of the function <p(/) is writ- 
ten as follows: 

<p(/)==<D(p) 
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10. Inversion theorem. If a function y(t) is the origi- 
nal function, and the function ^(p) serves as its image 
(transform), then 

V+ioo 

9(fl=-ST j eP«${p)dp,y>s (*) 

V — i oo 

where the integral is taken along the straight line Re p = y 
parallel to the imaginary axis and is understood in the 
sense of the principal value: 

V+ioo V+to 

\ eP* <P (p) dp = lim (" e^ <D (p) dp 

J. (0-»oo J 

V — i oo V — to 

Formula (*) is called the inversion formula of the Lap- 
lace transformation. If 

Q(p)= M(p) 
KP) N (p) 

where M (p) and jV (p) are polynomials in p, and the degree 
of the polynomial M (p) is less than that of the polynomial 
N (p), then for G> (p) the original function will be 

«p«) = £ T^rrn Jim ((/> - « ft ) n ^ (P) *"} 

where a ft are the poles of <D(p), n ft are their orders and 
the sum is taken over all poles of O(p). 
When all poles a k (& = 1, 2,..., /) of the function 

° (p) = "ffj" are sim P le - 



y M(a k ) (f) 
Nip) ■ Z* N'la k ) ^ ^ Kl) 



11. Product theorem (convolution theorem). Let the func- 
tions f (t) and y(t) be original functions, and let 

f(t) = F(p), 
<P(0 = «>(/>) 
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Then 



F(/7)(D (/>)=$ /(T)(p(*-T)dT 



(5> 







The integral on the right of (5) is called the convolu- 
tion of the functions / (t) and q> (t) and is denoted by the 
symbol f(t)*y(t). 

Thus, a product of transforms is also a transform, 
namely, the transform of the convolution of the original 
functions: 



12. Let the function f(x) be absolutely integrable over 
the entire axis — oo<jc< +oo. The function 



is called the Fourier transform of the function f(x). 

The inversion formula of the Fourier transformation is 
of the form 



In order to attain greater symmetry, the formulas of direct 
and inverse Fourier transformations are frequently written 
in the form 



F(p)<t>{p)=nt)*<p(t) 



f(X)= J f(x)e iXx dx 



— ao 



/(*) 



+« 




CHAPTER I 



VOLTERRA INTEGRAL EQUATIONS 



1. Basic Concepts 

The equation 

q> (*) = /(*) + t)y(t)dt (1) 

a 

where / (x), K (x, t) are known functions, <p(#) is the 
unknown function and K is a numerical parameter, is called 
Volterra's linear integral equation of the second kind. 
The function K (x, t) is the kernel of Volterra's equation. 
If 7(*)==0> then equation (1) takes the form 

X 

<?(x) = k^K(x, t)<?(t) dt (2) 

a 

and is called a homogeneous Volterra equation of the 
second kind. 
The equation 

X 

t)<p(t)dt = f(x) (3) 

a 

where q> (x) is the unknown function is called Volterra's 
integral equation of the first kind. Without loss of gene- 
rality, we can consider the lower limit a as equal to zero 
(in the sequel we shall assume this to be the case). 

A solution of the integral equation (1), (2) or (3) is a 
function (f(x), which, when substituted into the equation, 
reduces it to an identity (with respect to x). 

Example. Show that the function w(x)= 1 ,.„ is a 

. (l+**) 3/2 
solution of the Volterra integral equation 

X 
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Solution. Substituting the function 1 in place of 

6 (l + ^) 3/2 

y(x) into the right member of (4), we obtain 
_J P_J 1 A ; 

l+X* J l+X* (1+ /2)3/2 al ~ 

= jq^— T+f( _ (i+/t)>/2 )| <=fl 

75- = <P (*) 



= 

1,1 1 



1 (l +x I)3/2 !+«• (1+Jc2) 3/2 

Thus, the substitution of <p (*) = 1 ... into both 

sides of equation (4) reduces the equation to an identity 
with respect to x: 

l _ l 

According to the definition, this means that qp (x) = 
is a solution of the integral equation (4). 



(1+JC 2 ) 3/2 

Verify that the given functions are solutions of the 
corresponding integral equations. 

1. cp(x) = 



(l+^y 

. , 3x+2x 3 C3x+2x 3 — t ..... 



2. <p (x) = e* (cos e*— e* sine*); 

X 

q, (jc) = (1 — xe 2x ) cos 1 — e 2x sin 1 + J [l—(x—t) e tx ] <p(f) dt. 



X 

3. <p (x) = xe*; <p(x) = e*sinjt + 2 J cos(jc— t)q(t)dt. 



4. (p(x) = A; — -g-; y(x) = x— \ sinh(x— t)<p{t) dt. 
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X 

5. <f>(je)=l — jc; J e x ~ f (p(t)dt = x. 

o 

x 

6. q>(x) = 3; x 3 = J (x— tf <p (/)d/. 

o 

7. cp W =i.; 

o 

8. T (x) = -U; f^%*-l. 



Note. Volterra-type integral equations occur in problems 
of physics in which the independent variable varies in a 
preferential direction (for example, time, energy, etc.). 

Consider a beam of X-rays traversing a substance in the 
direction of the x-axis. We will assume that the beam 
maintains that direction when scattered. Consider a collec- 
tion of rays of specified wavelength. When passing through 
a thickness dx, some of the rays are absorbed; others un- 
dergo a change in wavelength due to scattering. On the 
other hand, the collection is augmented by those rays which, 
though originally of greater energy (i. e., shorter wave- 
length X), lose part of their energy through scattering. 
Thus, if the function f(X, x)dk gives the collection of 
rays whose wavelength lies in the interval from k to X + dX, 
then 



df (X, x) 



dx 



= -U./(A, X) + J?(X, T)/(T, X)dx 



where u. is the absorption coefficient and P (X, t)dx is the 
probability that in passing through a layer of unit thickness 
a ray of wavelength x acquires a wavelength which lies 
within the interval between X and X-\-dX. 

What we have is an integro-differential equation, i. e., 
an equation in which the unknown function f(X, x) is under 
the sign of the derivative and the integral. 



18 



INTEGRAL EQUATIONS 



Putting 

oo 

f(X, x) = $e-P*y(k, p)dp 



where p) is a new unknown function, we find that 

(X, p) will satisfy the Volterra integral equation of the 
second kind 

A. 



2. Relationship Between Linear Differential Equations 
and Volterra Integral Equations 

The solution of the linear differential equation 

S + M*)£S+--- +a n (x)y = F(x) (1) 

with continuous coefficients a t (x) (t'= 1, 2 n), given 

the initial conditions 

y(0) = C , y'(0) = C lt .... if-»(0) (2) 

may be reduced to a solution of some Volterra integral 
equation of the second kind. 

Let us demonstrate this in the case of a differential 
equation of the second order 

§ + aAx)f x + aAx)y = F(x) (V) 

y(0) = C„ </'(0) = C x (2') 

Put 

§ = <P(*) (3) 

Whence, taking into account the initial conditions (2'), we 
successively find 



% = ^<p(t)dt + C v y = ^(x-t)y(t)dt + C lX + C 



(4) 
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Here, we utilized the formula 

XXX X 

§dx§dx...§f(x)dx = J ^- i j T §(x-z)'>-if(z)dz 

X Xq X X 

n 

Taking into account (3) and (4), differential equat on (1) 
may be written as follows: 

<p (x) + J a x (x) q> (t) dt + Qa, (x) + J a, (x) (x—t) q> (I) dt + 



+ C,%a 2 (x) + C a 2 (x) = /=" (x) 

or 

X 

<p (x) + J [a, (x) + a 2 (x) (x— 0] 9 (t) dt = 

o 

= F(x)—C i a i (x)—C l xa 2 (x)— C a 2 (x) (5) 

Putting 

K(x, t)=-[a i (x) + a i (x)(x-t)} (6) 
f(x) = F (x)-C, ai (x)-C lX a 2 (x) -C a 2 (x) (7) 

we reduce (5) to the form 

X 

<p (*) = $*(*, t)<p(t)dt + f(x) (8) 





which means that we arrive at a Volterra integral equation of 
the second kind. 

The existence of a unique solution of equation (8) follows 
from the existence and uniqueness of solution of the Cauchy 
problem (l')-(2') for a linear differential equation with 
continuous coefficients in the neighbourhood of the point 
x = 0. 

Conversely, solving the integral equation (8) with K and / 
determined from (6) and (7), and substituting the expression 
obtained for <j>(x) into the second equation of (4), we get 
a unique solution to equation (!') which satisfies the initial 
conditions (2'); 
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Example. Form an integral equation corresponding to the 
differential equation 

y" + xy' + y = (1) 
and the initial conditions 

y(0) = 1, t/'(0) = (2) 

Solution. Put 
Then 

X XX 

^ = Jq)(/)d/ + /.(0)=f«p(0d/, y = j(*-/)<p(/)d*+l (4) 



Substituting (3) and (4) into the given differential equation, 
we get 

X X 



or 



q, (x) = — 1 — \ (2x— t) <p (/) dt 



Form integral equations corresponding to the following 
differential equations with given initial conditions: 

9. y" + y = 0; 0(0) = 0, 0'(O) = 1. 

10. 0'— = 0; 0(0)= 1. 

11. 0" + = cosx; i/(O) = 0'(O) = O. 

12. /_5/ + 64f = 0; 0(0) = 0, y' (0) = 1. 

13. tf + y = cosx; y(0) = 0, y'(0)=l. 

14. j/" — y' sin jc -(- 6 x y = jc; 0(O) = 1, 0'(O) =— 1. 

15. / + (l+tf)y = cosx; y(0) = 0, 0'(O) = 2. 

16. 0"'+x0" + (x a — x)y=xe*+l; 
y(O) = 0'(O) = l, /(0) = 0. 

17. -2x0 = 0; 0(0) = |, 0'(O) = 0"(O) = 1. 
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18. Show that a linear differential equation with constant 
coefficients reduces, under any initial conditions, to a Vol- 
terra integral equation of the second kind with kernel 
dependent solely on the difference (x — t) of arguments (integ- 
ral equation of the closed cycle or equation of the Faltung 
type, or convolution type). 

3. Resolvent Kernel of Volterra Integral Equation. 
Solution of Integral Equation by Resolvent Kernel 

Suppose we have a Volterra integral equation of the 
second kind: 

X 

<p(x) = f(x) + %lK(x,t)q>(t)dt (1) 



where K{x, t) is a continuous function for O^x^a, 
^ t ^ x, and / (x) is continuous for <: x ^ a. 

We shall seek the solution of integral equation (1) in the 
form of an infinite power in series X: 

op (x) = <p, (x) + Xq>! (x) + Pop, (x) + . . . + A»(p„ (x) + . . . (2) 

Formally substituting this series into (1), we obtain 

<Po (x) + hp l (x)+...+ X"cp n (x)+...= 

X 

= /(*) + * J /)[«p,(/) + X<p 1 (0+... 4-^,(0+...]* 



(3) 

Comparing coefficients of like powers of %, we find 

<Po (*) = /(*). 

x ■ X 

<p, (x) = S K (x, t) q>, (0 dt = \K (x, t) f (t) dt, (4) 
o o 

x x I 

<p 2 (x) = S K (x, t) (t) dt = lK(x,t)lK(t, t x ) f (tj dt.dt 





The relations (4) yield a method for a successive determi- 
nation of the functions y n (x). It may be shown that under 
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the assumptions made, with respect to f(x) and K(x, t), 
the series (2) thus obtained converges uniformly in x and 
X for any X and x£[0, a] and its sum is a unique solution 
of equation (1). 

Further, it follows from (4) that 



<p, (x) = J/C(jc, 0/(0* 



(5) 



cp 2 (x) = \K{x, t) \K(t,t l )f(t 1 )dt 1 



dt = 



= J / (t t ) dti J K (x, (t, t t ) dt = J K 2 (x, t t ) f (t t ) dt t (6) 

t, 

where 

X 

KAx,t,)=]K(x,t)K{t,t 1 )dt (7) 
<. 

Similarly, it is established that, generally, 

<P. 0/(0* (n = l,2, ...) (8) 
o 

The functions /(„(*, are called iterated kernels. It can 
readily be shown that they are determined with the aid of 
the recursion formulas 

KA*> t) = K(x, 0, 

x 

K n+1 (x,t) = ^K(x,z)K„(z,t)dz (n=l,2, ...) (9) 
t 

Utilizing (8) and (9), equality (2) may be written as 

<t(x) = f(x)+2k>)KAx,t)f(t)dt (10) 

v=I o 

The function R (x, t\ X) defined by means of the series 
R(x, t\%)= 2 X»K„Ax,t) (11) 

v=0 
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is called the resolvent kernel (or reciprocal kernel) for the 
integral equation (1). Series (11) converges absolutely and 
uniformly in the case of a continuous kernel K(x, t). 

Iterated kernels and also the resolvent kernel do not de- 
pend on the lower limit in an integral equation. 

The resolvent kernel R(x, t; X) satisfies the following 
functional equation: 

X 

R(x, t; l) = K(x, 1) + 'k\K(x, s)R(s, t\ %)ds (12) 

t 

With the aid of the resolvent kernel, the solution of 
integral equation (1) may be written in the form 

X 

cp {x) = f(x) + %\R{x, t; X)f(t)dt (13) 

(see [5], [15]). 

Example. Find the resolvent kernel of the Volterra 
integral equation with kernel K(x, /)=1. 

Solution.We have Ki(x, t) = K(x, t)=\. Further, by for- 
mulas (9) 



K 2 (x, ') = S K (x, 2)^(2, t)dz = \dz = x-t, 
t 

X 

K 3 (x, t) = ^\-{z-t)dz = [ - 

t 

KA*, t) = ^\. { ^^dz 



t 

(x-tf 



KAx, /) = ji .*._,(*. t)dz = l\.^^dz=^l 
i t 

Thus, by the definition of the resolvent kernel, 
R(x, t; X)= £ VK m+1 (x. i)= L ^j^W'*"" 

«= n = 
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Find the resolvent kernels for Volterra-type integral 
equations with the following kernels: 

19. K(x, t) = x—t. 

20. K(x, t) = e*-t. 

21. K(x, t) = e**-*\ 

22. K(x, = ^7? • 

23. tf(*./)-f±£f 

25. /((a;, = fl*"*(o>0). 

Suppose that the kernel K (x, t) is a polynomial of degree 
n — 1 in t so that it may be represented in the form 

K(x, t) = a % (x) + a l (x)(x-t) + . . . (14) 

and the coefficients a ft (*) are continuous ; n [0, a]. If the 
function g(x, t; X) is defined as a solution of the differen- 
tial equation 

g-X [«. (*)£^ + a, + . . . + W g] = (15) 

satisfying the conditions 

= 1(16) 



then the resolvent kernel # (*, K) will be defined by the 
equality 

R(x, (17) 

and similarly when 

/( (x, = MO + MO (*-*)+••• + ^rrjr (<-*)"• M 1 8) 
the resolvent kernel 

*(*. <; X) (19) 
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where g(x, t; X) is a solution of the equation 

g + k [b (t) +b n . 1 (t)g] = (20) 

which satisfies the conditions (16) (see [5]). 

Example. Find the resolvent kernel for the integral 
equation 

X 

<P(*) = /(*)+$(x-0<P(0d' 



Solution. Here, K(x, t) = x—t; 1=1; hence, by (14), 
a x (x) = 1 , and all the other a k (x) = 0. 

In this case, equation (15) has the form 

^y^-8(x. t; 1) = 

whence 

g(x, t; \)=g(x, t) = C 1 (t)e' + C t (t)e-* 



(21) 



Conditions (16) yield 

\ C^Oe'-CJOtr'^l 
Solving the system (21), we find 

and, consequently, 

g(x, = -§-(^ _t — e_u ~") = 8inh(jc— 

According to (17) 

R(x, t; l) = [sinh(jc — /)]i = sinh(;t — t) 

Find the resolvent kernels of integral equations with the 
following kernels (X = 1): 

26. K(x, t) = 2—(x—t). 

27. K(x, t) = ~2 + 3(x—t). 

28. K(x, t) = 2x. 
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30. Suppose we have a Volterra-type integral equation, 
the kernel of which is dependent solely on the difference 
of the arguments: 

X 

<p(x) = f(x)+\K(x—t)q>(t)dt (K=l) (22) 
o 

Show that for equation (22) all iterated kernels and the 
resolvent kernel are also dependent solely on the difference 
x— t. 

Let the functions f(x) and K(x) in (22) be original 
functions. Taking the Laplace transform of both sides of 
(22) and employing the product theorem (transform of a 
convolution), we get 

<D(p) = F(p) + K(p)<D(p) 

where 

q>(x) = a>(p), 
f(x) = F(p), 

K(x) = K( P ) 

Whence 

<D(p)= f i p) , K(p)^l (23) 

Taking advantage of the results of Problem 30, we can 
write the solution of the integral equation (22) in the form 

X 

<f(x) = f(x)+^R(x-t)f(t)dt (24) 
o 

where R (x — t) is the resolvent kernel for the integral 
equation (22). 

Taking the Laplace transform of both sides of equation (24), 
we obtain 

0>(p) = F(p) + R(p)F(p) 

where 

R(x)=±R(p) 

Whence 

W^^^fep 1 (25) 
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Substituting into (25) the expression for 0>(p) from (23), 
we obtain 

ft(p)= R ! p) (26) 
\-K(p) 

The original function of R(p) will be the resolvent kernel 
of the integral equation (22). 

Example. Find the resolvent kernel for a Volterra integral 
equation with kernel K(x, t) = s'm(x — /), K=l. 

Solution. We have R (p) = . By (26) 

l 



P 2 +l 



Hence, the required resolvent kernel for the integral equa- 
tion is 

R(x, t; \) = x—t 

Find the resolvent kernels for Volterra-type integral equa- 
tions with the kernels (k=l): 

31. K(x, /) = sinh(x— t). 

32. K(x, = *- u_o . 

33. K(x, = e" ( *" <) sin(A;— t). 

34. K( X , t) = cosh(x—t). 

35. K(x, /) = 2cos(x—0- 

Example. With the aid of the resolvent kernel, find the 

solution of the integral equation 



cp(x) = e* 2 +$e* 2 - l2 (p(/)d/ 

o 



Solution. The resolvent kernel of the kernel K (x, t) = 
=e* 2 -' 2 for X=\ is R(x, t; \) = e x ~ t e" 1 -* 2 (see No. 21). 
By formula (13), the solution of the given integral equation 
is 



cp ( X ) = e* 2 + \ e*-' e* 2 ~' 2 e { ' dt 
o 
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Using the results of the preceding examples, find (by 
means of resolvent kernels) solutions of the following inte- 
gral equations: 

X 

36. (f(x) = e x +le x -\(t)dt. 

o 

x 

37. <j> (x) = sin x + 2$ e* - ' <p (t)dt. 

o 

x 

38. q>(x) = x3 x — $3*-'q>(0d/. 

o 

30. <p W = e * s inx + j§±^<p(0d/. 
o 

40. q>(x) = 1— 2x— J q>(/)d/. 

o 

41. <p(x) = e* 1+2x + 2 Je* 2 -<*(p(/)d*. 

o 

42. <p(x)=l+tf + jl±^q>(0ttf. 

o 

x 

43. (p (x) = — l — + f sin (x— t) q> (/)d*. 
1 +* J 



o 



44. q> (x) = xe 2 + J e^*"" <p (/) d*. 

o 

45. (f (x) = e~* + J e -1 * - " sin (jc — <p (/) d*. 



Note 1. The unique solvability of Volterra-type integral 
equations of the second kind 

X 

<p(x) = f(x) + %\K(x, t)y(t)dt (1) 



VOLTERRA INTEGRAL EQUATIONS 



29 



holds under considerably more general assumptions with 
respect to the function f(x) and the kernel K(x, t) than 
their continuity. 

Theorem. The Volterra integral equation of the second 
kind (1), whose kernel K(x, t) and function f(x) belong, 
respectively, to spaces L, (£2 ) and L 2 (0, a), has one and 
only one solution in the space L 2 (0, a). 

This solution is given by the formula 

X 

<p(x)=f(x) + X$R(x, t; X)f(t)dt (2) 



where the resolvent kernel R(x, t; X) is determined by 
means of the series 

00 

R {x, t; X) = 2 (x, t) (3) 

v=o 

which is made up of iterated kernels and converges almost 
everywhere. 

Note 2. In questions of uniqueness of solution of an 
integral equation, an essential role is played by the class 
of functions in which the solution is sought (the class of 
summable, quadratically summable, continuous, etc., func- 
tions). 

Thus, if the kernel K(x, t) of a Volterra equation is 
bounded when x varies in some finite interval (a, b) so that 

\K(x, OKM, M = const, x g (a, b) 

and the constant term of f(x) is summable in the interval 
(a, b), then the Volterra equation has, for any value of X, 
a unique summable solution (p(x) in the interval (a, b). 

However, if we give up the requirement of summability 
of the solution, then the uniqueness theorem ceases to hold 
in the sense that the equation can have nonsummable so- 
lutions along with summable solutions. 

P. S. Uryson ([29]) constructed elegant examples of 
integral equations (see Examples 1 and 2 below) which have 
summable and nonsummable solutions even when the kernel 
K (x, t) and the function / (x) are continuous. 
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For simplicity we consider f(x) = and examine the 
integral equation 

<p(x) t)<e(t)dt (1) 



where K(x, t) is a continuous function. 
The only summable solution of equation (1) is <p (%)==(). 
Example 1. Let 



K(x, t) = { 



x, 
\ 0, 



t > x 



(2) 



The kernel K (x, t) is bounded in the square Q„ {0 < x, 
*<1}, since 0</C(x, *• What is more, it is 

continuous for O^/^x. In this case, equation (1) has 
an obviously summable solution <p(x) = Q, and, by virtue 
of what has been said, this equation does not have any 
other summable solutions. 

On the other hand, direct verification convinces us that 
equation (1) has an infinity of nonsummable solutions in 
(0, 1) in the form 

<P(*)= — 



where C is an arbitrary constant and x^O. 

Indeed, taking into account expression (2) for the kernel 
K (x, t), we find 



J/C(x, t)<p(t)dt= J te** 1 %~dt + 

o o 

+ J x^-dt=Cx+Cx\ne Tl ~ l =-^ 



i 

xe 
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I 

Thus we obtain 

X X v ^ ' 

Q 

This means that <p(*)=— is a nonsummable solution of 
equation (1). 

Example 2. Let 0^.i^x<a (a>0, in particular a= + oo), 

The function /C(->c, is even holomorphic everywhere, 
except at the point (0, 0). However, equation (1) with ker- 
nel (3) admits nonsummable solutions. Indeed, the equation 

X 

, . . 2 C xt 2 , n\ Jt 2 arctan* 2 ... 



has a summable solution since the function 

2 arctanx 2 



ji 



is bounded and continuous everywhere except at the point 
x = 0. 
The function 

C 0, x = 0, 

where ^(x) is a solution of (4) will now be a nonsummable 
solution of (1) with kernel (3). 
Indeed, for x > we have 

X XX 

\K(x, /MOd-fJ^tWctf + U^ (6) 



By virtue of equation (4), the first term on the right of (6) is 

. , , . 2 arctanx 2 
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The second term yields 

2- [ -rr-ji = — f~3 arctan ^] I'** = arctan (x > 0) 

n J x*+ / s n V* a * J \t=o nx* x 2 v ^ ' 

o 

Thus 

§K(x, t)(p(t)dt = Mp(x) + t arct *" - + ^ arctan ^ 

+ ^ = <P(*) 

which means that the function y(x) defined by (5) is a 
nonsummable solution of equation (1) with kernel (3). 
Example 3. The equation 

X 

(p(x)= i \ ) t x - t (p(t)dt (0<x, 



has a unique continuous solution (p(x)^0. By direct sub- 
stitution we see that this equation also has an infinity of 
discontinuous solutions of the form 

(f>(x)^Cx x ~ 1 
where C is an arbitrary constant. 

4. The Method of Successive Approximations 

Suppose we have a Volterra-type integral equation of the 
second kind: 

X 

+ $K(x, t)y(t)dt (1) 



We assume that f{x) is continuous in [0, a] and the kernel 
K (x, t) is continuous for ^ x ^ a, t ^ x. . 

Take some function <f> t (x) continuous in [0, c]. Putting 
the function <j> (*) into tire right side of (1) in place of 
(p(x), we get 

X 

q) 1 W = /(x) + X \K{x, t)%(t)dt 
o 
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The thus defined function (Pi(x) is also continuous in the 
interval [0, a]. Continuing the process, we obtain a sequence 
of functions 

<Po(*). <PiM <P»W. ••• 

where 

X 

<P„ (*)=/(*) + a. \K{x, 0<P„-i(0# 

o 

Under the assumptions with respect to / (x) and K {x, t), 
the sequence {<p„(#)} converges, as n — >oo, to the solution 
<p(x) of the integral equation (1) (see [13]). 

In particular, if for (p (x) we take f(x), then <p„(x) will 
be the partial sums of the series (2), of Sec. 3, which 
defines the solution of the integral equation (1). A suitable 
choice of the "zero" approximation <p (x) can lead to a rapid 
convergence of the sequence {<p„(x)} to the solution of the 
integral equation. 

Example. Using the method of successive approximations, 
solve the integral equation 

X 

q>(x) = \ + l<f>(t)dt 



taking q> (x) = 0. 
Solution. Since y (x)=Q, it follows that (p t (A:) = 1. Then 

X 

q>,(je)=l + S \-dt = \-\-x, 



X 

<Pa (*)= l+J.O + t)dt~l + X + ?, 
a 

x 

<p 4 (x)=l+j(l+* + y) dt=l+x + §+^ 
o 

Obviously 

<p nl x)=l +_+-+. ,. +? _ - 

Thus, cp„ (x) is the nth partial sum of " the series 

CO 

E — = e x . Whence it follows that <jp„ (x) — <- e x . It is easy 



n = o 
2-1583 
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to verify that the function (p(x) = e x is a solution of the 
given integral equation. 

Using the method of successive approximations, solve the 
following integral equations: 

X 

46. (p(x) = x — \(x— t)<p(t)dt, q> (x) = 0. 



47. <p(x)=l — $(jc— t)(f>(t)dt, <p„ (x) = 0. 



X 

48. <p(x) = 1 + l(x— t)<p(t)dt, (p t (x) = l. 



X 

49. <p(x)=x+ 1 — l(p(t)dt; 



(a) y (x) = l, (b) tp (x) = *+l. 

50. (p(x) = -j+Jt — V <p (/) d/; 

o 

(a) q> (x)=l, (b) <p (x) = x, (c) <p (x) = y + ;c. 

x 

51. q>(x)=l+X + t)(f(t)dt, (p (*) = l. 



52. (p(x) = 2x + 2— J <p(/)d*; 



(a) <p (x) = l, (b) <p,(x) = 2. 

53. <p(x) = 2x a + 2 — $x<p 

(a) «p,(x) = 2, (b) q>.(*) = 2x- 

X 

54. q>(x) = ^— 2x— jq>(<)df. (p (x) = X a . 

o 

55. Let K (x, t) satisfy the condition 

a x 
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Prove that the equation 

X 

<p(x)— k J K(x, t)<p(t)dt=0 



has, for any A,, a unique solution <p(x) = in the class 
MO, a). 

The method of successive approximations can also be 
applied to the solution of nonlinear Volterra integral equa- 
tions of the form 

X 

y(x) = y 9 +$F[t, y(t))dt (2) 



or the more general equations 

X 

<p(x) = /(*)+ t, <p(t))dt (3) 



under extremely broad assumptions with respect to the func- 
tions F (x, t, z) and f(x). The problem of solving the diffe- 
rential equation 

% = F(x, y), y\ x=0 = y<> 

reduces to an equation of the type (2). As in the case of 
linear integral equations, we shall seek the solution of equa- 
tion (3) as the limit of the sequence {<f„(x)} where, for 
example, <p (jc) = / (x) , and the following elements y k (x) are 
computed successively from the formula 

X 

^ k (x) = f(x) + ^F(x, t, y k _ t {t))dt (ft=l, 2, ...) (4) 
i> 

If f (x) and F(x, t, z) are quadratically summable and sa- 
tisfy the conditions 

\F{x, t, z 2 )—F(x, t, z t )\<a{x, t)\z t —z t \ (5) 
\F(x, t,f(t))dt (6) 

2* o 
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\ 

where the functions a(x, t) and n(x) are such that in the 
main domain (0 ^ t ^ x ^ a) 

a ax 

J n a AT 4 , Jd*Ja»(*, t)dt^A* - (7) 

o o o 

it follows that the nonlinear Volterra integral equation of 
the second kind (3) has a unique solution q>(x) €^(0, a) 
which is defined as the limit of q>„(x) as n — ► <»: 

<p (x) = lim q>„ (x) 



fn ' 



where the functions <p„ (x) are found from the recursion for- 
mulas (4). For % (x) we can take any function in L 2 (0, a) 
(in particular, a continuous function), for which the condi- 
tion (6) is fulfilled. Note that an apt choice of the zero 
approximation can facilitate solution of the integral equation. 

Example. Using the method of successive approximations, 
solve the integral equation 

„.(*)-|!±$2 dt 



taking as the zero approximation: (1) (p (x) — 0, (2) <p (x) = x. 
Solution. (1) Let %(x) = 0. Then 

<PiW = J TT72 = arctan x ' 


X 

<P a (x) = J ^'rctan'i ^ _ arc { an ^ _j_ i. arctan 3 a;, 



* 1 + f arctan ' + arctan 3 / 
<p s (X) = L ^ y+T* — dt = arctan x + 



o 



i 2 1 

+ y arctan 3 x + -^g- arctan 6 jc + arctan 7 x, 

X 

cp 4 (x) = J t = arctan x + -g- arctan 3 x + 
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+ 3X5 arCta " 5 X + 5775<9 arCtan? X + 5X7X& X + 

134 4 

r arctan 11 x + „- arctan 13 x + 



^9X 11x21x25 3x5x7x9x13' 

+ 7 2 X9 2 X i5 arctan 15 x, .... 

Denoting arctan x--u and comparing expressions for <p n (x) 
with the expansion 

, A 2 2 »(2 2 * — 1) ji 

tan« = £(_i)v-i ( ^ 2v)| B 2v u*~\ |u|< T 

V=l 

where B v are Bernoulli numbers,* we observe that 
cp„ (x) — * tan (arctan a:) = x 

It can easily be verified that the function (p(x) = x is a so- 
lution of the given integral equation. 
(2) Let <p (x) = x. Then 



f 1 -ht 2 

<Pi (*) =J T+t5 dt==x 



In similar fashion we find <p„ (x) = x (n = 2, 3, ...). 

Thus, the sequence {cp„(A:)} is a stationary sequence {x}, 
the limit of which is <p(x) = x. The solution of this integ- 
ral equation is obtained directly: 

cp (x) = x 

56. Use the method of successive approximations to solve 
the integral equation 



* The Bernoulli numbers S 2v+1 with odd index (odd degree) are all 
equal to zero, except B x = — ■ The number B =l, the numbers B 2v 
are obtained from the recursion formulas 

2V- 2 

B -- 1 . 1 V 2v(2v-l)...-(2v-2ft + 2) fl 
Bu ~ 2vTI + 2— 2- 51 

t = 2 
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57. Use the method of successive approximations to find 
a second approximation <p 2 (x) to the solution of the integ- 
ral equation 

q>(x)=l + JVW + fcp (t) + t*]dt 



58. Use the method of successive approximations to find 
a third approximation cp 3 (x) to the solution of the integral 
equation 

X 

q>(x)= J[*cp 2 (/)— \)dt 
o 



5. Convolution-Type Equations 

Let and <p 2 (#) be two continuous functions defined 

for x^O. The convolution of these two functions is the fun- 
ction cp 3 (x) defined by the equation 

X 

<M*)= \ <M*— t)^(t)dt (1) 



This function, defined for x^O, will also be a continuous 
function. If q> 1 (x) and <p 2 (x) are original functions for the 
Laplace transformation, then 

J?q> 1 = J?q»i-J2'q), (2) 

i. e., the transform of a convolution is equal to the pro- 
duct of the transforms of the functions (convolution theorem). 

Let us consider the Volterra-type integral equation of the 
second kind 

q> (x) = /(*)+ J /C(x-0«P(<)# (3) 
o 

the kernel of which is dependent solely on the difference 
x — t. We shall call equation (3) an integral equation of the 
convolution type. 
Let / (x) and K (x) be sufficiently smooth functions which, 
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as x — »-oo, do not grow faster than the exponential function, 
so that 

l/WKM/.', I K (x) |< M 2 e s ' x (4) 

Applying the method of successive approximations, we can 
show that in this case function <p(x) will also satisfy an 
upper bound of type (4): 

|<p(*)|< M 3 eHX 

Consequently, the Laplace transform of the functions f(x), 
K(x) and (p(x) can be found (it will be defined in the half- 
plane Re p = s>max(s 1 , s 2 , s a )). 

Let 

f(x) = F(p), <p(x) = a>(p), K(x) = K(p) 

Taking the Laplace transform of both sides of (3) and emp- 
loying the convolution theorem, we find 

<S>(p) = F(p) + K(p)<S>(p) (5) 

Whence 

0(p)= F( f ) (K(p)=£l) 
\—K(p) 

The original function q>(x) for O(p) will be a solution 
of the integral equation (3) (see [25]). 
Example. Solve the integral equation 

X 

(p (x) = sin x + 2 Jcos(x — t)q>(t)dt 

o 

Solution. It is known that 

sin^^qrp cosx = ^ 

Let (f(x)z=z(D(p). Taking the Laplace transform of both 
sides of the equation and taking account of the convolution 
theorem (transform of a convolution), we get 
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Whence 



+ 1 



or 



Hence, the solution of the given integral equation is 

<p (x) = xe x 

Solve the following integral equations: 

59. (p(x) = e x — ^e*- l (p(t)dt. 



X 

60. q> (*) = *— ^ e*~ \(t) dt. 



61. (p(x) = e 2x ^.^e t - x (p(t)dt. 



at 

62. <p (*) = *— t)y(t)dt. 



63. (p(x) = cosx— ^(*— *)cos(a:— t)cp(t)dt. 



64. <p(x) = l+x + Je" 2U -»(p(/)d/. 



a: 

65. cp(x) = A:+ ^ sin(x— /) q> (/)d/. 



66. cp(x) = sin*+^(.*;— t)<p(t)dt. 



67. <p(x) = x— £ sinh(x— t)y(t)dt. 
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68. 



<p(x)= 1— 2x— 4x 2 + J [3—Q(x—t) — 



—4(x — t)*]<p(t)dt. 

X 

69. cp(A:) = sinh x — ^ cosh(x — t)q>(t)dt. 



X 

70. 9 (x) = 1 + 2 J cos (x—t) <p (0 df. 



71. <p(x) = e* + 2 Jcos(x— t)y(t)dt. 



X 

72. cp(x) = cosa;+ ^<p(t)dt. 



The Laplace transformation may be employed in the solu- 
tion of systems of Volterra integral equations of the type 

<P, <*)=//(*) +2 {K u (x-t)v,(t)dt (t = l, 2, .... s) (6) 

3 = 1 „ 

where Kij(,x), f t (x) are known continuous functions having 
Laplace transforms. 
Taking the Laplace transform of both sides of (6), weget 

® i (p) = F i (p)+ J R u (p)®,(p) (t = l, 2, .... s) (7) 

This is a system of linear algebraic equations in <D y (p). 
Solving it, we find d> / (p), the original functions of which 
will be the solution of the original system of integral 
equations (6). 
Example. Solve the system of integral equations 



q>, (X) = 1 - 2 J e* «-'> q> 1 (<) d/ + $ 9, (0 



<p r (x) = 4x — J q> t (/) d/ + 4 J (x q>, (/.) dt 



(8) 
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Solution. Taking transforms and using the theorem on 
the transform of a convolution, we get 

®i (P) = 7 - p=2 ®i M + 7 °* W> 

Solving the system obtained for ^(p) and ® 2 (p), we find 

pi 1 



<D,(J>) = 



<D,(P) = 

3p+2 



(P+l) a P+l (P + l) 2 
8 1,1 1 



8 



1 



(p_2)(p+l)« 9 p-2 1 3 (p+l)* 9 p+l 

The original functions for <X>j (p) and <D 2 (/>) are equal, 
respectively, to 

q> t (x) = e~ x — xe~ x , 

The functions (f t (x), <p 2 (x) are solutions of the original 
system of integral equations (8). 
Solve the following systems of integral equations: 



73. 



74. 



75. 



q> t (x) = sinx+ J(p 2 (0 dt, 


X 

<p 2 (x) = 1 — cosx — 5 <Pi (0 ^* 

- 



q>, (x) = l-Se ,u -"q»i (0*- 


<p, (x) = e* + J «p, (0 d/ — J q>, (/) dt, 


q>, (jc) = — x— I (x—t) <p x (t) dt + $q> 2 (/) d/. 
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76. 



77. 



78. 



79. 



<Pl (x) = 


= «*— S<Pi (0# + 4$e*-*<P.(0#, 






* AT 




<M*) = 





<Pi w = 


Jt 

1 — Jq> x (t)di, 






* 

sin x + -g- J (x — t) cp x (<) dt. 

o 


<Pi (*) = 


o 


<Pi w - 


X 

COSX — 1 -+- 5 •PsCOdf, 




<Ps (*) = 


Jt 

cos x -\- J cp t (<) d<. 



<Pi (*) = 


< 

* + 1 + S (Pa (0*. 



(p 2 (X) = 


— x + 5 (* — <Pi if) dt, 




<Ps (*) = 


cos x — 1 — 5 *Pi (0 ^* 





6. Solution of Integro-DifFerential Equations 
with the Aid of the Laplace Transformation 



A linear integro-differential equation is an equation of 
the form 
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a (x) cp" {x) + aj (x) qp" -1 (x) + . . . + a n (x) cp (x) + 

+ 2 )K m (x, t)^(t)dt = f(x) (1) 

m=0 o 

Here a„(x), a n (x), f(x), K m (x, t) (m = 0, 1, . .., s) 
are known functions and <p(jc) is the unknown function. 

Unlike the case of integral equations, when solving in- 
tegro-differential equations (1), initial conditions of the form 

<P(0) = q>„, <p'(0) = (p' o <p»-i(0) = <pr u (2) 

are imposed on the unknown function <p(x). In (1), let the 

coefficients a k (x) = const (k = 0, 1 n) and let K m (x, t) = 

= K m (x — t) (m = 0, 1, .... s), that is, all the K, m depend 
solely on the difference x — t of arguments. Without loss 
of generality, we can take a = l. Then equation (1) assu- 
mes the form 

(p n (x) + a^-" (x) + . . . + a„9 (x) + 
+ 2 )K m (x-t)<f>™(t)dt=f(x) (a u a„ -const) (3) 

m=0 o 

Also, let the functions f(x) and be original func- 

tions and 

f(x) = F(p), K m (x) = R m (p) (m = 0, 1, .... s) 
Then the function cp (x) will also have the Laplace transform 

<p(x)=a>(p) 

Take the Laplace transform of both sides of (3). By virtue 
of the theorem on the transform of a derivative, 

(p l *> (x) = p*0(p)— p*" 1 q> — p*" 2 9; — . . . — «pf*- u (4) 

(k = 0, 1 «)■ 

By the convolution theorem 

X 

\K m {x-t)^{t)dt = 


rR a {p) [p m ® (p)-p"- 1 9o- ■ • • -q>r _1) ] (5) 
(m = 0, 1 s) 
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Equation (3) will therefore become 
(D(p) 



p" + a lP "- 1 + ...+a n + S R(p)p m 

m=0 
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A(p) (6) 



where A(p) is some known function of p. 

From (6) we find O(p), which is an operator solution 
of the problem (3)-(2). Finding the original function for 
<D(p), we get the solution <p(x) of the integro-differential 
equation (3) that satisfies the initial conditions (2). 

Example. Solve the integro-differential equation 

X 

q>" (x) + J e 2 u " * > q>' (/) dt = e %x (7) 

cp(0) = <p'(0) = (8) 
Solution. Let (p(x) = <D(p). By virtue of (8) 

<p'(x)= P a>( P ), 
<p"(x)= P *a>(p) 

Therefore, after taking the Laplace transform, equation (7) 
will assume the form 

pHD(p) + -i_0(p )a .. 1 



P-2 



or 



«D(p) 



p (p _ 1)2 _ 1 

p-2 p-2 



(9) 
(10) 



From (10) we find 

<D(p) = 



l 



p(p-1) : 



Hence, the solution cp(^) of the integro-differential equation 
(7) satisfying the initial conditions (8) is defined by the 
equality 

y(x) = xe x — e*-f 1 
Solve the following integro-differential equations: 

80, <p"(x)+ J e 2 «-'>q>' (t)dt = e lx ; «p(0) = 0, cp'(0)=l. 
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x x 

81. (p'(x)— <p(x)+ $(x— t)<p'(t)dt- $<p(/)d/ = x; 



<p(0) = - 1. 

X 

82. cp" (x) — 2q>' (x) + <p (x) + 2 J cos (x — q>* (t) dt + 

o 

x 

+2 J sin {x—t) «p' (/) d< = cos x; <p (0) = <p' (0) = 0. 



X 

83. (p"(x) + 2q>'(x) — 2 $ sin(x— t)<p'(t)di = cosx; 

o 

(p(0) = cp'(0)=0. 

84. <p"(x) + «p(x)+J sinh(x— /)<p(0# + 

o 

x 

+ $ cosh (x—t) <p' (/) d*=cosh x; <p (0)=<p' (0)=0. 

n 

85. <p"(x) + q>(x)+$ s!nh(x— /)<p(f)<tt + 

o 

+ J cosh(x— <p' (0 = cosh x; q>(0) = — 1, 

o 

<P'(0)= 1. 

7. Volterra Integral Equations with Limits (x, + oo) 

Integral equations of the form 

<p(xW(*)+$/C(x — <)<P(*)* (1) 

which arise in a number of problems in physics can also 
be solved by means of the Laplace transformation. For 



(3) 
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this purpose, we establish the convolution theorem for the 
expressions 

00 

$*(* — <)q>(0# (2) 
It is known that for the Fourier transformation 

I S S(x— *)ii>(0<ftj =V2nG(k)W(%) 

where G(k), W (k) are Fourier transforms of the functions 
g(x) and ty(x), respectively. 
Put g(x) = K. (x), i.e., 

( 0, x>0, 
( (p (x), x>0, 

i|>(*)=q>+(*) = | o, *<o < 4 ) 

Then (3) can be rewritten as 

r j J /C(x—0*(0*} = K^ (5) 

(here and henceforward the subscripts <F or will mean 
that the Fourier transform or the Laplace transform of the 
function is taken). 

To pass from the Fourier transform to the Laplace trans- 
form, observe that 

Fje(p) = V2H[f + (ip)] 3! (6) 

Hence, from (5) and (6) we get 

^\lK(x-t)if(t)dtj = V2a[K.{ip)),[0 + (p)]jg (7) 

We now express [V^n R_ (ip)\f in terms of the Laplace 
transform: 

00 

[V'2n R_(ip)] f = J K(x)e-P*dx= \K(— x)eP*dx 



4 g INTEGRAL EQUATIONS 

Putting K{—x) = 9C(x),..we get 

[V%i K_ (ip)\, = Xjg (- p) = J K (— x) dx 



And so 

j? | J /C ( x - 1) <p (0 d/ 1 = ^ ( - p) d> x (P) (8) 

Let us now return to the integral equation (1). Taking 
the Laplace transform of both sides of (1), we obtain 

®(p) = F(p) + K(-p)®(p) (9) 

(dropping the subscript £) or 

®(P)= ; i ( f : (W(-P)^l) (10) 

where 

00 

»( — p) = $ K(— (11) 



The function 

*W = i1 T=^ M " (I2) 

V-'oo 

is a particular solution of the integral equation (I). It must 
be stressed that the solution (9) or (12) is meaningful only 
if the domains of analyticity of 9C( — p) and F (p) overlap 
(see [17]). 
Example. Solve the integral equation 

40 

<p(x) = x+ $«»«*-*> <p(f)# (13) 
Solution. In this case, f(x) = x, K(x) = e ix . Therefore 

00 

F(P)^^-. W(-p) = §e-* x eP*dx = ^> Rep < 2 

o 
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Thus, .we obtain the following operator equation: 

so that 
Whence 

■y + i'oo 

VW-SZ J 1^=1) e " d P (°<V<2) (15) 

V— too 

Integral (15) may be evaluated from the Cauchy integral 
formula. The integrand function has a double pole p — 
and a simple pole p=l, which appears for v> 1; this is 
connected with including or not including in the solution 
of equation (13) the solution of the corresponding homoge- 
neous equation 

00 

<p(jc) = Je 2 u -'><p(/)ci/ 

X 

Let us find the residues of the integrand function at its 
poles: 

res ( /~ 2 ,, eP*) = 2x+ 1, res [ f,~ 2 u eP*) = — e x 

Consequently, the solution of the integral equation (13) is 
y(x) = 2x + I + Ce x (C is an arbitrary constant). 
Solve the integral equations: 

00 

86. <p(x) = e~*+ \<p(t)dt, 

* 
CO 

87. <p(x) = e-* + ^e*- t (p(t)dt. 

x 

00 

88. <p (*) = cos* e*-'<p 

X 

CO 

89. <p(jc)=l + $e B «*--* > q>(Odf (<x>0). 
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8. Volterra Integral Equations of the First Kind 

Suppose we have a Volterra integral equation of the first 
kind 

^K(x, t)<p(t)dt = f(x), /(0) = (1) 
o 

where q>(x) is the unknown function. 

Suppose that K(x, t), dK g X x ' - , /(*) and f (x) are conti- 
nuous for 0^je<;a, O^fs^x. Differentiating both sides 
of (1) with respect to x, we obtain 

X 

K (x, x) cp (x) + j 2*|ii> cp (t)dt = r (x) ■ (2) 



Any continuous solution cp(x) of equation (1), for a, 
obviously satisfies equation (2) as well. Conversely, any 
continuous solution of equation (2), for O^jc^a, satisfies 
equation (1) too. 

If K(x, x) does not vanish at any point of the basic 
interval [0, a], then equation (2) can be rewritten as 

o 

which means it reduces to a Volterra-type integral equation 
of the second kind which has already been considered 
(see [18]). 

If K(x, x) = 0, then it is sometimes useful to differen- 
tiate (2) once again with respect to x and so on. 

Note. If K(x, x) vanishes at some point x£[0, a], say 
at x — 0, then equation (3) takes on peculiar properties that 
are quite different from the properties of equations of the 
second kind. (Picard called them equations of the third 
kind.) Complications arise here similar to those associated 
with the vanishing of the coefficient of the highest deriva- 
tive in a linear differential equation. 

Example. Solve the integral equation 

X 

lcos(x—t)<p(t)dt = x (4) 
o 
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Solution. The functions f(x) = x, K(x, t) = cos(x — t) 
satisfy the above-formulated conditions of continuity and 
differentiability. 

Differentiating both sides of (4) with respect to x, we get 

X 

<p(x) cosO — ^ sin (x — /) q> (/) dt = 1 



or 

X 

(p(x) = l + Jsinfx— t)q>(t)dt (5) 
o 

Equation (5) is an integral equation of the second kind of 
the convolution type. 

We find its solution by applying the Laplace transfor- 
mation 

®(P) = J + ^®(P) 

whence 

<1) (P) = -^- = 7 + ^^= 1 +T 

The function cp(x)=l+-^- will be a solution of equation 

(5) and hence of the original equation (4) as well. This is 
readily seen by direct verification. 

Solve the following integral equations of the first kind 
by first reducing them to integral equations of the second 
kind: 

X 

90. $e*-'<p(0d' = sin*. 



91. J 3*~'(p (/)<#=*. 

o 

x 

92. $ a*"' <p (')<# = /(*), /(0) = 0. 

o 

x 

93. J(l-* , + /*)<p(0<tt— f • 
o 
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x 

94. $(2 + x 2 — t 2 )<p(t)dt = x\ 
o 

x 

95. J sin(x— t)(p(t)dt = e*'/*— 1. 
o 

9. Euler Integrals 

The gamma function, or Euler's integral of the second 
kind, is the function T (x) defined by the equality 

00 

r (x) = \e- t t x ~ 1 dt (1) 

o 

where x is any complex number, Re x > 0. For x = 1 we 
get 

r (1)= J = 1 (2) 

o 

Integrating by parts, we obtain from (1) 

OS 

Y{x) = ±le-n*dt=^±V (3) 
o 

This equation expresses the basic property of a gamma 
function: 

T{x+\) = xT(x) (4) 

Using (2), we get 

r(2) = r(i + i) = i.r(i) = i, 
r(3) = r(2+i) = 2-r(2) = 2i, 
r(4) = r(3+ i) = 3-r (3> = 3i 

and, generally, for positive integral n 

I» = (n-1)! (5) 

We know that 
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Putting x = t 2 here, we obtain 




o 



Taking into account expression (1) for the gamma function, 
we can write this equation as 



Whence, by means of the basic property of a gamma 
function expressed by (4), we find 



Generally, it will readily be seen that the following equa- 
lity holds: 



(n a positive integer). 

Knowing the value of the gamma function for some value 
of the argument, we can compute, from (3), the value of 
the function for an argument diminished by unity. For 
example, 






(6) 




For this reason 




(7) 



2 



Acting in similar 



fashion, we find 




2 
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2 

r ^ — t) = — tb"^" anc ' so on ' 

It is easy to verify that r(0) = T(— 1)= . . . =T(— n) = 
= ...=oo. Above we defined r (x) for Rex > 0. The indi- 
cated method for computing T (jc) extends this function into 
the left half-plane, where T (x) is defined everywhere except 
at the points x — — n (n a positive integer and 0). 
Note also the following relations: 

r(x)T(l—x)=^—, (8) 

v ' v ' sin n x ' v ' 

r(x)T(x + ^=2 1 -**n l />r(2x) (9) 

and generally 

r W r(,+±)r(,+±)...r (,+!=!)- 

= (2it) 2 rt 2 ~"*r(/u) 

(Gauss-Legendre multiplication theorem). 

The gamma function was represented by Weierstrass by 
means of the equation 



where 



{(!+*)•-*} (l0) 



7= lta (l+| + i.+ ... + l._l„ m )_0. 

m -* od * / 



57721. 



is Euler's constant. From (10) it is evident that the func- 
tion T(z) is analytic everywhere except at z = Q. z — — 1, 
2 = — 2, .... where it has simple poles. 

The following is Euler's formula which is obtained 
from (10): 

r «-ra{( l+ x)'('+T)"'} <"> 
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It holds everywhere except at 2 = 0, z = — 1, z — — 2 

96. Show that T'(l)=—y. 

97. Show that for Rez>0 



98. Show that 

r'(D f '(t) { 
r(i) r ^ z inz 

99. Prove that 

T (z) = lim lx2 --- (w - 1) n* 

We introduce Euler's integral of the first kind B(p, q), 
the so-called beta function: 
i 

B(p, q) = ^xP- 1 (l-x)o- 1 dx (Rep>0, Re<7>0)(l2) 



The following equality holds (it establishes a relationship 
between the Euler integrals of the first and second kinds): 

«(/>. — ro , +9) (13) 

100. Show that 

B(p, q) = B(g, p) 

101. Show that 

B(p, q) = B{p+l,q) + B(p,q+\) 

102. Show that 

B(p+l,q) = fB(p, q+l) 

103. Show that 

i 

jj ( 1 + xy- 1 ( 1 — x)i- 1 dx = 2P+1- 1 B (p, q) 
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104. Evaluate the integral 

n 

2 

/ = J cos m - 1 xsm"- 1 xdx (Rem>0, Ren>0) 



10. Abel's Problem. 
Abel's Integral Equation and Its Generalizations 

A particle is constrained to move under the force of 
gravity in a vertical plane (%, r\) along a certain path. It is 
required to determine this path so that the particle, having 
started from rest at a point on the curve (path) with ordi- 
nate x, reaches the £ axis in time t — f 1 (x), where ■f l (x) is 
a given function (Fig. 1). 




Fig. I 



The absolute velocity of a moving particle isv=y r 2g(x — r\). 
Denote by p the angle of inclination of the tangent to the 
£-axis. Then we will have 



§=-K2g(x-T!)sinp 



whence 



dt = 



V2g(x-x\) sinf( 



Integrating from to x and denoting —i—g- = <p 0l). we get 

sin p 

Abel's equation 



mil- VTelM 

V X — T] 
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Denoting —V^gfiix) by f(x), we finally obtain 

f-^L-dWw (i) 

J V x—y\ 



where q>(x) is the required function and f(x) is the given 
function. After finding q>(r\) we can form the equation of 
the curve. Indeed, 

whence 

ri = 0>(P) 

Further -7 

rfr) _ (D' (P) rfft 



whence 



^~ tan p tanfJ 



and consequently, the required curve is defined by the pa- 
rametric equations 

T1 = Q(P) f ( ^ 

Thus, Abel's problem reduces to a solution of the integral 
equation 

X 

f(x) = ^K(x, t)<f{t)di 



with given kernel K{x, t), given function f(x) and unknown 
function q>(%); in other words, it reduces to finding a solu- 
tion of the Volterra integral equation of the first kind. 

The following somewhat more general equation is also 
called Abel's equation: 

X 

l<^dt = f(x) - (3) 
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where a is a constant, 0<a< 1 (Abel's generalized equa- 
tion). We will consider that the function f(x) has a conti- 
nuous derivative on some interval [0, a]. Note that for 

a^-i- the kernel of equation (3) is quadratically noninte- 

grable, i.e., it is not an L 2 -function. However, equation (3) 
has a solution which may be found in the following manner. 

Suppose equation (3) has a solution. Replace x by s in 
the equation and multiply both sides of the resulting equa- 

ds 

lity by : and integrate with respect to s from 

(X S) 

to x: 

\ dx f £ig df_.f f & ffc (4) 



Qianging the order of integration on the left, we obtain 

/ 

where 



In the inner integral make the substitution s = t +y(x — t): 

X 1 

C ds C dy n 

J (x— s)>-«(s— <)" J y) 1-a — sin an 

/ 

Then from equation (5) we have 



or 

, . sin an r , , . sin om { C f (s) , \ 

«P(*)=— ^'W = — (7) 
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Thus, the only solution of equation (3) is given by for- 
mula (7), which, via integration by parts, can also be re- 
written in the form 



. , sin an 

<pW=-t- 



This solution has physical meaning only when its absolute 
value is not less than 1 ^since ep(x)= ^ ^ ^ . 

We will show that in the case / (x) = C = const, the so- 
lution of Abel's problem is a cycloid. (The tautochrone 
problem: to find the curve along which a particle moving 
under gravity without friction reaches its lowest position 
in the same time, irrespective of its initial position.) 

In this case a = y • Hence, by formula (8) 

1 C 



<pW = 



And therefore 



whence 



Further 



sin p = 



n y x 



*l = -§-sin 2 p 



, t dn C 2 2 sin P cos P ,„ C 2 ,, , M . .„ 



Finally, 



S = 5(P+T sin2 P)+ C ' 

* = (l*+T sin Sphe- 
re £(1- cos 2p) 



2ji a ' 

(parametric equations of the cycloid). 

105. Show that when f(x) = CV~x the solution of Abel's 
problem will be straight lines. 



60 INTEGRAL EQUATIONS 

Solve the following integral equations: 

106. l^Bw = xn (°< a < 1 >- 


X 

107. \ y' = sin;t. 



,08. {sm^-*. 

J Vx-t 



,09. A 

J 



110. Solve the two-dimensional Abel equation 

ff <p(x,y)dxdy . 

JJ V(y»-y)*-(x a -xy - nx "' y <> ) 

D 

Here, the domain D is a right isosceles triangle with hypo- 
tenuse on the x-axis and vertex at the point (x , y ). 
Consider the integral equation (see [21]) 

K 

\{x — tfy(t)dt = x* (9) 



(X^O, p> — 1 is real), which in a sense is a further 
generalization of the Abel equation (3). 

Multiply both sides of (9) by {z—xY (jx > — 1 ) and in- 
tegrate with respect to x from to z: 

Z X ? 

5 (2 — xf ($(* — 0Pq>(f)*)<Lt=S* x (z— (10) 



Putting x — pz in the integral on the right side of (10), 
we obtain 

Z 1 

^ x x (z — xYdx = z*+v-+ J ^ p x (l — p)ndp = 



= z^+> B (X. -f 1 , p + 1 ) = z^ r %H^ l) 

(X + H+l>X^0) (11) 
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Changing the order of integration on the left side of 
(10), we get 

Z X 

J ( J (z — xY (x— tf <p (t) dt)dx = 

o o 

z z 

= l(l(*—xY(x—tfdx)(f>(t)dt (12) 

o / 

In the inner integral on the right of (12) put 
x=t + p(z — /) 

Then 

2 1 

J (2 — x)r (x—t)$ dx = (z — ty+*+ l jj pf ( 1 — p)v dp = 

* o 
« ( 2 _^+f>+> fl (0 + 1 , ,1+ 1) = r( r P +^^+ 1) (2- *) P+ * +1 

(13) 

Taking into account (11), (12), (13), we obtain from (10) 
r (p+ 1) 



np+n+2) 



j {z -ty+^ (p (o dt = r ^+^ 2) z^ +l (H) 



Choose u. so that u. + |J -|- 1 = n (a nonnegative integer). Then 
from (14) we will have 



r(P + i) 

T(rt+1) 





\(z-tYW)dt = V( r^ +X) ^ 



or 

Ul=lYlvU)dt= T{X + X) 2 x + »-|J (15) 



Differentiating both sides of (15) n+l times with res- 
pect to 2, we obtain 

r(X+l)(X + n-P)(X + /»-P-l)...(X-P) ^b., Mfi . 

'l 2 ' r(p + i)r(X+«-p+i) 

or for X— p + fe^O (fe = 0, 1, .... «) 
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W (Z)= F (*-+!) 2 \-P-i (17} 

This is the solution of integral equation (9). 

Observe that if the quantity X — p" — 1 is equal to a 
negative integer, we get q>(2) = 0. In this case, equation 
(9) does not have a solution in the class of ordinary func- 
tions. Its solution is a generalized function (see p. 67). 

Example. Solve the integral equation 

X 

\ ) (x—t)y(t)dt=x i 



Solution. In the given case, |3 = 1, X = 2. Since X — p + 

+ k^=0 (k = 0, I, 2 n), it follows, from formula (17) 

that 

«pw=iW(n^- 1 - 1 = 2 

Solve the integral equations: 

* JL — 

111. $(*— t) 3 cp(t)dt = x s — x\ 
a 

x j_ 

112. J (x—t) 2 <p(t)dt=nx. 



X 1 

113. J (x— 4 <p(/)d/ = x + je». 



at 

114. *)* <P (0 # = 



a: 

115c yj"(*— 2 <J>(0<# = cosa:— 



11. Volterra Integral Equations of the First 
Kind of the Convolution Type 

An integral equation of the first kind 

X 

lK(x-t)<?(t)dt = f(x) (1) 
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whose kernel K (x, t) is dependent solely on the difference 
x — t of arguments will be called an integral equation of 
the first kind of the convolution type. 

This class of equations includes, for instance, the gene- 
ralized Abel equation. 

Let us consider a problem that leads to a Volterra in- 
tegral equation of the convolution type. 

A shop buys and sells a variety of commodities. It is 
assumed that: 

(1) buying and selling are continuous processes and pur- 
chased goods are put on sale at once; 

(2) the shop acquires each new lot of any type of goods 
in quantities which it can sell in a time interval T, the 
same for all purchases; 

(3) each new lot of goods is sold uniformly over time T. 
The shop initiates the sale of a new batch of goods, 

the total cost of which is unity. It is required to find the 
law q>(t) by which it should make purchases so that the 
cost of goods on hand should be constant. 

Solution. Let the cost of the original goods on hand at 
time t be equal to K (t) where 

\ 0, t>T 

Let us suppose that in the time interval between t and 
x + dx goods are bought amounting to the sum of qp(r)dx. 
This reserve diminishes (due to sales) in such a manner 
that the cost of the remaining goods at time t > t is equal 
to K(t — x)q>(x)dx. Therefore the cost of the unsold part 
of goods acquired via purchases will, at any time t, be 
equal to 

t 

^K(t—x)y(x)dx 



Thus, <p(/) should satisfy the integral equation 

t 

l—K{t)=lK(t-t)(((x)dx 
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We have thus obtained a Volterra integral equation of 
the first kind of the convolution type. 
Let f(x) and K(x) be original functions and let 

f (x) = F(p), K(x) = K(p), q>(*) = Q(p) 

Taking the Laplace transform of both sides of equation (1) 
and utilizing the convolution theorem, we will have 

R(p)Q>(p) = F(p) (2) 

whence 

®(p) = ?M (K(P)^O) (3) 
K KP) 

The original function q>(jc) for the function O(p) defined 
by (3) will be a solution of the integral equation (1). 
Example. Solve the integral equation 

X 

J e*-\ (t) dt = x (4) 



Solution. Taking the Laplace transform of both sides 
of (4), we obtain 

^®(P) = J* (5) 

whence 

yr p 2 p p 2 

The function (p(x)=\— x is a solution of equation (4). 
Solve the integral equations 

X 

116. ^cos(jc — t) cp (t) dt = sin x. 



X 

117. $e*-'q>(/)d/ = sinh*. 



118. J (x— t) 2 <p(t)dt = x 2 . 



X 

119. < j j e ux - t) <p(t)dt = smx. 
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120. J e*"' <p (f) d< = 



121. ^ cos (jc— <p(t)dt = xs\nx. 

o 
x 

122. $ sinh(x— <P (0 ^ = 



123. ^ J (x—t)y(t) dt = sin x. 



124. Jcosh(x— t)y{t)dt = x. 


X 

125. $cos(jt — t)<p(t)dt = x + x 2 . 


X 

126. £(*■— t 2 )<f>(t)dt = j. 



j: 

127. j(* 2 — 4x< +3/ , )«P(«*=5- 



128. j j (a: 3 — 4x/ + 3/ 2 ) 9 (t) dt = x*J t (2 Vx). 



129. §{x—2t)<((t)dt = — £. 



Note. If *) = /C(0)=7&0, then equation (1) defini- 

tely has a solution. In Problem 122 the kernel K(x, t) 
becomes identically zero for t = x, yet the equation has 
a solution. 

As has already been pointed out before, a necessary con- 
dition for the existence of a continuous solution of an 
integral equation of the form 

X 

Ji^?rr'p(o*=/w (6) 


3-1583 
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consists in the function f (x) having continuous derivatives 
up to the nth order inclusive and in all its n — 1 first 
derivatives vanishing for x = 

This "model" equation (6) points to the necessity of 
matching the orders of vanishing of the kernel for t = x and 
of the right side f (x) for x — (the right side must exceed 
the left by at least unity). 

Consider the integral equation 

X 

\{x — t)y(t)dt = x (7) 



Here, f(x) = x, n = 2. Obviously, f (x) has derivatives of 
all orders, but its first derivative /' (x) = 1 0; that is, 
the necessary condition is not fulfilled. 

Taking the Laplace transform of both sides of (7) in 
formal fashion, we get 

whence 

0)(p) = l 

This is the transform of the c-function 6(x). 
Recall that 

6(x)= 1 

where m is an integer ^0. 

Thus, the solution of the integral equation (7) is the 
6-function: 

<p(x) = 6(x) 

This is made clear by direct verification if we take into 
account that the convolution of the fi-function and any 
other smooth function g (x) is defined as 
g (x) * 6 (*) = g (at) 
ft*>{x)*g(x) = g ll "(x) (*= 1, 2, ...) 
Indeed, in our case g (x) = K (x) = x and 

X 

< \ j K(x — l)bU)dt = K {x) = x 
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Thus, the solution of equation (7) exists, but now in 
the class of generalized functions (see [4j, [22\). 
Solve the integral equations: 

X 

130 i \ j {x — t)<p(t)dt = x* + x—l. 

X 

131. l(x— t) <p (I) dt = smx 



X 

132. ^(x — tfy(t)dl = x*-\-x 3 . 



X 

133. ^s\n(x — t)y(t)dt = x+\. 



X 

134. ^ s'm(x— t)(p(t)dt = 1 — cos >c. 

ft 

Integral equations of the first kind with logarithmic 
kernel 

J <p(t)\n(x — t)dl =/(*), /(0) = (8) 



can also be solved by means of the Laplace transformation. 
We know that 

... • Hv + l) 



(Rev>-1) (9) 



Differentiate the relation (9) with respect to v: 



dv ^ p» + 1 p 

or 

"dr (v+1) 

, lnjt!?5! 



p 

For v = we have (see p. 54) 

r (1) 



L r(v+i) + in p, 



T(I)~ 7 
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where y is Euler's constant and formula (10) takes the form 
ln* = l(- T -ln P ) = -!M±! (11) 

Let cp(x) = (D(p), f (x) = F (p). Taking the Laplace trans- 
form of both sides of (8) and utilizing formula (11), we get 



whence 



_4, (p) i5it2 = f W 



Let us write <I>(p) in the form 

p/pi- p^(p)-/ , (Q) no) n3 v 

ww P(lnp + v) Pdnp+Y) U ' 

Since /(0) = 0, it follows that 

p 1 F(p)-r(0) = f{x) (14) 



Let us return to formula (9) and write it in the form 

JC" . 1 

r (v+i)" 7_ p i,+i 



(9') 



Integrate both sides of (9') with respect to v from to oo. 
This yields 

c x* . r d\ __ i 
J r(v+i) flv -^J pv + i- plnp 



By the similarity theorem 

? * Vg ~ v dv= 1 = ! 

Jr(v+l) uv - pin (op) p(lnp+lna) 
o 

If we put a = eT, then 

J T(v + 1) ~p(lnp + v) ^ 
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Take advantage of equality (13). By virtue of (15) 



p(lnp + Y ) • ' v 'J T(v+1) 
o 



dv 



Taking into account (14) and (15), the first term on 
the right of (13) may be regarded as a product of transforms. 
To find its original function, take advantage of the convo- 
lution theorem: 

MftzLB ='[r(t)( f dv) dt 

p(ln p + v) J w I J r(v+l) ) 
M> / 

Thus, the solution <p(x) of the integral equation (8) will 
have the form 



where y is Euler's constant. 

In particular, for f(x) — x we get 



<p(*) = — inv+T) dv 



The convolution theorem can also be used for solving 
nonlinear Volterra integral equations of the type 

X 

<p (*) = /(*) + A, J<p(0<P(*— t)dt (16) 
o 

Let 

<p(x\=<b(p), f(x) = F(p) 

Then, by virtue of equation (16), 

a)(p) = F(p) + X0*(p) 

whence 

1 ± }T\-4lF(p) 
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The original function of <t>(p), if it exists, will be a so- 
lution of the integral equation (16). 
Example. Solve the integral equation 

X 

j«p(/)<P (*-*)<« = ■!- M7) 
o 

Solution. Let (p(x) (p). Taking the Laplace transform 
of both sides of (17), we get 

^(P)^ 

whence 

*(p) = ±jr 

The functions q>, (x) = x, (f i (x) — ~x will be solutions of 
the equation (17) [the solution of equation (17) is not uni- 
que]. 

Solve the integral equations: 

X 

135. 2cp (x) — J q> (t) cp (x— t) dt = sin x. 



X 

1 r> 1 

136. (})(a:)=2- \ q>(/)q)(;t — t)dt — ysinhx 
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12. Fredholm Equations of the Second Kind. Fundamentals 

A linear Fredholm integral equation of the second kind is 
an equation of the form 

b 

<p(x) — X^K(x,t)<p(t)dt = f(x) (1) 

a 

where cp(x) is the unknown function, K (x, t) and f (x) are 
known functions, x and / are real variables varying in the 
interval (a, b), and X is a numerical factor. 

The function K(x, t) is called the kernel of the integral 
equation (1); it is assumed that the kernel K (x, I) is defined 
in the square Q {a^.x^b, a^l^Lb) in the (x, t) plane 
and is continuous in Q, or its discontinuities are such that 
the double integral 

b b 

J \ \K(x,t)\ 2 dxdt 

a a 

has a finite value. 

I f f (x) ^ 0, equation ( 1 ) is nonhomogeneous; but if / (x) e= 0, 
then (1) takes the form 

b 

(f(x) — l < \ l K(x,t)(p(t)dt=0 (2) 

a 

and is called homogeneous. 

The limits of integration, a and b in equations (1) and (2), 
can be either finite oi infinite. 

A solution of the integral equations (1) and (2) is any 
function 9 (x) which, when substituted into the equations, 
reduces them to identities in x£(a, b). 

Example. Show that the function (p(jK) = sin y is a solu- 
tion oi the Fredholm-type integral equation 

<Ux)-£§K(x, t)<f>(t)dl=± 
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where the kernel is of the form 

/ *(2-/) 

K(x, t) =\ t(2 2 _ x) 



2 

Solution. Write the left-hand side of the equation as 
i 

( 

= q>W- j || K (x, /) q) (0 dt + j (X, f) <p (0 <tt J = 
= ¥ W "? (1 ^ <P «> * + 1 ^ 9 «> * j = 

, X I ^ 

I o x > 

SIX 

T 



Substituting the function sin^ in place of q>(x) into this 



expression, we get 

* ..... nt I .._ 
sin £ — £ { (2 - x) \ i dt + x I (2 - 1)—^- dt } = 



2 r sm "9 r sin ' 

5in f-^ j(2-x) J * + (2-0— 



. h n' L . / t nt . 2 . n<\ |<=* , 

= sin T-T {< 2 ~ *> V - « cos T + J? sin T j U + 

T 2—t nt 2 . ntl |<=M a; 

+ X cos -o iSin-s- > =- Fr 

1 L 11 2 ji 2 2 J 2 

Thus, we have y = y, which, by definition, implies 

that (p (jt) = sin t£ is a solution of the given integral equa- 
tion. 

Check to see which of the given functions are solutions 
of the indicated integral equations. 

i 

137. (p(x)=l, <p(x) + \x{e xt — l)q>(t)dt = e x — x. 
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138. (f>(x) = e x (2x— I 



<p(x) + 2 < j j e*- t (p(t)dt = 2xe x 





Jl 

139. (p(x)=l— ^ilLi, <i,(x)—[cos(x + t)<p(t)dt = l. 

1 FT J 

-5 

1 

140. cj> (.*;) = Kx, (p(x) — $K(x, t)q>(t)dt = 

o 

= KI + ^(4x'/ 2 -7). 



2 

1 



141. <p (*) = «*, cp(x) + X J sinxt <p(t)dt = 1. 

o 

142. cp(x) = cosx, <p (x) — (* 2 + cos tq> (t)dt = sin x. 

o 

143. q>(x)=*xe-*, <p(x) — i\e- lx+t) <p(t)dt = (x—l)e- 



o 



144. <p (x) = cos 2x, cp(x) — 3^ K(x, t)<p(t)dt = cos x, 



K(x, t)-. 





sinxcos^, ^ x ^ 
sin t cos x, t ^ x ^ ji. 



145. (p(x) = -^-sinx, where C is an arbitrary constant, 

00 

<p (x) — — J sin x-^-^-<p(/)d/=0. 



13. The Method of Fredholm Determinants 

The solution of the Fredholm equation of the second kind 

6 

(p(x)— l[K(x, t)<i>(t)dl = f(x) (1) 
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is given by the formula 

b 

<p (*) = /(*) + *,$/?(*, /; %)f(t)dt 



(2) 



where the function R(x, t; X) is called the Fredholm resol- 
vent kernel of equation (1) and is defined by the equation 



Rix , 



(3) 



provided that D(X)=£0. Here, D(x, t; X) and D(ty are 
power series in X: 



D(x, t; X) = K(x, t) + Y t L^B n (x. t)X», 

oo 



(4) 
(5) 



whose coefficients are given by the formulas 

K(x, t) K(x, t t ) ...K(x, t n ) 
K^, t)K(t v tJ.-.K^, t n ) 
K{t % , t)K(t t , tJ.-.K^, t n ) 



B n (x 



a a 
n 



K{t n , t)K(t n ,t 1 )...K(t n ,t n ) 



dt x ...dt n 
(6) 



and 




B (x, t) = K(x, t) 

Kit,, h) Kit,, Q...KU,, t n ) 
Kit 2 , t,) Kit,, Q ...Kit,, t n ) 
Kit a , t,) Kit 3 , t t ) ... Kit 3 , t n ) 



dt,...dt n (7) 



Kit n , tj Kit n , t,) ...Kit n 
The function Dix, t; X) is called the Fredholm minor, 



FREDHOLM INTEGRAL EQUATIONS 



75 



and D(X) the Fredholm determinant. When the kernel 
K(x, t) is bounded or the integral 

b b 

\\K 2 (x, t) dx dt 

a a 

has a finite value, the series (4) and (5) converge for all 
values of X and, hence, are entire analytic functions of X. 
The resolvent kernel 

,£>(*, t; X) 



R (x, t; X) ■■ 



D(X) 



is an analytic function of X, except for those values of X 
which are zeros of the function D(X). The latter are the 
poles of the resolvent kernel R (x, t; X). 

Example. Using Fredholm determinants, find the resolvent 
kernel of the kernel K(x, t) = xe u , a = 0, 6=1. 

Solution. We have B (x, t) = xe' Further, 



B t (x, = 



xe' xe' 1 
t,e' t,e'> 
xe' xe 1 ' xe u 
t^' t x e x > t 1 e u 
t,e' t 9 e u t.,e u 



dt, = 0, 



d/, dt 2 = 



since the determinants under the integral sign are zero. 
It is obvious that all subsequent B n (x, i)—0 Find the 
coefficients C n : 

i i 
C, = J K(t t , t l )dt i = \t l e , 'dt 1 = 1, 



1 I 



Le*' Le'> 



dt, dt, = 



Obviously, all subsequent C n are also equal to zero. 
In our case, by formulas (4) and (5), we have 

D(x, t; X)=K(x, l) = xe l \ D(X)--- \—X 

Thus, 

D(x, t; X) _ xe- 



R (X, I, A) = 



76 INTEGRAL EQUATIONS 

Let us apply the result obtained to solving the integral 
equation 

i 

<p(x)—klxe t tp(t)dt = f(x) (Jt#l) 



By formula (2) 

i 

<p(x) = f(x) + X$-^f(t)dt 

o 

In particular, for f(x) = e~ x we get 

(p( x ) = e -*_l__A_ x 

Using Fredhohn determinants, find the resolvent kernels 
of the following kernels: 

146. K(x, t) = 2x—t; 0<%< 1, 0<*<1. 

147. K (x, t)=xH—xt\ 0<x< 1, 0<* < 1. 

148. K(x, t) = smxcost; 0<x<2n, 0<^<2n. 

149. K(x, = sin*— sin*; < 2n, 0</<2n. 

Practically speaking, only in very rare cases is it possible 
to compute the coefficients B n (x, t) and C n of the series (4) 
and (5) from formulas (6) and (7), but from these formulas 
it is possible to obtain the following recursion relations: 

b 

B n (x, t) = C n K(x, t)-n\K{x, s)B n . 1 (s, t)ds, (8) 

a 

b 

C n = lB n .^ s)ds (9) 



Knowing that the coefficient C =l and B (x, t) = K{x, t), 
we can use formulas (9) and (8) to find, in succession, C lt 
B 1 (x, t), C 2 , fl 2 (x, /), C 3 and so on. 

Example. Using formulas (8) and (9), find the resolvent 
kernel of the kernel K(x, t) = x — 2t, where O^x^l, 
0</<l. 
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Solution. We have C = l, B (x, t) = x — 2t. Using for- 
mula (9), we find 

i 

C 1 = J(-s)ds = -i- 

By formula (8) we get 
i 

B, (x, t) = — £=21 __ j (x— 2s) (s— 20 ds = —x—/ + 2x/ + -| 



We further obtain 



C 2 = J(-2s + 2s 2 + |-)ds = l 



1 

B,(x, t) = Z^—2 2s) (— s— / + 2s/ -f f)ds = 



C, = C 4 =....= 0, B 3 *) = -64 (*. = • • • = 
Hence, 

D(^) = l+-£- + ^; D(x, X,)=jc— 2zf + (a: + ^ — 2xt— -|U 



The resolvent kernel of the given kernel is 

- x— 2t+ (yX + t — 2xt — j^j-K 

R(x, t; %) = r — 

> + T + T 

Using the recursion relations (8) and (9), find the resol- 
vent kernels of the following kernels: 

150. K(x, t) = x + t + l; — — 1</-<L. 

151. K(x, t)=\+3xt; 0<x< 1, 0</< I. 

152. K(x, t) = 4xt— x 3 ; 0<x<l, 0<<<1. 

153. K(x, t) = e x - t ; . 0<x< 1, 0</<l. 

154. K(x, t) = sin (x + t); 0<x^2jx, 0<f<2n. 

155. K(x, t) = x — s'mht; —l^x^l, — l</«£l. 
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Using the resolvent kernel, solve the following integral 
equations: 

156. y(x) — X) sin(x + /)(p(/)d/ = 1. 



1 

157. <p(Jt) — A. §(2x— t)q>(t)dt = ^ . 

a 

in 

158. (p(x) — ^ sin xcos / q> (t)dt — cos 2x. 



1 

159. <p(x)+le x -\(t)dt = e*. 



1 

160. y(x) — X ^(4xt— x 2 )<p(t)dt=x. 



14. Iterated Kernels. Constructing the Resolvent Kernel 
with the Aid of Iterated Kernels 

Suppose we have a Fredholm integral equation 

b 

y(x) — 'k\K(x,t)<f>(t)dt = f(x) (1) 

a 

As in the case of the Volterra equations, the integral equa- 
tion (1) may be solved by the method of successive appro- 
ximations To do this, put 

ce 

q> (*) = /{*)+ 2 *„(*) ^" (2) 

1=1 

where the ^ n (x) are determined from the formulas 

b 

(x) = J K (x, t)f(t)d(, 

a 

b b 

Mp,(x) = \K(x, t)$,(t)dt=\KAx> t)f(t)dt, 

a a 
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b b 

r|> s (x)=\K (x, t) t|) 2 (/) dt = J K 3 (x, t) f (t) dt 



and so on. 
Here 



K* (x, t) = \ K (x, z) /(,(>, t)dz, 



a 

and, generally, 

K n (x, t) = \K(x, z)K n .Az, t)dz (3) 

a 

n = 2, 3, .... and K x {x, t)==K(x, t) The functions 
K n (x, t) determined from formulas (3) are called iterated 
kernels For them, the following relation holds: 

A 

K n (x, t) = \K m (x, s)K n . m (s, t)ds (4) 

a 

where m is any natural number less than n. 

The resolvent kernel of the integral equation (1) is deter- 
mined in terms of iterated kernels by the formula 

CD 

R(x, t; X)=.2 K n (x, O^" 1 (5) 

n=\ 

where the series on the right is called the Neumann series 
of the kernel K(x, t). It converges for 



F b b 

where S = ]/ J J K 2 (x, t) dxdt 

a a 

The solution of the Fredholm equation of the second 
kind (1) is expressed by the formula 

<${x) = f(x) + x\R{x, t;X)f(t)dt (7) 
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The boundary (6) is essential for convergence of the se- 
ries (5). However, a solution of equation (1) can exist for 

values of | k | > as well. 

Let us consider an example: 

i 

<P(jc)— X$ <p(t)dt = l (8) 



Here K (x, = 1. an d hence 

it ii 
B* = J J /(» ( Xi t)dxdt = \\dxdt = \ 

oo oo 



Thus the condition (6) gives that the series (5) converges for 

IM<i- 

Solving equation (8) as an equation with a degenerate 

i 

kernel, we get (1 — K)C=l, where C=jj<p(/)d/. 



For X=l, this equation is unsolvable and hence for 
X = 1 the integral equation (8) does not have any solu- 
tion. From this it follows that in a circle of radius greater 
than unity, successive approximations cannot converge for 
equation (8). However, equation (8) is solvable for |X|> 1. 

Indeed, if %=£\, then the function <p(*) = t~ ' s a so ' u " 
tion of the given equation. This may readily be verified by 
direct substitution. 

For some Fredholm equations the Neumann series (5) 
converges for the resolvent kernel for any values of X. Let 
us demonstrate this fact. 

Suppose we have two kernels: K(x, t) and L(x, t). We 
shall call these kernels orthogonal if the following two 
conditions are fulfilled for any admissible values of x and 

b b 

\K(x,z)L(z,t)dz = 0, J L(x, z)K(z, t)dz = (9) 

a a 

Example. The kernels K{x, t)=xt and L(x, t)=xH 2 are 
orthogonal on [ — 1, 1]. 
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Indeed , 

i i 

J (xz)(zH l )dz = xt* J z 3 dz = 0, 
-i -i 
i i 

5 (x 2 z 2 ) (zt) dz = xH ^ z s dz = 
-i -i 

There exist kernels which are orthogonal to themselves. 
For such kernels, K 2 (x, t) = 0, where K 2 (x, t) is the second 
iterated kernel. It is obvious that in this case all subsequent 
iterated kernels are also equal to zero and the resolvent 
kernel coincides with the kernel K(x, t). 

Example. K(x, *) = sin (x— 2t); 0<x<2ji, 0<^<2ji. 

We have 

211 

^ sin (x — 2z)sin (z — 2t)dz = 



211 

= y [ cos (x -f 2t — 3z) — cos (x — 2t — z)]dz = 

o 

= y [— ysin (* + 2* — 3z) + sin(x— 2/— z)j |*"'* = 

Thus, in this case the resolvent kernel of the kernel is equal 
to the kernel itself: 

R(x, t; K) = sm{x—2t) 

so that the Neumann series (5) consists of one term and, 
obviously, converges for any X. 

The iterated kernels K„ (x, t) can be expressed directly in 
terms of the given kernel K(x, t) by the formula 

b b b 

*.(*,/)= SJ - S 2 ) ... 

a a a ' 

K (s n _ it t) ds lt ds 2 ...ds n _j - (10) 

All iterated kernels K„ (x, t), beginning with K 2 (x, t), 
will be continuous functions in the square a^.x^.b, 
a<*<6 if the initial kernel K{x, t) is quadratically sum- 
mable in this square. 
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If the given kernel K(x, t) is symmetric, then all itera- 
ted kernels K„(x, t) are also symmetric (see [15]). 

The following are some examples in f.nding iterated ker- 
nels. 

Example I. Find the iterated kernels for the kernel 
K(x, t) = x—t if a = 0, b= 1. 
Solution. Using formulas (2), we find in succession: 
t) = x — t, 



i 

K t (x, = jV-s) (s-t)ds- 



■ x+t -xt- 1 



2 3 ' 

x — t 
12 ' 



KAx. t) = ^(x-s) (i±i-st-^ds = 

o 

i 

K t (x, t) = - ^ j (x— s) (s - ds = 

K s (x, t) = -~^(x-s) (t±i- st -±y s = 



12" V 2 3 y • 

From this it follows that iterated kernels are of the form: 

(1) for n = 2fe— 1 

(2) for n = 2k 

^.o-t^f-f-*-*) 
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where ft= 1, 2, 3, . . . 

Example 2. Find the iterated kernels K, (x, () and K 2 (x, t) 
if K (x, t) = e mi " " '>, a = 0, 6= 1. 

Solution. By definition we have 



min (x, t) = { 

\ t, if t ?C x?C \ 



and for this reason the given kernel may be written as 

I e x , if 0<x</, 
I e , if * <x< 1 

This kernel, as may easily be verified, is symmetric, i. e., 

K(x, t) = K(t, x) 

We have (x, t) = K(x, t). We find the second iterated 
kernel- 

i i 
K, (x, t) = \ K (x, s) /(, (s, f) cfs = J K (x, s) K (s. t) ds 

n 

Here 

\ e x , if 0<x<s, * 

| e\ if s<x< 1, 

I e 5 , if OsCss^/, 
* (S '' ) = U if 

Since the given kernel K(x,() is symmetric, it is sufficient 
to find K*(x, t) only for x> t. 
We have (see Fig 2) 

K 2 (x, t)=^K(x, s)K(s, t)ds + [K (x, s)K(s, l)ds + 

/ 

I 

+ \ K (x, s) K (6, i)ds 
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H 1 •- 

t X S 

Fig. 2 



In the interval (0, t) we have s < / < x, and therefore 
t t t 

§K(x, s)K(s, t)ds={ j e s e s ds = ^e 2S ds-- 



e 2 ' — I 



2 

In the interval (/, x) we have / < s < x, and therefore 

X X 

\K{x, s)K(s, t)ds=\e s e l ds = e x+t —e zt 
t t 

In the interval (x, 1) we have s>x> t, and therefore 
i i 
J K(x, s)K(s, t)ds=[e x e t ds = (\—x)e x+t 

X X 

Adding the integrals thus found, we obtain 

K % {x, t) = {2-x)e*+ t — (x>t) 

We will find the expression for K 2 (x, t) for x < / if we 
interchange the arguments x and / in the expression K 2 {x, i) 
for x > t: 

Kz(x, t) = (2-t)e x+i -^(x < t) 

Thus the second iterated kernel is of the form 

(2 — *)e*+'_!±^, if 0<x<^ 
{2—x)e x+i —-^-, if 
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Note. If the kernel K (x, t), which is specified in the 
square a^.x^.b, a^.t^.b by various analytic expressions, 
is not symmetric, then one should consider the case x < / 
separately. For x < t we have (see Fig. 3) 

b x t b 

K t {x, t) = \K(x, s)K(s, t)ds= J + S + S 

a a x t 



a x t s b 

Fig. 3 

Example 3. Find the iterated kernels /C x (x, t) and K.< (x, t) 
if a = 0, 6=1 and 



K(x, t) = 



x + t, if 0^x<t, 
x—t, if t <xs^\ 
Solution. We have K t (x, t) = K (x, t), 



K % (x, *)=$*(*, s)K(s, t)ds 

o 

where 

I x + s, 0<*<s, I s + t, 0<s</, 

K (x, s) = < K (s, t) = < , 

Since the given kernel K(x, t) is not symmetric, we con- 
sider two cases separately when finding K 2 (x, t): (1) x < t 
and (2) x> t. 

( 1 ) Let x<t. Then (see Fig. 3) 

K 2 (x, t) = I 1 + I i + I 3 

where 

X 

I 1 =\(x-s)(s + t)ds = ^--\ 
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t 

l. l = §(x + s)(s + t)ds = 5 -£-™+jXt*-^xH, 

x 
1 

l, = $(x+s)(s-.t)ds = ^+%-xt + ±—L + ± 
t 

Adding these integrals we obtain 

KAx, = t 3 —-jX 3 —xH + 2xl 2 — xt + + j (x<0 
(2) Let x > Then (see Fig. 2) 

K t (x, t) = l 1 + l 1 + l 3 

where 

i 

I 1 = §(x-s)(s + t)ds = ^xl*-^, 



/ I = J(x-s)(s-/)ds = 1 — -g — g" + -g-. 
i 

l 3 = §(x + s)(s-t)ds= -^x* + ^xH + ^-—xt + j 

Adding these integrals, we obtain 
KAx, t)= -jx 3 -t° + xH+2xt 2 — xt + ^ + j (x>t) 

Thus, the second iterated kernel is of the form 
K t (x, t) = 

j — jX* + t 3 -xH + 2xl 2 -xt+-^--\-j, 0<x<^ 
( —jX 3 — t 3 + xH+2xt* — xt + ^- + j, t<x^\ 

The other iterated kernels K„(x, t)(n = 3, 4, ...) are 
found in similar fashion 

Find the iterated kernels of the following kernels for 
specified a and b. 

161. K(x, t) = x— t; a=—\, b=\ 

162. K(x, t) = sin (x—t); a = 0, = {n = 2, 3). 
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163. K(x, t) = (x—t)*; a= — 1, b= 1 (n = 2, 3). 

164. K(x, t) = x+s'mt; a=—n, b = n. 

165. K(x, t) = xe u , a = 0, 6=1. 

166. /C(x, /) = e*cos/; a = 0, 6 = n. 

In the following problems, find K^ix, t): 

167. K(x, = 6"-"; a = 0, 6= 1. 

168. K(x, t) = e^^'\ a= — 1, 6=1. 

We now give an instance of constructing the resolvent 
kernel of an integral equation with the aid of iterated kernels. 
Consider the integral equation 

<p (x) — X \xtif{t)dt = f{x) . (11) 
o 

Here K (x, t) = xt; a = 0, 6=1. In consecutive fashion we 
find 

Ktix, t) = xt, 
i 

KAx, t) = ^(xz)(zt)dz = ^-, 



I 

K 3 (x, t) = ^{xz)(zt)dz = ^. 



3xt 
3—X 



According to formula (5) 

R(x, t; \) = ^K n (x, O^-'-^S (tP ,= 

where | X | < 3. 

By virtue of formula (7) the solution of the integral 
equation (11) will be written as 

i 

<P (*)=/(*) + * J 3^X /CM' 
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where X 3. 

Construct resolvent kernels for the following kernels: 



If M (x, t) and N (x, t) are two orthogonal kernels, then 
the resolvent kernel R(x, t\ X) corresponding to the kernel 
K(x, t) = M-\-N, is equal to the sum of the resolvent 
kernels Rx{x, t; X) and R 2 (x, t; X) which correspond to 
each of these kernels. 

Example. Find the resolvent kernel for the kernel 



Solution. As was shown above, the kernels M (x, t) = xt 
and N (x, /) = x 2 / 2 are orthogonal on [ — 1, 1] (see p. 80). 
For this reason the resolvent kernel of the kernel K (x, t) is 
equal to the sum of the resolvent kernels of the kernels 
M(x, t) and N (x, t). Utilizing the results of problems 173 
and 174, we obtain 



R K (x, t; X) = R M (x, t; X) + R N (x, t; X) + T 



Find the resolvent kernels for the kernels: 

175. K(x, t) = sin x cos / + cos 2x sin 2t; a = 0, b = 2n. 

176. K(x, 0= 1 +(2x— l)(2i — 1); a = 0, 6=1. 

This property can be extended to any finite number of 




a = 0, 6=1. 
fl = 0, 6 = f. 
c= —1, 6 = 1. 
a= — 1, 6 = 0. 
a= — 1, 6=1. 
a= — 1, 6=1. 



K{x, f)=xt + x*t\ a=— 1,6=1 



3 

where | X | < -j . 



kernels. 



If the kernels M (1 > (x, t), M <2 > (x, t), .... M (n) (x, t) are 
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pairwise orthogonal, then the resolvent kernel corresponding 
to their sum, 

K(x, /)= 2 M™(x, t) 

m= 1 

is equal to the sum of the resolvent kernels corresponding 
to each of the terms. 

Let us use the term "nth trace" of the kernel K (x, t) for 
the quantity 

b 

A n =\K n {x, x)dx, (n=l, 2, ...) (12) 



where K n (x, t) is the nth iterated kernel for the kernel 
K(x, t). 

The following formula holds for the Fredholm determi- 
nant D(X): 

™ ±A n %°-> (13) 

The radius of convergence of the power series (13) is equal 
to the smallest of the moduli of the characteristic numbers. 

177. Show that for the Volterra equation 

X 

y(x)—k J K(x, t)y(t)dt=f(x) 



the Fredholm determinant D (X) = e~ A ^ and, consequently, 
the resolvent kernel for the Volterra equation is an entire 
analytic function of X. 

178. Let R(x, t; X) be the resolvent kernel for some 
kernel K(x, t). 

Show that the resolvent kernel of the equation 

<p(x) — |i lR(x, t; l)y(t)dt = f(x) 

a 

is equal to R(x, t; k + n). 
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179. Let 



J J/C 2 (*. t)dxdt = B\ 



a a 



b b 



J J/CJU*. t)dxdt = B* 



a a 



where K n (x, t) is the nth iterated kerne) for the kernel 
K (x, t). Prove that if fi 2 = B 2 , then for any n we will 
have fi„ = fi". 



15. Integral Equations with Degenerate Kernels. 
Hammerstein Type Equation 

The kernel K(x, t) of a Fredholm integral equation of 
the second kind is called degenerate if it is the sum of a 
finite number of products of functions of x alone by func- 
tions of t alone; i.e., if it is of the form 



We shall consider the functions a k (x) and b k (t) (k = l, 
2, ...,«) continuous in the basic square a^.x, i^.b and 
linearly independent. The integral equation with degenerate 
kernel (1) 



n 



K(x, = 2 a k (x)b k (i) 



(1) 




is solved in the following manner. 
Rewrite (2) as 



n 




(3) 



and introduce the notation- 




(4) 



a 



FREDHOLM INTEGRAL EQUATIONS 



91 



Then (3) becomes 



q> (*) = /(*) + A, 2 C k a k (x) 

k = 1 



(5) 



where C k are unknown constants, since the function cp (x) is 
unknown. 

Thus, the solution of an integral equation with degene- 
rate kernel reduces to finding the constants C k (k=\, 
2, n). Putting the expression (5) into the integral 
equation (2), we get (after simple manipulations) 

b 



MO 



f(t)+x£c& k (t) 



dt}a m (x) = 



Whence it follows, by virtue of the linear independence of 
the functions a m (x) (m = 1, 2, n), that 



K (0 



f(t) + Xy C k a k {t) 



or 



C m J ( C k J a k (t) b m (t) dt = S b m (t)f(t)dt(m- 1 , 2 n) 



For the sake of brevity, we introduce the notations 

b b 

a km =\a h (t)b m {t)dt, f m =\b m (t)f(t)dt 

a a 

and find that 

C m -X 2 a km C k = f m (m = l, 2, n) 

or, in expanded form, 

( 1 — Xa u ) C, — Xa, 2 C S — . . . — Xa ln C„ = f v 

— Xa 2l C l + (\-'ka ii ) C 2 — . . .—Ka tn C n = f 2 , 

- Xa nl C, - - . . . + ( 1 - Xa nn ) C n = /„ 



(6) 
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For finding the unknowns C k , we have a linear system of n 
algebraic equations in n unknowns. The determinant of this 
system is 

1— Ka„ —Xa„ ... — Xa ln 
— %.a, 



A(X) = 



— %a 21 l—Xa M 



\—%a„ 



(7) 



If A(a.)=^0, then the system (6) has a unique solution 
Cj, C 2 , .... C„, which is obtained from Cramer's formulas 

1— Xa n ... — "ka lk _J x — Ka lk+i ... — hi ln 
— %a 2l ... — ka ik _J i — Xa 2ft+1 ... — Xa 2 „ 



AW 



(8) 



*«i ••• — ha nk _ 1 f„— ha nk+i ... 1— Xa„ 
(k=\, 2, n) 

The solution of the integral equation (2) is the function 
<p(x) defined by the equality 

n 

<p(x) = f(x) + % S C ft a ft (x) 
k— i 

where the coefficients C ft (fe=l, 2, n) are determined 
from formulas (8). 

Note. The system (6) may be obtained if both sides of 

(5) are consecutively multiplied by a t (x), a 2 (x) a n (x) 

and integrated from a to b or if we put (5) into (4) for 
(p(x), replacing x by 

Example. Solve the integral equation 
n 

(p(x)—K J (x cos t + 1 2 sin x + cos x s'm t) q> (t) dt = x (9) 
Solution. Write the equation in the following form: 

31 Jl 

<p(x) = kx^ q> (t) cos tdt + K sin x \ t 2 y(t)dt + 

-31 -31 
31 

-J-a-cosx ^ y(t) sintdt+x 
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We introduce the notations 

ji n n 

C 1= 5 <p (t) cos tdt; C 2 = J t\{t)dt; C 3 = J <p(0sin<d/ 

-ji -ji -ji 

(10) 

where C x , C 2 , C 3 are unknown constants. Then equation (9) 
assumes the form 

(p (x) = C x Xx + C 2 X sin x + C 3 X cos x + x (11) 

Substituting expression (11) into (10), we get 
ji 

C t = J (C^ + C^sinf + CjXcosf-f-Ocos/d/, 
-ji 
it 

C 2 = J (C 1 M + C 1 Xsin/ + C t A,cos/ + 0* , d/, 
ji 

C 3 = ^ (C X M + C 2 k sin f + C 3 K cos / + /) sin / dt 

-n 

or 

/ ji \ ji 

CJ1— A, J f cos tdt )— C 2 X. J sin * cos tdt— 
\ -jt ' -it 

ji ji 

— C 3 a. J cos s /d* = J t cos tdt 
-ji -it 

— C,a J / 3 d/ + C 2 ^l — X. J t 2 sin tdt)— 

ji ji 

— C 3 X J Pcostdt= J / 3 df, 
-ji -n 
ji ji 

— CjX J t sin tdt—C^K jj sin 2 *d/ + 

— ji -ji 

(JI v Jt 

1— X J cos/sin/d/] = $ rfsin/d* 
-n / -ji 
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By evaluating the integrals that enter into this system we 
obtain a system of algebraic equations for finding the un- 
knowns C,, C 2 , C s : 

C, — XnC 3 = 0, 
C, + 4XnC 3 = 0, 
— 1 AJtC, — ajiC, + C 3 = 2ji 



(12) 



The determinant of this system is 
I — rik 
A (X) = 1 4nX 
—2nX — Xn 1 



= 1 + 2XW 



The system (12) has a unique 'olution 

Ll _ 1 +2A«n« ' 2_ l+2A 2 n 2 * 3_ I-|-2X 2 ji 2 

Substituting the values of C v C 2 , C 3 thus found into (11), 
we obtain the solution of the given integral equation 

(p(x) = 1 4. 2X 2 ji 2 — 4Xn s'm x + cos x) -t- x 

Solve the following integral equations with degenerate 
kernels: 



180. cp(JC)— 4$ sin 2 ;t(p (/) dt = 2x— n. 



1 

181. <p(x)— J e arc5ln * <p (/) d/ = tan jc. 

-1 

ji 
4 

182. <f>(x) — X J tan/cp(Od/ = cot;t. 

4 

I 

183. <p(x)— X ^ cos {q\nt)q>(t)dt = I. 


1 

184. (f(x)— X J arc cos /q> (/) = -p===- . 
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1 

185. <p(x) — I j (ln-j-J%(0^ = 1 (p>-l). 



I 

186. <p(x) — A, j(xln* — Mnx)<p(/)d/ = |-(l— 4x). 



o 



187. (p(x) — X J sin xcos /q> (/)d/ = sinx. 

o 

188. <p(x) — a, ^l 11 — * I sin xcp (t)dt = x. 



189. <p(x) — X J sin(x— t)y(t)dt = cosx. 



190. g>(x)- 



— X ^ (sin x cos < — sin 2x cos 2* 4- sin3xcos3/) <p(t)dt =cosx. 



191. q>(x) — 



- T J (*-|(3' 2 -l) + y'(3x*- 1) 



q>(/)d/ = 1 



Many problems of physics reduce to nonlineai integral 
equations of the Hammerstein type (see [24], [28]). 

The canonical form of the Hammer stein-type equation is 



(1> 



where /((x, /), /(i, u) are given functions and <p(x) is the 
unknown function. 

The following equations readily reduce to equations of 
type (1): 

b 

q>U)= $/<(*, /)/«. <p(/))d/ + i|>(x) (!') 
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where ^(x) is the known function, so that the difference 
between homogeneous and nonhomogeneous equations, which 
is important in the linear case, is almost of no importance 
in the nonlinear case. We shall call the function K(x, t) 
the kernel of equation (1). 
Let K(x, t) be a degenerate kernel, i. e., 

m 

K(x,t)=2a i (x)b i (t) (2) 
/ = i 

Then equation (1) takes the form 

m b 

<p(*)=2a,(x) lb,(t)f(t,<p(t))dt (3) 

' = 1 a 

Put 

ft 

Ct=lb t (t)f(t,<p(t))dt (t = l, 2 m) (4) 

a 

where the C,- are as yet unknown constants. Then, by virtue 
of (3), we will have 

m 

fW=2C ft W (5) 
/= 1 

Substituting into (4) the expression (5) for (f>(x), we get m 
equations (generally, transcendental) containing m unknown 
quantities C lt C 2 C m : 

C^Vtfv C 2 CJ (t = l, 2 m) (6) 

When /(/, u) is a polynomial in u, i. e., 

f(t,u) = p (t) + p l (t)u+...+p n (t)W (7) 

where p (t), . . ., p n (t) are, for instance, continuous functions 
of t on the interval [a, b], the system (6) is transformed 
into a system of algebraic equations in C u C 2 , C m . 
If there exists a solution of the system (6), that is, if there 
exist m numbers 

CO pO /->0 
1) l -'2> • • • > 

such that their substitution into (6) reduces the equations 
to identities, then there exists a solution of the integral 
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equation (3) defined by the equality (5): 

m 

<p(*)= 2 Qa t (x) 
i= i 

It is obvious that the number of solutions (generally, 
complex) of the integral equation (3) is equal to the number 
of solutions -of the system (6). 

Example. Solve the integral equation 

i 

q>(x) = xlxt<p*(t)dt (8) 
o 

where % is a parameter. 
Solution. Put 

i 

C = J t<p* (0 dt (9) 
o 

Then 

<p(x) = Ckx (10) 

Substituting q>(x) in the form (10) into the relation (9), we get 

i 

C = J dt 
o 

whence 

C=^-C» (11) 
Equation (11) has two solutions 



Cj — 0, C 2 — -jj 



Consequently, integral equation (8) also has two solutions 
for any X 0: 

(Pi (x) = 0, q> t (x) — 

There exist simple nonlinear integral equations which do not 
have real solutions at all. 
Consider, for example, the equation 

1 x+t 

<p(x)=^e~(l + <p*(t))dt (12) 



4-15 83 
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Put 

i r - 

C== TJ e2 (1+9* (0)* (13) 



Then 

<p(x) = Ce T (14) 
For a determination of the constant C we obtain the equation 

(e^ — l)c a — 3C + 3(e^"— 1)=0 (15) 

It is easy to verify that equation (15) does not have real 
roots and, hence, the integral equation (12) has no real 
solutions. 

On the other hand, let us consider the equation 
i 

<p(x) = $a(x)a(t) v (t)s\n($fc)di (16) 



(a(0>0 for all <g [0, 1]) 

In order to determine the constant C, we arrive at the 
equation 

i 

1 = J o» (/)■*■ sin C (17) 
i) 

i 

If ^a*(t)dt>l, then equation (17) and, hence, the original 



integral equation (16) as well have an infinite number of 
real solutions. 
Solve the following integral equations: 
i 

192. <p(x) = 2$x/g> s (0#. 



1 



193. q>(*)= l(xt + xH 2 )<p 2 (t)dt. 

-l 
i 

194. <f(x)*=l&Pip*{t)dt. 
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r xt 

195. q>(*)= ) i + 9 i(/) at. 

-l 
l 

196. ff>(x)=l(l+q> 3 (t))dt. 

i) 

197. Show that the integral equation 

<p(x) = jja(x)a(t)(l+<p*(t))dt 



(a(x) > for all x£ [0, 1]) 
i 

has no real solutions if ^a 2 (x)dx> 1. 



16. Characteristic Numbers and Eigenf unctions 

The homogeneous Fredholm integral equation of the se- 
cond kind 

6 

<p(x) — X J K (x, t) <p (t)dt = Q (1) 

a 

always has the obvious solution <f>(x) = 0, which is called 
the zero (trivial) solution. 

The values of the parameter X for which this equation 
has nonzero solutions <p(x)^0 are called characteristic num- 
bers* of the equation (l)or of the kernel K(x, t), and every 
nonzero solution of this equation is called an eigenf unction 
corresponding to the characteristic number X. 

The number X = is not a characteristic number since 
for a. = it follows from (1) that <jp(x) = 0. 

If the kernel K(x, t) is continuous in the square fi{a^x, 
t^b) or is quadratically summable in Q, and the numbers 
a and b are finite, then to every characteristic number X 



* Some authors use "eigenvalues" in place of the term "characte- 
ristic numbers". We use eigenvalue for the quantity <* = ^< where X, is 
the characteristic number. 
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there corresponds a finite number of linearly independent 
eigenfunctions; the number of such functions is called the 
index of the characteristic number. Different characteristic 
numbers can have different indices. 
For an equation with degenerate kernel 

<P (x) - 1 J [ J> ft W h (f )] 9 (0 = (2) 

the characteristic numbers are roots of the algebraic equation 
I— %a u — taz 12 ... — %a ln 
— Xa 



4W = 



1— -tai 2i! . . . — Xa 2 „ 







(3) 



the degree of which is p^n. Here, A (X) is the determinant 
of the linear homogeneous system 



(l-Xa 11 )C l -ka lt C t -...-*a lH C n = 0, 
— ^a 21 C 1 + ( 1 — Xa M ) C 2 — . . . — Xa 2 „C„ = 0, 



(4) 



■Ka nl C 1 — Xa n2 C i 



+ (1-K„)C„ = 



where the quantities a mk and C m (fe, m=l, 2 n) have 

the same meaning as in the preceding section. 

If equation (3) has p roots (l^p<n), then the integ- 
ral equation (2) has p characteristic numbers; to each cha- 
racteristic number % m (m = l, 2, p) there corresponds 
a nonzero solution 

1 » 2 ' •••> W 1» 
/°(2) /°<2> ^(2) . 1i 



c«">, c</>, .... c<»— ^ 

of the system (4). The nonzero solutions of the integral 
equation (2) corresponding to these solutions, i.e., the eigen- 
functions, will be of the form 

(PiW = 2 c * 1>a *W. q>.W = 53Ci"a*(*), .... 
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An integral equation with degenerate kernel has at most n 
characteristic numbers and (corresponding to them) eigen- 
functions. 

In the case of an arbitrary (nondegenerate) kernel, the 
characteristic numbers are zeros of the Fredholm determi- 
nant D(X), i.e., are poles of the resolvent kernel R(x, t; 
X). It then follows, in particular, that the Volterra integral 

X 

equation q> {x) — X ^ K (x, t) q> (/) dt = where K (x, t) g L 2 (Q B ) 



has no characteristic numbers (for it, D(X) =e~^«\ see 
Problem 177). 

Note. Eigenfunctions are determined to within a multi- 
plicative constant; that is, if cp(x) is an eigenfunction cor- 
responding to some characteristic number X, then Cq> (x), 
where C is an arbitrary constant, is also an eigenfunction 
which corresponds to the same characteristic number X. 

Example. Find the characteristic numbers and eigenfunctions 
of the integral equation 



(p (x) — X J (cos 3 x cos 2t + cos 3x cos 3 /) <p (t) dt = 



Solution. We have 

n n 

(f>(x) — X cos 2 x^(p(t) cos 2td(+X cos 3x $ <p (t) cos 3 / dt 

. 

Introducing the notations 

it n 
C, = $ <p(t) cos 2/ dt, C 2 = J <*> (0 cos 3 1 dt 0) 



we get 

q) (x) = C i X cos 2 x + C a ^ cos Zx (2) 
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Substituting (2) into (1), we obtain a linear system of ho- 
mogeneous equations: 

C 1 (\— X J cos 2 /cos 2/ d/) — C 2 X $ cos 3/ cos 2/ d/ = 0, 



—C,X I cos 6 1 dt + C 2 (l — X $ cos 3 1 cos 3/ dt = 



(3) 



But since 

it n 

Jcos 2 /cos2/d/ = -j, Jcos3/cos2/d/ = 



ji n 

jcos 6 /d/ = 0, j cos 3 / cos 3/d/ = -J 



it follows that system (3) takes the form 
l_^\ Cl =0/| 



4 

T 



= 



(4) 



The equation for finding characteristic numbers will be 



Kn 



4~ 



1 — 



8 



= 



The characteristic numbers are X, —— , >,, = —. 

1 it ' n 

For X =— , system (4) becomes 

0-^=0, 
y-C 2 = 

whence C 2 = 0, C t is arbitrary. The eigenfunction will be 
(jp t (x) = Cj X. cos 2 x or, setting Cj X = 1 , we get tpj (x) = cos 2 x. 
g 

For X= — , system (4) is of the form 
i (-1) ^ = 0, 
\ 0C 2 = 
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whence C 1 =0, C t is arbitrary and, hence, the eigenf unc- 
tion will be <p 2 (x)=C 2 Xcos3x, or, assuming C 2 K = l, we 
get cp 2 (x) = cos 3x. 

Thus, the characteristic numbers are 




and the corresponding eigenfunctions 



A homogeneous Fredholm integral equation may, gene- 
rally, have no characteristic numbers and eigenfunctions, 
or it may not have any real characteristic numbers and 
eigenfunctions. 

Example 1. The homogeneous integral equation 



has no characteristic numbers and eigenfunctions. Indeed, 
we have 



<f>! (x) = cos 2 x, 



q> 2 (x) = cos 3x 




n 



<p(x) = X(Sx—2) J t<f>(t)dt 







Putting 




(1) 








(2) 




(3) 



i 



But since \ (3/ 2 




equation (3) yields C = and, 







hence, cp (x) = 0. 
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And so for any X, this homogeneous equation has only 
one zero solution <p(x) = and, hence, it does not have 
any characteristic numbers or eigenfunctions. 

Example 2. The equation 
i 

<p(x) — X J (Vxt — VTx) q>(t)dt=0 



does not have real characteristic numbers and eigenfunc- 
tions. 
We have 

q>(x) = C 1 lVx— C 2 Xx (1) 

where 

i i 
C t = J fq> (0 ett, C 2 = S Vt 9 (/) <tt (2) 



Substituting (1) into (2), we get (after some simple mani- 
pulations) the system of algebraic equations 



(i_^)c 1+ Ac 2 = o, 
-*c 1 + (i+^-)c 2 = o 



(3) 



The determinant of this system is 



1- 
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For real X it does not vanish, so that from (3) we get 
C x = and C 2 = and, hence, for all real X the equation 
has only one solution, namely, the zero solution cp (x) = 0. 
Thus, equation (1) "does not have real characteristic numbers 
or eigenfunctions. 

Find the characteristic numbers and eigenfunctions for 
the following homogeneous integral equations with dege- 
nerate kernels: 



198. (p(x)—X^ sin 2 x<p {t)dt = 0. 
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2n 

199. <p(x) — sin x cos t<p (t) dt = 0. 

o 

200. tp(jc) — %\ sin x sin <(p (t)dt = 0. 

o 

ii 

201. qj(A;)- x\ cos (x + <P (<) d/ = 0. 

o 
i 

202. <p (x)—X J (45x a In < — 9* 2 In x) <p (/) dt = 0. 

o 

i 

203. <p(x) — X^ (2xt — 4x 2 ) q> (t)dt = 0. 

o 

204. <p(x) — k J (5^ 3 + 4x 2 O<p(/)d/ = 0. 

-l 

205. (f(x) — X J (5xP + 4a; 2 / + 3^)(P(0^ = 0. 

-i 
i 

206. <p(x) — l I (jecosh* — t smhx)(p(t)dt = Q. 

-l 

207. <p(x)— % \ (a; cosh/ — * 2 sinh x) <j> (/) = 0. 

-l 
i 

208. <p(x) — k J (xcosh* — tcoshx)(p(t)dt = 0. 

-i 

If the nth iterated kernel K n (x, t) of the kernel K (x, t) 
is symmetric, then it may be asserted that K(x, t) has at 
least one characteristic number (real or complex) and that 
the nth degrees of all characteristic numbers are real num- 
bers. In particular, for the skew-symmetric kernel K(x, t) = 
——K(t, x) all characteristic numbers are pure imaginary 
X = where p is real (see Problem 220). 

The kernel K(x, t) of the integral equation is called 
symmetric if the condition K(x, t) = K{t, x) (a^x, t ^.b) 
is fulfilled. 
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The following theorems hold for the Fredholm integral 
equation 

b 

<P(x)—%\K(x, t)<p(t)dt = (1) 

a 

with symmetric kernel K (x, t): 

Theorem I. Equation (1) has at least one real charac- 
teristic number. 

Theorem 2. To every characteristic number X there corres- 
ponds a finite number q of linearly independent eigenfunc- 
tions of equation (1), and 

sup q < X 2 B 2 

where 

b b 

B 2 = J J /C 2 (x, t)dxdt 

a a 

The number q is called the index or multiplicity of the 
characteristic number. 

Theorem 3. Every pair of eigenf unctions «p 1 (x), cp 2 (x), 
corresponding to different characteristic numbers X t X it is 
orthogonal; i.e., 

S <Pi (*) W2 (x) dx = 

a 

Theorem 4. There is a finite number of characteristic 
numbers in every finite interval of the X-axis. The upper 
bound for a number m of characteristic numbers lying in an 
interval — / <.X < / is defined by the inequality 

m < l 2 B l 

When the kernel K(x, t) of the integral equation (1) is 
the Green's function of some homogeneous Sturm-Liouville 
problem, finding the characteristic numbers and eigenfunc- 
tions reduces to the solution of the indicated Sturm-Liou- 
ville problem. 
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Example. Find the characteristic numbers and eigenfunc- 
tions of the homogeneous equation' 



where 



n 

<p(x) — X$/C(jc, t)<p(t)dt = 
( 

K(x, = { 



cos x sin t, Q^x^Lt, 
cos / sin x, t ^ x <; n 



Solution. Represent the equation in the form 

x n 

cp(jc) = A, $/((*, t)y(t)dt + KlK(x, t)y(t)dt 



or 



<p(x:) = X-sin x J (p(t) cos / dt + X cos x ^ 9 (t) sin / dt (1) 

X 

Differentiating both sides of (1), we get 

X 

q/ (x) = X cos x J (p (0 costdt + X sin a: cos x(p (a;) — 


n 

— ^ sin x ^ <p(/) sin tdt — X sin a: cos xy(x) 



or 



<p' (a:) = X cos a: cp (/) cos * dt — X sin a: $ <p (/) sin / dt (2) 

* 

Differentiating again, we get 

X 

q>" (x) = — X sin x J q> (/) cos / d/ + X cos 2 *(p (a;)— 


— a. cos x J <P (0 sin / dt + X sin 2 x cp (*) = 



= X(p (x) — 



x n "| 

X. sin a: $ q> (0 cos / ctt + X cos x J (p (/) sin / dt 
X J 
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The expression in the square brackets is equal to <p (x) 
so that 

<f>" (x) = Xy (x) — <p (x) 

From (1) and (2) we find that 

<P(ji) = 0, (p'(0) = 

Thus, the given integral equation reduces to the following 
boundary-value problem: 

q>" (x)-(X- l)(p(%) = 0, (3) 
<p(n) = 0, <r/(0) = (4) 

The three following cases are possible: 

(1) X — 1=0 or X=\. Equation (3) takes the form 
(p"(x) — 0. Its general solution will be <p (x) =CjX + C 2 . 
Utilizing the boundary conditions (4), we obtain (for find- 
ing the unknowns C t and C 2 ) the system 

\ Q = o 

which has a unique solution: C 1 — 0, C 2 = 0, and hence the 
integral equation has only the trivial solution 

<p (x) s= 

(2) X — 1 > or X > 1. The general solution of equation 
(3) is of the form 

<p (x) = C 1 cosh Y%— \x + C 2 sinh VX—lx 

whence 

<p' ( X ) = VX^-i. (C t sinh Vx^lx + C 2 cosh Vx~^lx) 

For finding the values of C, and C 2 , the boundary condi- 
tions yield the system 

/ C t cosh it K^— 1 -f-C 2 sinh:rcjA— 1 =0, 

\ c„ = o. 

The system has a unique solution: C l = 0, C 2 = 0. The in- 
tegral equation has the trivial solution (p (x) = 0. Thus, 
for X^l the integral equation has no characteristic num- 
bers and, hence, no eigenf unctions. 
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(3) X — 1 < or X < 1. The general solution of equation 
(3) is ~ 

<p(x) = C x cos]/l— Xx-\-C i sin]/"l— Ax 
Whence we find 

q>' (x) = Kl— J- (— C t sin X x + C 2 cos KIT^X x) 



In this case, for finding C x and C 2 the boundary conditions 
(4) yield the system 

C, cos ji j/l— ^ + C 2 sin ji V \—% = 0, 



-XC 2 = | 



(5) 



The determinant of this system is 

cosjij/l — X, sinnj/l — X 
KT^I 



Setting it equal to zero, we get an equation for finding the 
characteristic numbers: 



cosnjA — X sinnj/l — X 

VT=x 



= 



(6) 



or V\ — AccsjiKl — X — 0. By assumption V\ — X^=Q and 
so cosnl/l — X = 0. Whence we find that nj/l— X = 
= -2- + nn, where n is any integer. All the roots of equa- 
tion (6) are given by the formula 

For values X = X n the system (5) takes the form 

| 0, = 0, 



C 2 =0 

It has an infinite number of nonzero solutions 

C, = C, 

'2 _ 



I c x = c, 
I c, = o 
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where C is an arbitrary constant. Hence, the original inte- 
gral equation also has an infinity of solutions of the form 



<p(x) = C cos + 

which are eigenfunctions of this equation. 

Hence, the characteristic numbers and eigenfunctions of 
the given integral equation will be 

k= i— ( n +\y< <p n (*)= c ° s ( n +y) x 

where n is any integer. 

Find the characteristic numbers and eigenfunctions of the 
homogeneous integral equations if their kernels are of the 
following form: 

:(t+l), 



210. *(,,/)-{ H * +1) ' ° <X<t ' 

211 K(x t) = l {X+1) «- 2) > °< x < { > 
' \ (t+l)(x-2), 

(sin x cost, 0^.x^i, 
sin t cos X, y. 

213. K(x,t) = { S[nxc0st > ^ x < f ' 

I sin t cos ^ ^ x ^ it. 

214. K{x,t)={ s ! n f ir \ ( '-!>- 

I sin* sin (a:— 1), *<x<jt. 

(sin (.* + -£-) sin ( /— ) , 0<x< /, 
sin ( ^ + J sin ( x — -j) , tt^x^n. 

216. K(x, = e-l*-", 0<x< 1, 0</< 1. 

217. K( X ,t)={ - e :l sinh w °,$ x it' 

7 I — e *sinh*, 
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218. Show that if k lt X l z f (zX 2 are characteristic num- 
bers of the kernel K(x, t), then the eigenfunctions of the 
equations 

b 

(p(x) — % l \ j K(x, t)<p(t)dt = 

a 
b 

i|) (x)—X 2 J K (t, x) i|j (/) dt = 

a 

are orlhogonal, i.e., 

b 

^ <p (x) i|> (x) dx = 

a 

219. Show that if K(x, t) is a symmetric kernel, then 
the second iterated kernel /C 2 (a:, t) has only positive cha- 
racteristic numbers. 

220. Show that if the kernel K(x, t) is skew-symmetric, 
that is, K(t, x) = — K(x, t), then all its characteristic num- 
bers are pure imaginaries. 

221. If the kernel K(x, t) is symmetric, then 

£^ = 4. (m = 2, 3, ...) 

n= 1 

where X.„ are characteristic numbers and A m are the mth 
traces of the kernel K(x, t). 

Taking advantage of the results of Problems 209, 212, 
216, find the sums of the series: 

n=l 

( b ) 2- (4n 2 — I) 2 

n=l 

"L 1 

(c) y 1 ^ ^ > where (x„ are the roots of the equation 

n- 1 " 

2 cot n = (i — -. 
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The resolvent kernel of a symmetric kernel is a mero- 
morphic function of X, for which the characteristic numbers 
of the integral equation are simple poles. Its residues with 
respect to the poles X,- yield eigenfunctions (for any value 
of t) corresponding to these values of A,,-. 

Find the eigenfunctions of the integral equations who.e 
resolvent kernels are defined by the following formulas: 

222. R(x, t; A) = 3-X + 3 (1-^-1) ( 2< -l) 

116. K (X, t, K) 4X»— 16X+1S ' 

224. R(x, t; X) = 

4(5— 2A,) [3— 2X + (3— 6A,)*f] + 5 (4X-*— 8A.+3) (3x 2 — 1) (3/ a — 1) 
~ 4(1— 2A.)(3— 2X)(5— 2\) 

Fredholm integral equations with difference kernels. 
Suppose we have the integral equation 

71 

<p(*) = A, \ K(x— t)<p(t)dt (1) 

where the kernel K(x) ( — n^x^K) is an even function 
which is periodically extended to the entire x-axis so that 

K(x—t) = K(t—x) (2) 
It can be shown that the eigenfunctions of equation (1) are 
q>™(x) = cosnx (n=l, 2, ...), 
<P« ' (x ) = sin n* (n = 1 , 2, . . . ) 

and the corresponding characteristic numbers are 

= i (' l = 1 ' 2 >---) ( 4 > 
where a„ are the Fourier coefficients of the function K (x): 

n 

a„ = J K (x) cos nxdx (n — 1 , 2, . . . ) (5) 

Thus, to every value of %„ there correspond two linearly 
independent eigenfunctions cosnjc, sin/zx so that each k n 



} (3) 
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is a double characteristic number. The function <p (x) = 1 
is also an eigenf unction of equation (1) corresponding to 
the characteristic number 

at 

x °=i' where a °=ir 

We shall now show that, for example, cosnx is an eigen- 
f unction of the integral equation 



«P = §K(x-t)<p(t)dt (6) 



where 



a n = J /C (x) cos nx dx 
-ji 

Making the substitution x — / = 2, we find 

Jl JC-Jl 

^ K (x—t) cos nt dt = — ^ K (z)cosn(x— z)dz = 

-Jl K + JI 

JC + Jl AT + Jl 

= cosnx J K(z)cosnzdz + sin nx J /( (z) sin nzdz = :rta n cosnx 



since by virtue of the evenness of K (x) the second integral 
is zero, and the first integral is a Fourier coefficient a n 
multiplied by n in the expansion of the even function K (x). 
Thus, 



JL 

cos rtx = — !— [ K (x — /)cosn/ 



dt 



and this means that cosnx is an eigenfunction of equa- 
tion (6). 

Similarly, we establish the fact that sinnx is an eigen- 
function of the integral equation (6) corresponding to the 

same characteristic number — — . 
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225. Find the eigenfunction and the corresponding cha- 
racteristic numbers of the equation 

<p(x) = k jj cos 2 (a:— t) (p(t)dt. 
■ -jt 

226. Show that the symmetric kernel 

*(*. < )=aT 1 _ 2ftcos(x _ /)+A t (-»<*, t<n) 

has for |/i|< 1 the eigenfunctions 1, s'mnx, cos fix, which 
correspond to the characteristic numbers 1, l/h", l/h". 

227. Find the characteristic numbers and eigenfunctions 
of the integral equation 

it 

<p(x) = A, ^ K(x — t)q>(t)dt 
-n 

where K (x) = x 2 ( — jt < x ^1 it) is a periodic function with 
period 2n. 

Extremal properties of characteristic numbers and eigen- 
functions. 

The absolute value of the double integral (Hilbert's 
integral) 

|(/Cq>, <P)| = |SS*(*. t)<f>(x)<?(t)dxdt\ (1) 

a a 

where K(x, t) = K(t, x) is a symmetric kernel of some inte- 
gral equation, on the set of normalized functions <p(x), 
i. e.,. such that 

b 

(q>, q>) = J (p 2 (X) dx = 1 
a 

has a maximum equal to 

max|(/Ccp. <P)l = -jT7f ( 2 ) 

where X. t is the least (in absolute value) characteristic number 
of the kernel K(x, t). The maximum is attained forcp(x) = 
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=(p 1 (jc), which is the eigenfunction of the kernel correspon- 
ding to X 1 . 
Example. Find the maximum of 

n n 

l(Kq>, «p)| = |SS^(*. t)<p(x)<p(t)dxdt\ 



provided 

n 

(<P> <P) = $ <P 2 (x)dx=\ 

o 

if 

K (x, t) = cos x cos 2< + cos t cos 2*+ 1 
Solution. Solving the homogeneous integral equation 

<p (a;) = X ^ (cos ^ cos 2* + cos / cos 2x + 1 ) cp (^) 



as an equation with a degenerate kernel, we find the cha- 

1 2 
racteristic numbers ^j=— and A- 2i 3 = ±— and the corres- 
ponding eigenf unctions cp, (x) = C,, cp s (jc) = C 2 (cos x + cos 2x), 
(p 3 (x) = C 3 (cosx— cos 2x), where C lf C, and C 3 are arbitrary 
constants. 

The smallest (in absolute value) characteristic number is 

X 1 = -jp to which corresponds the eigenfunction (p,(x) = C,. 

From the normalization condition (q), <p)=l, we find C, = 

=± ,1_ . Hence 
\f2n 

n n 

max | J J ( cos cos 2^ + cos < cos 2x + 1 ) <p (/) d/ ( = 2n 



and it is attained on the functions (p(x) = ± 1 



228. Find the maximum of 



b b 



| \\K(x, t)q>(x)y(t)dxdt 
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b 



provided that \ cp 2 (x) dx = 1 , if 



a 



(a) K(x, t) = xt, 0<jc, *<1; 

(b) K(x, t) = xt + xH\ — 




Bifurcation points. Suppose we have a nonlinear integral 
equation 



By analogy with linear integral equations, the nonzero so- 
lutions (p(x)0O of equation (1) are called eigenf unctions 
and the corresponding values of the parameter X are called 
characteristic numbers of the equation. 

Ordinarily, the integral equations (1) do not have nonzero 
small solutions for small \X\\ that is, for small | X | equation 
( 1 ) has no eigenf unctions with small norm. Small eigenfunc- 
tions can appear in the case of increasing \X\. Let us intro- 
duce the following concept. 

The number X t is called a bifurcation point of the non- 
linear equation (1) if for any e>0 there is a characteristic 
number X of equation (1) such that \X — X |<e, and to 
this characteristic number there corresponds at least one 
eigenfunction (p (x) (tp (x) 0) with norm less than e: lKp||<e. 
Roughly speaking, a bifurcation point is that value of the 
parameter X, in the neighbourhood of which a zero solution 
of equation (1) branches; i.e., there appear small (in norm) 
nonzero solutions of equation (1). For linear problems, the 
bifurcation values coincide with the characteristic numbers. 

In problems of technology and physics involving condi- 
tions of stability, bifurcation points determine critical forces. 
Thus, the problem of the bending of a rectilinear rod of 
unit length and variable rigidity p (x) under the action of 



b 




(1) 



a 



Let (f(x) 



be a solution of the equation, and 
K(x, t, 0)^0 
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a force P leads to the solution of the following nonlinear 
integral equation: 



<p(jc) = Pp(*) S 



\k(x, t)<p(t) yi 



iK'Ax, t)y{t)di 



a 

dt 



(!') 

where cp(%) is the unknown function. 

For small P, equation (1') has a unique zero solution in 
the space C [0, 1] (by virtue of the contraction- mapping 
principle). This means that for small P the rod does not 
bend. However, a deflection occurs for forces greater than 
the so-called critical force of Euler. Euler's critical force 
is the bifurcation value. 

By way of an illustration of finding bifurcation points, 
let us consider the nonlinear equation 



Put 



Then 



q>(x) = A.$[(p(0 + <P 3 (0]* (2) 



1 



<p(*) = CX (3) 
and equation (2) reduces to the algebraic equation 

C = XC + k 3 C 3 (4) 

From (4) we get 
whence, by (3), 



^=0, c 2iS = ± 



Thus, for any < k < 1, equation ,(2) admits real nonzero solu- 
tions. For X = 1 it has only the zero solution q> = (three- 
fold). 
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Thus, for any 0<e< 1, the number X=l — e is a cha- 
racteristic number of equation (2) to which there correspond 
two eigenfunctions: 

Y l—e y 1 — e 

where e=l— X. Hence, the point a. = 1 is a bifurcation 
point of equation (2). One can also speak of bifurcation 
points of nonzero solutions of nonlinear integral equations. 

Find the bifurcation points of the zero solutions of the 
integral equations: 
■ 

229. <p(x) = X §xt(<p(t) + <( 3 (t))dt. 



1 

230. <p (x) = X J (3x — 2) t (q> (/) + q> 3 (/)) d/. 



(Bifurcation points are discussed in more detail in [10], 
[9], [28].) 



17. Solution of Homogeneous Integral Equations 
with Degenerate Kernel 

The homogeneous integral equation with degenerate kernel 




where the parameter X is not its characteristic number (i. e., 
A (X) 0) has a unique zero solution: (p(x) = 0. But if X 
is a characteristic number (A (X) = 0), then, besides the zero 
solution, equation (1) also has nonzero solutions — the eigen- 
functions which correspond to that characteristic number. 
The general solution of the homogeneous equation (1) is 
obtained as a linear combination of these eigenfunctions. 
Example. Solve the equation 

<p (x) — X J (cos* x cos 2t + cos 3 1 cos 3x) <p (/) dt = 
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Solution. The characteristic numbers of this equation are 
'K 1 =-^, % i = ^-\ the corresponding eigenf unctions are of the 
form 

tpj (x) = cos 2 x, <p 2 (x) = cos 3x 
The general' solution of the equation is 
q> (x) = C cos 2 x if k =-i, 
q>(x) = Ccos3x if *.=-|' 
<p(x) = if X a-=^- 

where C is an arbitrary constant. The last zero solution is 
obtained from the general solutions for C = 0. 

Solve the following homogeneous integral equations: 
ji 

231. (p(x) — K I cos (x+t)<p(t) dt = Q. 



1 

232. (f>(x)— X $ arc cos xy (t) dt = 0. 



n/t 

233. «PW-2 JrJ&d/^O. 



2 

234. <f>{x)— J- §\x\<p(t)dt = 0. 

-2 
1 

235. <p(x) + 6$(;t a — 2xt)<p\t)dt = 0. 



18. Nonhomogeneous Symmetric Equations 

The nonhomogeneous Fredholm integral equation of the 
second kind 

b 

<f>(x)-X^K(x, t)y(t)dt = f(x) (1) 
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is called symmetric if its kernel K(x, t) is symmetric: 
K(x, t) = K(t, x). 

If f(x) is continuous and the parameter X does not co- 
incide with the characteristic numbers X„ (n = 1 , 2, . . . ) of 
the corresponding homogeneous integral equation 

b 

<p(x)—X^K(x,t)<f>(t)dt = (2) 

a 

then equation (1) has a unique continuous solution, which 
is given by the formula 

00 

<p(x) = f{x)-kH x=T n <t> n (x) (3) 
where <f n (x) are eigenfunctions of equation (2), 

b 

a n =lf(x)(p n (x)dx (4) 

a 

The series on the right side of formula (3) converges ab- 
solutely and uniformly in the square as^.x, t^.b. 

But if the parameter X coincides with one of the cha- 
racteristic numbers, say X = X k , of index q (multiplicity of 
the number X k ), then equation (1) will not, generally 
speaking, have any solutions. Solutions exist if and only 
if the q conditions are fulfilled: 

b 

(/, (pj = or \f(x)<p m (x)dx = (5) 

a 

(m=l, 2, q) 

that is, if the function f{x) is orthogonal to all eigenfun- 
ctions belongirg to the characteristic number X k . In this 
case equation (1) has an infinity of solutions which con- 
tain q arbitrary constants and are given by the formula 

00 

<((x) = f(x)-X £ r zx<P„W + 

n=q+i " 

+ Ci<P» W + C 2 cp 8 {x)+...+ C g <p g (x) (6) 
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where Cj, C 2 , . .., C q are arbitrary constants. 
In the case of the degenerate kernel 

m 

K(x, 0=2 a* WMO 

formulas (3) and (6) will contain finite sums in place of 
series in their right-hand members. 

When the right-hand side of equation (1), i. e., the fun- 
ction f(x), is orthogonal to all eigenfunctions <p n (x) of 
equation (2), the function itself will be a solution of equ- 
ation (1): cp(x) = f(x). 

Example 1. Solve the equation 
i 

(p(xj—k^K(x,t)<p(t)dt=x, (1) 

o 

where 

x(t— 1), if 0<x<* 
< x < 1 



j x(t— 1), if 



Solution. The characteristic numbers and their associated 
eigenfunctions are of the form 

%„ = — n 2 n 2 , q>„ (x) <=> sin nnx, n=l, 2, ... 
If k^X n> then 

n= 1 



will be a solution of equation (1). We find the Fourier 
coefficients a n of the right side of the equation: 

i i 

a„ = \ x sin rmxdx =\xd[ = v — '- — 

" J. J V 7 n" 



Substituting into (2), we get 

, , X v ( — l) n+1 
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For X = — n 2 n 2 equation (1) has no solutions since 

|_nn + l 

a» = — ^0 

Example 2. Solve the equation 
i 

cp(x) — X ^ K (x, t) cp (t) dt = cos nx 



where 



/ (x+\)t, 0<x<f, 
K(x, /) = < 



(t+\)x, /<x< 1 

Solution. The characteristic numbers are 

X =l, X n = — nW (n=l, 2, ...) 

Their associated eigenf unctions are 

(p (x) = e x , cp„(x) = sin rcnx + mi cos nnx (n = 1, 2, ...) 

If and X=^= — n 2 n 2 , then the solution of the given 

equation will have the form 

<p (x) = cos nx — x\ + ^ \±nt n 2 ( sln n3i;A: + nit cos /?JlA; ) 
and since 



a. 



1 * j l+e 

e* cos nx ax = — , , „ 
1 + ji 2 



» I 0, n#l, 

a„ = \ cos nx (sin nnx + nn cos nnx) dx = < n _ , 





it follows that 



cp (x) = cos nx+X [1±£ ^__JL_ (s in nx + n cos nx)] 

For X.= l and X = — n a (n = l) the equation has no solu- 
tions since its right-hand side, that is, the function cosnx, 
is not orthogonal to the corresponding eigenfunctions 

<Po (*) = e*. 

fj (x) = sin nx + n cos nx 
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But if X = — n 2 jx 3 , where n = 2, 3, then the given 

equation has an infinity of solutions which are given by 
formula (6): 

q> (x) = cos nx+k \j^j^r x — 2 (X " n ^ (sinnx + Jtcosnx)] + 
+ C (sin nnx + tin cos nnx) 

where C is an arbitrary constant. 

In certain cases, a nonhomogeneous symmetric integral 
equation can be reduced to a nonhomogeneous boundary- 
value problem. This is possible when the kernel K(x, t) of 
the integral equation is a Green's function of some linear 
differential operator. Let us illustrate how this is done. 

Example 3. Solve the equation 
i 

, <f>(x)—X^K(x, t)<p(t)dt = e* (1) 



where 

sinh x sinh (t— 1) n -„^, 

sinhl ' 

sinh t sinh (x — 1) . ^ - . 
sinh , , 

Solution. Rewrite the equation as 

X 

tp (x) = e* + X si f n l) j sinh t y (t) dt + 



1 

+^^Jsinh(/-l)<p(/)<tt (2) 
Differentiating twice, we obtain 

X 

<p' (x) = e* + J sinh t<p (t) dt + 



K(x, /) = ■ 



A. sinh (x — 1) 
sinh 1 



sinh xy (x) + X sinh t 
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x 

<*> = e " + ^sinhT^ I Sinh ^ f > * + 
o 

1 

. Xsinhxf . , ,, /j\ jj i X cosh (*— 1) sinh x , . 

+ UnhT J Sinh ^- 1 >«P(0*+ khT fW" 

X cosh x . , . , . , . 

-ifinrr^ mh < JC - 1 ) < PW 

or 

q,"(je) = e* + (p(A;) + A.<p (*) 

Putting jc = and x=l in (2), we get q>(0) = l, cp (1) = e. 
The required function q>(x) is a solution of the nonhomo- 
geneous boundary-value problem 

q>» (jc) - (A, +.l)q> (*) = «* (3) 
<P(0)= 1, cp(l) = e (4) 

Let us consider the following cases: 

(1) X.+ 1=0, or X = — 1. Equation (3) is of the form 
(p"(x) = e x . Its general solution is 

q>(x) = C l x + C 2 + e x 

Taking into account the boundary conditions (4), we get 
the following system for finding the constants C, and C 2 : 

C,+ l = l. 

C^ + C. + e-e 

Its solution is of the form C t = 0, C 2 = 0, and, hence, 

9 (x) = e* 

(2) a.+ 1>0, or k> — 1, X =^0. The general solution 
of equation (3) is 

<p (x) = C 1 cosh VT+Tjc + C 2 sinh VT+lx — j 

The boundary conditions (4) yield the following system for 
finding C, and C 2 : 

C -1-1 

C t cosh + C t sinh /l +1 — j- = e 
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whence 

C e-™*yj+t (l + ') 
1 2 sinhK 1 + X, V / 

After simple manipulations, the unknown function q> (x) is 
reduced to 

. . /, , 1 \sinh/T+~Ml--*) e x 



(5) 



sinhl^l+X 

(3) X+l <0, or A,< — 1. Denote X + l= — n 2 . The ge- 
neral solution of equation (3) is 

q> (x) = C x cos u-x + C 2 sin |xx + j-p^- 
The boundary conditions (4) yield the system 

1 ' 1+|X 2 

C x cos u. + C 2 sin u. = e y-^j- 

In turn, two cases are possible here: 

(a) (i is not a root of the equation sin |i = 0. 
Then 

/-> __Ml_ n _ (e — cos n) n 2 
1 — l + |x 2 ' 2 — (l+p. 2 )sinn 

and, hence, 

/ \ r i e — cos u . "I . e x 

<p (*) =rfr [ c os ^ + sin H + t+f 

where n = K — a. — 1- 

(b) (A is a root of the equation sin u-=0, i. e., 
u. = wi(n=l, 2, ...)• System (5) is inconsistent and, con- 
sequently, the given equation (1) has no solutions. 

In this case, the corresponding homogeneous integral 
equation 

(f(x) + (l+n 2 n 2 ) < \K(x, t)y(t)dt = (6) 
o 

will have an infinity of nontrivial solutions, that is, the 
numbers ^„ = — (1 + n 2 n 2 ) are characteristic numbers and 
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their associated solutions <p„ (x) = sin nnx are eigenfunctions 
of equation (6). 

Solve the following nonhomogeneous symmetric integral 
equations: 

236. «p(x)— £ §K(x, t)y(t)dt = ±, 

*2f2L,o <x< t, 




sinh.* sinh (t — 1) ~ ^ ~ . 
sTnhl ' ^ XS5 = * 



sinh / sinh (x— 1) ,' . 



K(x, /) = 

sinh 1 

238. <p(x)— I lK(x, /)q>(/)d/ = x— 1, 



o 

x — /, ^ x< /, 



K(x, /) = , 

/ — X, t ^ 1. 

Jt/2 

239. q>(x)— 2 J K(x, t) cp(/)d/ = cos 2x, 



K(x, t) - 



sin x cos/, <;£</, 



sin/ cos /^X^y. 
240. <p(x) — X J/C(x, /)cp(/)d/=l, 



K(x, /) = / sin * 
[ sin / 



x cos/, 0^x</, 
cos / ^ a; ^ n. 







241. <p(x) — ^ J /C (at, t)y(t)dt = x, 

(x+[)(t-3) 
(t + l)(x— 3), * <x< 1 



, )== | (x+l)(/-3), 0<*</, 
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242. y{x)— \K(x, t) cp (t) dt = sin x, 

o 

| sin (x+x) sin (* — -J-).-0<*</, 
K(x, t) = \ i . , , 

y sin ( t+ -j J sin ( x — -j), t ^ x ^.n. 

i 

243. (p(x)— i \K(x, t)(f> (t)dt =--s'mhx, 

o 

( — e"'sinhx, O^x^/, 



244. cp (x) — X J K (x, t) cp (/) dt = cosh 
o 

cosh x cosh (t — 1) 



K(x, t) = 



sinh 1 



cosh t cosh (x — 1 ) , . . . 
sinh 1 - ^ 

245. q>(jc) — A. J |jc — /|cp(/)d/ = 1. 



19. Fredholm Alternative 

For Fredholm integral equations we have the theorems: 
Theorem 1 (Fredholm alternative). Either the nonhomo- 
geneous linear equation of the second kind 

b 

q>(x)—%\K(x, t)<p(t)dt=f(x) (1) 

a 

has a unique solution for any function f(x) (in some suf- 
ficiently broad class) or the corresponding homogeneous equa- 
tion 

b 

(p(x)—K i \ ) K(x, t)<p(t) dt = (2) 

a 

has at least one nontrivial (that is, not identically zero) 
solution. 
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Theorem 2. // the first alternative holds true for equa- 
tion (1), then it holds true for the associated equation 

b 

*(x)-klK(t, x)y(t)dt=g(x) (3) 

a 

as well. The homogeneous integral equation (2) and its 
associated equation 

b 

Kf(x) — %\K{t, x)ty(t)dt = Q (4) 

a 

have one and the same finite number of linearly independent 
solutions. 

Note. If the functions ^{x), q> 2 (x), <p„(x)are solu- 
tions of the homogeneous equation (2), then their linear 
combination 

n 

<P (*) = C i<Pi (x) + C 2 cp 2 (*)+...+ C„<p„ (x) = 2 C k <Vk (x) 

where the C k {k=\, 2 n) are arbitrary constants, is 

also a solution of the equation. 

Theorem 3. A necessary and sufficient condition for the 
existence of a solution q>(x) of the nonhomogeneous equa- 
tion (1) in the latter case of the alternative is the condition 
of orthogonality of the right side of the equation, i.e., 
of the function f(x), to any solution Mp(x) of the homoge- 
neous equation (4) associated with (2): 

b 

lf(x)ty{x)dx = (5) 

a 

Note. When condition (5) is fulfilled, equation (1) will 
have an infinite number of solutions, since this equation 

will be satisfied by any function of the form cp (x) + <p (x), where 
q>(x) is some solution of equation (1) and y(x) is any 
solution of the corresponding homogeneous equation (2). 
Besides, if equation (1) is satisfied by the functions ^(x) 
and (jpj(jc), then by virtue of the linearity of the equation 
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their difference, y l (x) — (p 2 (x), is a solution of the corres- 
ponding homogeneous equation (2). 

The Fredholm alternative is particularly important in 
practical situations. Instead of proving that a given integral 
equation (1) has a solution, it is often simpler to prove 
that the appropriate homogeneous equation (2) or its asso- 
ciated equation (4) has only trivial solutions. Whence it 
follows, by virtue of the alternative, that equation (1) 
indeed has a solution. 

Remarks. (1) If the kernel K(x, t) of the integral equa- 
tion (1) is symmetric, that is, K(x, t)^=K(t, x), then 
the associated homogeneous equation (4) coincides with 
the homogeneous equation (2) which corresponds to equa- 
tion (1). 

(2) In the case of the nonhomogenebus integral equation 
with degenerate kernel 



the orthogonality condition (5) of the right side of this 
equation yields n equalities 



br n 



a L* = i 



b 



[f(t)b k (t)dt = (fc=l, 2, . 



, n) 



a 



Example 1. 



<P(jc)— I J (5jc 2 — 3)P<p(t)dt 



o 



Solution. We have 



q>{x) = Ck(5x i — 3) + e x 



(1) 



where 




(2) 







Substituting (2) into (1), we get 



C = CI $ (5* 4 — 3/ 2 ) dt + J dt 







5-1583 
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whence 

C = e—2 

For any X, the given equation has a unique solution: 

<p(x) = k (e— 2) (5x 2 — 3) + e* 

and the corresponding homogeneous equation 

i 

<p(*)— % $(5x 2 — 3)t\(t)dt = 



has a unique zero solution <p(x) = 0. 
Example 2. 

i 

cp(jc) — % J sin lnx<p(0d^ = 2x 



Solution. We have 

<p (x) = Ck sin In a; + 2x 

i 

where C = $ <p (/) Substituting the expression 9 (i) into 


the integral, we obtain 

C = CK J sin In t At + 1 


whence 

c(i+A) = , 

If — 2, then the given equation has a unique solution 

2X 

V W = 2+X s ' n ' njt + 2Ar; the corresponding homogeneous 
equation 

<p(x) — a. J sin In xq> (t) it = 



has only the zero solution <p (x) = 0. 
But if K=— 2, then the given equation does not have 
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any solutions since the right side / (x) = 2x is not orthogonal 
to the function sin In x; the homogeneous equation has an 
infinity of solutions since it follows from the equation defi- 
ning C, 0C = 0, that C is an arbitrary constant; all these 
solutions are given by the formula 

q)(x) = C sinlnx (C=— 2C) 

Example 3. 

JT 

<p (x)~ X J cos (x+ t) q> (t) dt = cos 3x 

o 

Solution. Rewrite the equation in the form 

n 

(f(x) — k ^ (cos x cost — sin x sin t) <p (t) dt — cos 3* 

o 

Whence we have 

<p (x) = cos x— Cjk sin * -(-cos 3.x; (1) 

where 

n 

Cj = J <p (t) cos / dt, 

JT 

C 2 = ^ cp (t) sin tdt 



Substituting (1) into (2), we get 

C 1 = \ (Cjk cos t— C 2 l sin t + cos 3/) cos / dt, 



Jl 

C„ = J (CjK cos * — C 2 a. sin t + cos 3*) sin * d< 



5* 
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C, f 1 — k J cos 2 / dt) 4- C,X I sin t cos t dt = 

n 

= ^ COS 3/ COS ^ d*. 
o 

— C,?,$ cos / sin f dt + C 2 f 1 + X j sin 2 / dA = 

\ u / 

71 

= ^ cos 3/ sin / dt 



c,(i-x£) = o. 

C 5 (l+Xf)=0 
The determinant of this system is 



(3) 



A(X) 



1— X- 







= 1 -^A 2 
4 



(1) If a- =5*= ±-£(A (A.)#0), then system (3) has a unique 

solution C x = 0, C 2 = and, hence, the given equation has 
the unique solution cp (x) = cos 3* and the corresponding 
homogeneous equation 



<p(x)— X I cos (x + t)q>(t)dt = 



(4) 



only has the zero solution (p(jt) = 0. 

2 

(2) If X = — , then system (3) takes the form 

lC,.0 = 0, 
10,. 2 = 
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Whence it follows that C 2 = and C 1 = C, where C is an 
arbitrary constant. The given equation will have an infinity 
of solutions which are given by the formula 

2 

<j> (x) = — C- cos x + cos 3x 

or 

(p (x) = C • cos x -f- cos 3x = ; 

the corresponding homogeneous equation (4) has an infinity 
of solutions: 

(f)(x) = CcOSX 

(3) If %=——, then system (3) takes the form 
t C 2 = 

whence C x = 0, C 2 — C, where C is an arbitrary constant. 
The general solution of the given equation is of the form 

2 

<p (x) = — C • sin x + cos 3x 

or 

<p(x) = C-sinx + cos3x ^C = ^^ 

In this example, the kernel K{x, t) = cos(x + t) of the 
given equation is symmetric: K(x, t) — K(t, x); the right 
side of the equation [that is, the function / (x) = cos 3x] is 
orthogonal to the functions cosx and sinx on the inter- 
val [0, n]. 

Investigate for solvability the following integral equations 
(for different values of the parameter A.): 

246. cp(x)— l^cos 2 xy(t)dt = I. 



1 

247. cp(x) — K J xe'q> (t) dt = x. 

- 1 

in 

248. <p(x)— X \ \x— n\(p(t)dt=x. 
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1 

249. <p{x)— X J (2x/ — 4x*)<((t)dt= 1 — 2x. 



I 

250. <p(x)— X J (jc 2 — 2xt)<p(t)dt = x 3 — x. 

-i 

251. <f{x)~ cos x cos t + -^- sin 2x sin 2t^j (p(t)dt = 



a 

■■ sin x. 
i 



252. q>(x) — X 0<P(0<ft-l 



where 

j coshxsinh/, O^x^/, 
K(*. ^ = \ CO sh^-sinhx, 

20. Construction of Green's Function 
for Ordinary Differential Equations 

Suppose we have a differential equation of order n: 

L [y] - P. (*) «/" + ft W f/"" 1 +.,.+p.(x)y = (1) 

where the functions p t (x), p 1 (x), p n (x) are continuous 
on [a, b], p {x) =^= on [a, b], and the boundary conditions are 

V* (y) = «*y (a) + ail/' (a) + ■ ■ • + «2" V 1 ^) + 

+ p*p + Pbr (ft) + . . . + pr V 1 W 

(fe=l. 2 n) (2) 

where the linear forms V x , . . . , V„ in y (a), y' (a), . . . , t/ (n_1) (a), 
t/(6) </ ( " -1) (ft) are linearly independent. 

We assume that the homogeneous boundary-value problem 
(l)-(2) has only a trivial solution y(x) = 0. 

Definition. Green's function of the boundary-value pro- 
blem (l)-(2) is the function G(x, |) constructed for any 
point I, a < | < b, and having the following four properties: 

(1) G(x, |) is continuous and has continuous derivatives 
with respect to x up to order (n — 2) inclusive for a<je<6. 
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(2) Its (n — l)th derivative with respect to x at the point 
x = % has a discontinuity of the first kind, the jump being 

equal to -4-; , i.e., 

a»-»G(x, l) \ d"~ 1 G (x, l) \ 1 

dx"-i 1=1+0 dx»~i \ x =t-o Po(l) 1 ; 

(3) In each of the intervals [a, £) and (|, b] the function 
G(x, |), considered as a function of x, is a solution of 
equation (1): 

L[G] = (4) 

(4) G(a;, I) satisfies the boundary conditions (2): 

V k (G) = (*=1, 2 n) (5) 

Theorem 1. // //te boundary-value problem (l)-(2) has 
only the trivial solution y(x) = 0, then the operator L has 
one and only one Green's function G(x, £). 

Proof. Let y i (x),y i (x), . . . , y„ (x) be linearly independent 
solutions of the equation L [y] — 0. Then, by virtue of Pro- 
perty (3), the unknown function G(x, £) must have the 
following representation on the intervals [a, t) and (|, b]: 

G (x, |) = a x y t (x) + a 2 y t (x)+...+ a n y n (x) f or a < x < \ 
and 

G (x, I) = (x) + b^ (*)+...+ b n y n {x) for % <x < b 

Here, a lt a 2 a„, b lt b t b„ are some functions of £. 

The continuity of the function G(x, £) and of its first n— 2 
derivatives with respect to x at the point x = \ yields the 
relations 

[b iyi (1) + . . . + b n y n (I)] - [a iyi (1)+...+ a n y n (£)] = 0, 
M (£)+••• + bjn (I)] - Mi (!)+•••+ a n y n {\)\ = 0, 



[W 2 ' (!)+••• +«- s> (i)]-k</r 2) (£) + 
...+a n yr a) a)]=o 

and condition (3) takes the form 

[Mr 1 ' (!)+•■• + W" (i)] -N"-" (y + 
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Let us put c k (t) = b k (l)-a k (t) (k=\, 2, 
get a system of linear equations in c ft (|): 

Cii/i (I) + c 2 y 2 (£)+...+ c n y n (1) = 0, 
c$'x (i) + c 2 y' 2 (|) + • • • + c n y' n (I) = 0, 



c 1 <- 2) (I) + c 2 y 2 n ~» (I) + . . . + c n <-*> (E) = 0, 

c^r 11 a) +c t yr u (!)+••• (5) = 



, n); then we 



Po(l) 



(6) 



+ oj»-V-»(a), 
+ pl»-»y»-»(&). 



The determinant of system (6) is equal to the value of the 
Wronskian W (y lt y 2 y n ) at the point x = £ and is there- 
fore different from zero. For this reason, system (6) uniquely 
defines the functions c k (1-)(k= ■ 1, 2, ...,n). To determine 
the functions a k (%) and let us take advantage of the 

boundary conditions (2). We write V k (y) in the form 

V k (y) = A k (y) + B k (y) (7) 

where 

A k {y)=a l jf(a) + aj t li y {a)+ . 
B k (y) = $ k y(b) + W>y'(b)+. 

Then, by conditions (5), we get 

V k (G) = a x A k (yj + a 2 A k (y 2 ) +...+ a„A k (y„) + 

+ hB k (yj + b 2 B k (y 2 ) + ...+ b n B k (y n ) = (k = 1 , 2 n) 

Taking into consideration that a k = b k — c k , we will have 
(&i —Ci) A k (yj + (b 2 -c 2 ) A k (y 2 ) +...+(b n -c n ) A k (y n ) + 
+ hB k (yj + b 2 B k (y 2 ) +...+ b n B k (y n ) = 
(k=l, 2, ...,n) 

Whence, by virtue of (7), 
KV k (yj + b 2 V k (y 2 ) +...+ b„V k {y n ) = 
= fAW + fAte)T...+fAW(*='. 2, ...,n) (8) 
Note that system (8) is linear in the quantities b lt .... b n . 
The determinant of the system is different from zero: 

VM VM ••• vm 

VM V 2 (y 2 ) ... V t (y a ) 



VM v n (y 2 ) ... v n (y n ) 



¥=0 



(9) 
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by virtue of our assumption concerning the linear indepen- 
dence of the forms V lt V 2 V n . 

Consequently, the system of equations (8) has a unique 
solution in b l (£), b s (l), b n {%), and since a k (l) = b k (l) — 
— c ft (|), it follows that the quantities a k (l)(k= 1, 2, n) 
are defined uniquely. Thus the existence and uniqueness of 
Green's function G(x, £) have been proved and a method 
has been given for constructing the function. 

Note 1. If the boundary-value problem (l)-(2) is self- 
adjoint, then Green's function is symmetric, i.e., 

G(x, l) = G(l, x) 

The converse is true as well. 

For the conditions of self-adjointness of the boundary- 
value problem for differential operators of the second and 
fourth orders see [20], Vol. I. 

Note 2. If at one of the extremities of an interval [a, b] 
the coefficient of the highest derivative vanishes, for example 
p (a) = 0, then the natural boundary condition for bounded- 
ness of the solution at x — a is imposed, and at the other 
extremity the ordinary boundary condition is specified (see 
Example 2 below). 



An Important Special Case 

Let us consider the construction of Green's function for 
a second-order differential equation of the form 

(p (x) y')' + q (x)y = 0, 
p (x) =^ on [a, b], p(x)£C n) [a, b] (10) 

with boundary conditions 

y(a) = y(b) = (11) 

Suppose that t/xC*:) is a solution of equation (10) defined 
by the initial conditions - 

yx(a) = 0, y[(a) = a^=0 (12) 

Generally speaking, this solution need not necessarily satisfy 
the second boundary condition; we will therefore assume 
that y, (6)^=0. But functions of the form C 1 y 1 (x), where 
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C x is an arbitrary constant, are obviously solutions of equa- 
tion (10) and satisfy the boundary condition 

y(a) = 

Similarly, we find the nonzero solution y 2 (x) of equation 
(10), such that it should satisfy the second boundary con- 
dition, i. e., 

y,(b) = (13) 

This same condition will be satisfied by all solutions of the 
family C 2 */ 2 (x), where C 2 is an arbitrary constant. 

We now seek Green's function for the problem (10)-(11) 
in the form 

G ^-{ C 2 y A x) for£<*<6 < 14 > 

and we shall choose the constants C 1 and C 2 so that the 
Properties (1) and (2) are fulfilled, i. e., so that the function 
G(x, |) is continuous in x for fixed £, in particular, conti- 
nuous at the point x = %: 

c 1 y 1 (l) = c 2 y,a) 

and so that G' x (x, £) has a jump, at the point x = %, equal 
to — Lr : 

Rewrite the last two equalities as 

-C 1 y 1 (l) + C t y 2 (l) = 0, | 

-CJAQ + Cd'Ai)--^ j (15) 

The determinant of system (15) is the Wronskian W [y 1 (x), 
y 2 (x)] = W{x) computed at the point x = £ for linearly 
independent solutions y 1 {x) and y 2 (x) of equation (10), and, 
hence, it is different from zero: 

so that the quantities C 1 and C 2 of system (15) are deter- 
mined at once: 

r y * ® r - yi(^) 
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Substituting the expressions for C, and C % into (14), we 
finally get 

yi (*) yi (S) n< - v< "t 



j/i (I) Vi, (x) 



Note 1. The solutions y t (x) and y 2 (x) of equation (10) 
that we have chosen are linearly independent by virtue of 
the assumption that y 1 (b) 0. 

Indeed, all solutions linearly dependent on y 1 (x) have 
the form (x) and, consequently, for C t 0, do not 
vanish at the point x = b at which, according to our choice, 
the solution y 2 (x) vanishes. 

Note 2. The boundary-value problem for a second-order 
equation of the form 

y" (x) + Pi (x) y' (x) + p 2 (x) y(x) = (18) 

and boundary conditions 

y(a) = A, y{b) = B (19) 

reduces to the above-considered problem (10)-(11) as follows: 

(1) The linear equation (18) is reduced to (10) by multi- 
plying (18) by p(x) = e$ PlMdx [we have to take p(x)p 2 (x) 
for q(x)]. 

(2) The boundary conditions (19) reduce to zero conditions 
(11) by a linear change of variables 

z(x) = y {x)—jZA( x —a) — A 

The linearity of equation (18) is preserved in this change, 
but unlike equation (10), we now obtain the nonhomogeneous 
equation L[z]=f(x), where 

f(x) = - [A+^(x-a)]q(x)-^p(x) 

However, we construct Green's function for the homogeneous 
boundary- value problem L[z]=0, z(a) = z(b) = 0, which 
fully coincides with the Problem (10)-(11). 
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Example 1. Construct Green's function for the homo- 
geneous boundary-value problem 

</ IV (x) = 0, (1) 
*/(0) = </'(0) = 0, \ 

y(l) = */'(l) = I (Z > 

Solution 1. We shall first show that the boundary-value 
problem (l)-(2) has only a trivial solution. Indeed, the 
fundamental system of solutions for equation (1) is 

y 1 (x)=l, y 2 (x) = x, y 3 (x) = x\ </ 4 (x) = x 3 (3) 

so that its general solution is of the form 
y ( X ) = A + Bx + Cx 2 + Dx 3 

where A, B, C, D are arbitrary constants. The boundary 
conditions (2) give us four relations for determining A, B, 
C, D: 

t,(0) = A = 0, 
y> (0) = B = 0, 

y(l) = A + B + C + D = 0, 
«/'(l) = B + 2C + 3D = 

Whence we have A = B = C = D = 0. 

Thus, the problem (l)-(2) has only a zero solution 
y(x)=0, and, hence, for it we can construct a (unique) 
Green's function G {x, \). 

. . 2. We now construct Green's function. Utilizing the fun- 
damental system of solutions (3), represent the unknown 
Green's function G{x, |) in the form 

G (x, 1) = a t ■ 1 rf- a 2 ■ x + a s - Jc 2 + a 4 • x 3 for ^ x < |, (4) 
G(x, 5) = V1+V* + V* 2 + V* 3 for|<x<l (5) 
Where Oj, a 2 , a 3 , a 4 , b lt b 2 , b 3 , 6 4 are as yet unknown 
functions of £. Put c k (l) = b k (l)— a k {l){k= 1, 2, 3, 4) and 
write out the system of linear ; equations for finding the 
functions c k (%) [see system (6) on p. 136]: 

' ~ ~ c 1 + c£ + c 3 l* + c i % 3 = 0, | 
c 2 + c 3 -2i + c 4 .3i a = 0, I 

c a ;2-f6 V 6i-0, f , ~W 
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Solving the system, we get 

Cl (i) = -|i 3 , 



C 3 (l) = — C 4 (6) = rg- 



(7) 



We further take advantage of Property (4) of Green's fun- 
ction, namely, that it must satisfy the boundary conditions 
(2), i.e., 

G(0, |)=o, g;(0, ?) = o, i 
G(i,i) = o, g;(1,&)=o j 

In our case these relations take the form 

a, = 0, 



a 2 = 0, 
6 2 + 26 s + 36 4 = 



(8) 



Taking advantage of the fact that c k = b k —a k (k= 1, 2, 3, 4), 
we find from (7) and (8) that 



a, = 0; 6 t = a 2 = 0; 6, = -J-t 2 ; 



(9) 



G(jc, £) = 



Putting the values of the coefficients a lt -a t , : . . , b t from (9) 
into (4; and. (5), we obtain the desired Green's function: 

This expression is readily transformed to 

G (x, |) = [~x — x % + 1 r>) £■ - (1 -1 *■ + 1 g» 
for I ? 
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so that G(x, $) = G(£, x )> l - e > Green's function is sym- 
metric. This was evident from the start since the boundary- 
value problem (l)-(2) was self-adjoint. 

The reader is advised to establish this by himself. Also 
we suggest checking to see that the Green's function which 
we have found satisfies all the requirements (1 to 4) given 
in the definition. 

Example 2. Construct Green's function for the differential 
equation 

xy" + y' = (1) 

for the following conditions: 

y (x) bounded as x — »- 0, 

y(\) = a y'(l), a^O (2) 

Solution. First find the general solution of equation (1) 
and convince yourself that the conditions (2) are fulfilled 
only when 

y(x) = o 

Indeed, denoting y'(x) — z(x) we get xz' + z = 0, whence 
lnz = lnc, — \nx, z = -j- and, hence, 

y(x) = c l \nx + c 1 (3) 

It is clear that y{x) defined by formula (3) satisfies the con- 
ditions (2) only for c, = c, = 0, and, hence, Green's function 
can be constructed for the problem (l)-(2). 
Let us write down G(x, |) formally as 

i a x -f a t In x for < x <; i, 
G(x, £) = j bi+bt i nx for 1< x < 1 (4) 

From the continuity of G(x, I) for x = \ we obtain 

fti-r-hj ln|— a x — a t ln| = 
and the jump G' x (x, I) at the point * = £ is equal toy so that 

Putting 

= — a v c 2 = 6 4 — o, (5) 
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we will have 



whence 



( Cj + c 2 In | = 0, 
I * c, = l 

c x = — In I, c 2 = 1 (6) 



Now let us use conditions (2). The boundedness of G(x, g) 
as * — <• gives us a 2 = 0, and from the condition G (x, |) = 
= aG' x (x, |) we get b 1 = ab i . Taking into account (5) and 
(6), we get the values of all coefficients in (4): 

a ± = a -f In g, a i = 0, b t = a, b 2 = 1 

Thus 

a + lng, 0<*<£, 
a -f- In x, | ^ jc ^ 1 



G(AT, |) = 



Example 3. Find Green's function for the boundary-value 
problem 

y"(x) + k*y = 0, 
{/(0) = t/(l) = 

Solution. It is easy to see that the solution y t (x) — sin kx 
satisfies the boundary condition y 1 (0) = 0, and the solution 
y i (x)=sink(x — l)satisfies the condition y 2 (l)=0; they are 
linearly independent. Let us find the value of the Wronskian 
for sin ft* and sin ft (x — 1) at the point x = %: 

sinfti sin ft (|— 1) 
ftcosft£ ft cos ft (| — 1) 
= ft[sinft| cos ft (| — 1) — sin ft (1 — l)cos ft£] = ftsinft 

Noting, in addition, that in our example p(x)—\, we get, 
by (17), 

/ sin ft (6-1) sin <Jt<i, 
Civ %\— ) ft sin ft o. 
v{X, sinA g sinft ( x _ 1 ) 

\ ft sin £ ' ^ ^ 

In the following examples, establish whether a Green's 
function exists for the given boundary-value problem and 
if it does, construct it. 



144 INTEGRAL EQUATIONS 

253. y" = 0; y(0) = y'(l), r/' (0) = # ( 1 ). 

254. </" = 0; </(0) = </(l), y'(0) = y'(l). 

255. y" + y = 0; y(0) = y(n) = 0. 

256. #'v = o ; y(0) = y'(0) = y"(\) = y"'(\) = 0. 

257. */"' = 0; </(0) = y' (1) = 0; (0) = 

258. y"' = 0; y(Q) = y(\) = 0; t/' (0) = y'(l). 

259. y" = 0; t/(0) = 0, = y' (1). 

260. y" + y' = 0; f/(0) = </(l), */ (0) = </' (1). 

261. y" — k 2 y = 0(k^=0); y(0)=y(l) = 0. 

262. </" + </ = 0; «/(0) = z/(l), */'(0) = {/'(l). 

263. y"' = 0; y (0) = = 0, (0) + y' (1) = 0. 

264. y" = 0; (0) = % (0), y' (1) = - Hy{\). 

265. x 2 «/" -f 2xi/ = 0; y(x) is bounded for 

</(l) = a*/'(l). 

266: v + 6x 2 y' " + 6xy" =0; y (x) is bounded as x — ► 0, 

*/(l) = #'(l) = 0. 

267. x 2 y" + xy' — 1/ = 0; is bounded as x— <-0, 

1/(1) = 0. 

268. + ^t/ = 0; i/(0) is finite, i/(l) = 0. 

269. x 2 y" + xy' — n 2 y = 0; y(0) is finite, #(1) = 0. 

270. x 2 (ln.x— l)y" — xy' + y = 0; y(0) is finite, 

</(i) = o. 

271. 1 =0; t/(0) = 0, ^(1) is finite. 

272. xy" + y' = 0; y(0) is bounded, y(l) = 0. 

273. */"-</ = (); y(0) = y'(0), </(/) + Xt/' (/) = 0. 
(Consider the cases: %=\, X = — 1, |X|=^1.) 

21. Using Green's Function in the Solution 
of Boundary-Value Problems 

Let there be given a nonhomogeneous differential equation 
L [y]=p (x) t/<"> (x)+ Pl (x) y"-" (x)+...+p n (x)y(x)=f (x) (1) 
and the boundary conditions 

v 1 (y) = o, v 2 (y) = o v n (y)=o (2) 
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As in Sec. 20, we consider that the linear forms V lf V it . . . , 
V„ in y{a), y' (a), . . ., */<»-» (a), y(b), y' (b), . ..,y«-»(b) are 
linearly independent. 

Theorem. // G(x, £) is Green's function of the homoge- 
neous boundary-value problem 

V k (y) = 0, (k=l, 2 n). 

then the solution of the boundary-value problem (l)-(2) is 
given by the formula 

b 

y{x)=\G(x,l)f{\)dl (3) 

a 

(see [20]). 

Example 1. Using Green's function, solve the boundary- 
value problem 

y"(x)-y(x) = x, (1) 
</(0) = </(l) = (2) 

(a) Let us first find out whether Green's function exists 
for the corresponding homogeneous boundary-value problem 

y"(x)-y(x) = 0, (V) 
y(0) = y(\) = (2') 

It is obvious that y x (x) = e x , y 2 (x) = e~ x is the fundamental 
system of solutions of the equation (1'). Hence, the general 
solution of this equation is 

y (x) = Ae x + Be~* 

The boundary conditions (2) are satisfied if and only if 
A = B = 0, i. e., y(x) = 0. Thus, Green's function exists. 

(b) It can readily be verified that 

■ sinhs sinh (1 — 1) rj<JC< p 

G(x, £)H • u Si °u n ' ^ ° (3) 
° » sinh 1 sinh (x— 1) , . K ' 

sinhl ' k 5 ^* 5 * 5 
is Green's function for the boundary- value problem (l')-(2'). 
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(c) We write the solution of the boundary-value problem 
(l)-(2) in the form 

i 

y{x) = \G(x,\)Ul (4) 



where G (x, £) is defined by formula (3). 

Splitting up the interval of integration into two parts 
and substituting from (3) into (4) the expression for Green's 
function, we obtain 

X 1 

, , p \ sinh j sinh (x—\) M f j sinh x sinh (|— 1) ^ e 

y( x > = ) ien dl+ ) — ion d6 = 

X 

•^^-^jisinhidi + ^jisinhd-Ddi (5) 



sinh 1 

o 

But 



j £ sinh % d\ = x cosh x — sinh x, 



1 

Jlsinh (£ — l)d£=l— xcosh (x— l) + sinh(x — 1) 



and therefore 
1 

sinh 1 



y (x) = . . {sinh (x — 1 ) [x cosh x — sinh x) + 



+ sinh x [1 — xcosh (x— 1) + sinh (x— 1)]} = ^jjjjj * — * 

Here we take advantage of the formula 

sinh (a ± (J) = sinh a • cosh p ± cosh a sinh |5 

and also the oddness of the function sinh a:. 

Direct verification convinces us that the function 

, , sinh x 
^W=iinhT-* 

satisfies equation (1) and the boundary conditions (2). 
Example 2. Reduce to an integral equation the follow- 
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ing boundary- value problem for the nonlinear differential 
equation: 

y" = f(x,y(x)), (1) 

y(0) = y(l) = (2) 

Constructing Green's function for the problem 

y" = o, (3) 

y(0) = y(l) = (2) 

we find 



G( *- 6 > = \-(x-l)6. 



(g-i)*, o<x<i, 



Regarding the right side of equation (1) as the known func- 
tion, we get 

i 

y(x) = \G(x,l)f(t,ym)d$ (4) 



Thus the solution of the boundary- value problem (l)-(2) 
reduces to the solution of a nonlinear integral equation of 
the Hammerstein type (see Sec. 15), the kernel of which is 
Green's function for the problem (3)-(2). The significance 
of the Hammerstein-type integral equations lies precisely in 
the fact that the solution of many boundary-value problems 
for nonlinear differential equations reduces to the solution 
of integral equations of this type. 

Solve the following boundary-value problems using Green's 
function: 

274. y" + y = x; y(0) = y(j)=0. 

275. </ IV = l; 0(O) = */'(O) = /(1) = */"'(1) = O. 

276. xy" + y' = x; t/(l) = #(<?) = 0. 

277. y" + n'y = cos nx; y(0) = y(l), </' (0) = */' (1). 

278. /— */ = 2sinh 1; i/(0) = */(l) = 0. 

279. if- y = -2e*\ y(0) = y'(0), y(l) + y'(l) = 0. 

280. y" + y=:x*; y{0) = y(j)=0. 
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22. Boundary- Value Problems Containing a Parameter; 
Reducing Them to Integral Equations 

Many situations require the consideration of a boundary- 
value problem of the following type: 

L[y] = ky + h(x), (1) 
V k {y) = (ft=l,2 n) (2) 

where 

L [y } = Po (*) y w (x) + p t (x) tp-»{x)+...+p„ (x) y (x) 
V k (y) ^a k y (a) + c#> y' (a) + . . . + <~ u t/ (n-1, a + 

+$ k y(b) + W > y'(b)i-... +pr ,, y ,n " 1, (ft) (* = i.2, .... «) 

(the linear forms V t , V it . . . , V a are linearly independent); 
h (x) is a given continuous function of x; k is some nume- 
rical parameter. 

For h(x) = we have the homogeneous boundary-value 
problem 

l [y] = ty* I 

V„(y) = 2, ...,n) j < 3 > 

Those values of X for which the boundary-value problem (3) 
has nontrivial solutions y(x) are called eigenvalues of the 
boundary-value problem (3); the nontrivial solutions are 
called the associated eigenf unctions. 

Theorem. // the boundary-value problem 

L[y] = o, v 

V k (y) = (k=l, 2, ...,n) } ( 4 ) 

lias the Green's function G(x, £), //ien boundary-value 
problem (l)-(2) is equivalent to the Fredholm integral 
equation 

b 

y(x) = %\Gix,%)y{l)d\ + f{x) (5) 

a 

where 

* 

/(*) = -$G(*,6)ft(6)da (6) 
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In particular, the homogeneous boundary- value problem 
(3) is equivalent to the homogeneous integral equation 

b 

y(x) = x\G(x,l)y(l)dl (7) 

a 

Note. Since G(x, £) is a continuous kernel, the Fredholm 
theory is applicable to the integral equation. Therefore the 
homogeneous integral equation (7) can have at most a coun- 
table number of characteristic numbers X lt X 2 , ...,X n , ... 
which do not have a finite limit point. For all values of X 
different from . the characteristic values, the nonhomogeneous 
equation (5) has a solution for any continuous right side 
f(x). This solution is given by the formula 

b 

y(x) = k\R(x,l;X)fa)dl + f(x) (8) 



where R (x, £; X) is the resolvent kernel of the kernel G (x, £). 
Here, for any fixed values of x and | in [a, b] the function 
R (x, I; X) is a meromorphic function of X; only characte- 
ristic numbers of the homogeneous integral equation (7) may 
be the poles of this function. 

Example. Reduce the boundary-value problem 

y" + Xy = x, (1) 

y(0) =y(j)=o (2) 

to an integral equation. 

Solution. First find the Green's function G(x, I) for the 
corresponding homogeneous problem: 

y°(x) = 0, | 

,(0) = ,(f)=oj (3) 

Since the functions y l (x) = x and y t (x) = x — y are, respec- 
tively, linearly independent solutions of the equation y" (x) = 
that satisfy the conditions i/(6) = and t/(-jj=0, we seek 
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Green's function in the form 

Vl (*) 02 (I) ()<- r ^t. 



where 
Thus, 



1 1 



2 



G(x, E) = J >? ( „ (4) 



Further, taking advantage of Green's function (4) as the 
kernel of an integral equation, we get the following inte- 
gral equation for y(x): 





where 

n 
: 





Thus, the boundary-value problem (l)-(2) has been redu- 
ced to the integral equation 



•A 
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Reduce the following boundary-value problems to integral 
equations: 

281. y" = %y + x*; y (0) = y (-f ) = 0. 

282. y" = Xy + e*; ^(0) = ^(1) = 0. 

283. y" + ^y = Xy+cos^; y(~\) = y {\), 

y' (-i) =«/'(!)■ 

284. y" + Xy = 2x+l; y(0) = y'(\), r/'(0) = </(l). 

285. y" = Ky+l; y (0) = y' (0) = 0, iT(l) = jr (1) = 0. 

286. #"' + Ai/ = 2x; </(0) = #(l) =0, */'(0) = </'(l). 

287. y" + %y = e*; </ (0) = </' (0), y{\) = «/' (1). 



23. Singular Integral Equations 

We shall call the following integral equation 

<p{x) = f(x) + % \ K(x, t)<p(t)dt (1) 

a 

singular if the interval of integration (a, b) is infinite or 
the kernel K(x, t) is nonintegrable [for example, in the 
sense of L 2 (Q)]. 

In the case of singular integral equations, things may 
occur which do not have any analogy in the case of a finite 
interval (a, b) and a well-behaved kernel K(x, t) (continu- 
ous or lying in L 2 (Q)). 

thus, if the kernel K (x, t) is continuous in Q {a ^ x, t^b} 
and a and b are finite, then the spectrum of the integral 
equation, that is, the set of characteristic numbers, is disc- 
rete and to every characteristic number there corresponds 
at most a finite number of linearly independent eigenfunc- 
tions (the characteristic numbers have a finite multiplicity). 

In the case of singular integral equations the spectrum 
may be continuous, that is, the characteristic numbers may 
fill whole intervals, and there may be characteristic num- 
bers of infinite multiplicity. 

Some examples will illustrate the situation. 
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We consider the Lalesco-Picard equation 

+ OS 

<p(x) = X J e-i*-'i<p(t)dt (2) 

— co 

The kernel of this equation, K(x t t)~e~ , possesses 
an infinite norm. Indeed, 

+ 00 +00 +00 +00 + 00 

J J K*(x, t)dxdt= J 5 e- l]K ~ t] dxdt= \ dx 



— 00 — 00 



If the function <p(%) is twice differentiate, then the integ- 
ral equation (2), which can be written in the form 



(p (x) = X 



e~ x $ e t (p(t)dt + e x \e~\(t)dt 



is equivalent to the differential equation 

<p"(x) + (2X— l)q>(x) = (3) 
The general solution of equation (3) is of the form 

(p(x) = C ie rx + C 2 e- rx (4) 

(Cj, C 2 are arbitrary constants), where 

r = Vl — 2X (5) 

Here, for the integral in the right-hand member of (2) to 
exist, it is necessary that |Rer|<l, that is, that X- be 
greater than zero for real X. Hence, in the domain of real 
numbers the spectrum of equation (2) fills the infinite inter- 
val 0<A.<-t-oo. Every point of this interval is a cha- 
racteristic number of equation (2) of multiplicity 2. Howe- 
ver, the associated eigenfunctions do not belong to the 
class L 2 (— oo, + oo). 

For X > y , sin V2X— Ix, cos V2X — \x are, by (4), eigen- 
functions; for X=-^ we get q>(x) = C l + C i x. Thus, for 
there exist eigenfunctions bounded in ( — oo, +oo). 
However, if the real part Vl — -2X is positive and less 
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than unity, then formula (4), for any choice of the con- 
stants C lt C 2 (C\ + CI 0), yields a solution of the integral 
equation (2) unbounded on ( — oo, +oo). 

This example illustrates the essential role of the class of 
functions in which the solution of the integral equation is 
sought. 

Thus, if we seek the solution of equation (2) in the 
class of bounded functions, then, as we shall see, all the 

values of X > are characteristic. 

But if the solution of equation (2) is sought in the class 
of L 2 -functions ( — oo, + oo), then for any value of X equa- 
tion (2) has only the trivial solution (p(x)^0, i.e., not 
one of the values of X is characteristic for solutions in 

M— °°. + «>)• 

Let F(x) be a continuous function absolutely integrable 
on [0, -foo] and having a finite number of maxima and 
minima on any finite interval of the x-axis. - 

Let us construct the Fourier cosine transform of this 
function: 

— +co 

F 1 {X)=y -2- j F(x)cosXxdx 



Then 

— +m 

P(x)=y A £ F x (X) cosXxdX 



Adding these two formulas, we get 

— +CD 

F 1 (x) + F(x)^ y ^ I [F , (t) + F(t)] cos xtdt 

o 

that is, for any choice of the function F(x) satisfying the 
above-indicated conditions, the function <p (x) = (x) + F (x) 
is an eigenfunction of the integral equation 

+ 00 

(p(x) = X J q>(t) cos xtdt (6) 



corresponding to the characteristic number X= y — . 



and take 
Then 
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Since F(x) is an arbitrary function, it thus follows that 
for the indicated value of X equation (6) has an infinite 
number of linearly independent eigenf unctions. 

This peculiarity of equation (6) is due to the fact that (6) 
is singular [the interval of integration in (6) is infinite]. 

Example. Consider the integral equation 

CO 

(p(x) = X J (f (/) cos xt dt (7) 



F(x) = e- ax (a>0) 
M*)= V i]e-'* cos xtdt= j/f^ 



Further, 

<p{x) = F(x)+F 1 (x) = e-'+\/'±-^ (8) 
Substituting cp (x) into equation (7), we have 

T f n a*+x* 

= K \]e-<* cos xt dt + /I ] dt] (9) 

As has already been pointed out 

00 

\ e~ at cos xt dt = „ a „ 



The second integral on the right of (9) may be found by 
using Cauchy's theorem on residues: 

P COS£/.,__Jl 


From (9) we thus obtain 
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It is clear from this that if X = then the function 

will be a solution of the integral equation (7). Hence, 
X — is a characteristic number of (7), and the function 

is the corresponding eigenfunction; now, since a is any 
number greater than 0, the characteristic number >„= J/ — 

is associated with an infinity of linearly independent eigen- 
f unctions (8). 

In similar fashion, we can show that equation (7) has 

a characteristic number X= — y — associated with the 
eigenfunctions 

288. Show that the integral equation 

<p (x) = X J q> (t) sin xt At 
o 

has characteristic numbers X, = ± of infinite multip- 

licity, and find the associated eigenfunctions. 

289. Show that the integral equation with a Hankel kernel 

<f(x) = X $ J, (2 Vxi)(f(t)dt 



[where «/, (z) is a Bessel function of the first kind of order v] 
has characteristic numbers A, = ± 1 of infinite multiplicity, 
and find the associated eigenfunctions. 
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290. Show that for the integral equation 

CO 

«PW = xJ^=j^«p(0* 

X 

any number K for which one of the values " + |/ k has a 
positive real part is a characteristic number. 
291. Show that the Volterra integral equation 



X 

q,(*) = *j(i-JL)q)(0d/ 



has an infinity of characteristic numbers A, = | + ' T 1» where 
the point (£, t|) lies outside the parabola l + 'H^O. 

The solution of certain singular integral equations with 
the aid of Efros' theorem (generalized product theorem). 

Let 

«p(x)==<D(p), 
u(x, x) = U (p)e- % i^ 

where U (p) and g(p) are analytic functions. Then 

00 

0>(q, (p))U(p)=\<p(x)u(x, x)dx (1) 



This is the generalized product theorem (theorem of Efros). 
If u(x, xy=u(x — x), then q(p) = p and we obtain the 
ordinary product theorem: 

00 

<D (p) U (p) = J <p (t) « (x— t) d-t 



If t/( P )= ' , q{p) = Vp, then 

W (x,x) = - 7 L-e"^ (2) 

Therefore if it is known that O (p) = cp (x), then by the 
Efros theorem we find the original function for - ^fj^ : 
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®(V p) 



go 



CJ)(T)<? "dt 



Example. Solve the integral equation 



CD Z" 

jix ! 



V XX J 





(3) 



(4) 



Solution. Let cp (x) = <£> (p). Taking the Laplace transform 
of both sides of (4), we get, by formula (3), 

®(V~p) i 



whence 



<t>(P). 



V) 



-V, or <D(p) = -=l 



Hence, (p(i)=l is a solution of (4). 
Solve the following integral equations: 



292 



293 



294 



oo 

nx J 

n 

05 

jix J 



cp (t)dt — 2x— sinh x. 



It is known that 



/ 2 ./ n (2 ]f t) = -L. e ~* (n = 0, 1, 2, ...) 

where J n (z) is a Bessel function of the first kind of order n. 
In particular, 



Jo (2Vt)=±e 
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By virtue of the similarity theorem 

X 

J a {2 V7t) = je~" 

whence it is seen that for the Efros theorem we should then 
take 

¥<P)-7 
Example. Solve the integral equation 

00 

<p(x) = xe-* + k $./.„ (2 Vxt)y(t)dt (\X\=£l) (5) 
o 

Solution. Let <p (x) = O (p). Taking the Laplace transform 
of both sides of (5) and taking into account the Efros 
theorem, we find 

0(p) = ( -^4-X>(i) (6) 
Replacing p by -j, we get 

From (6) and (7) we find 
or 

Whence 

(ttt+T^) 
Solve the following integral equations (A.=^±l): 

00 

295. <f(x) = e* + X J ]/ j J t (2 V7t)y(t)dt. 

o 

x 

296. q>(x) = C0SJC + A,$ J (2 KH) q> (/) dt. 

o 
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297. <p (x) = cos x + X j y J 2 (2 Vxt) cp (/) d*. 

o 

« , 

298. cp(x) = sinx + xf |/ jj, (2K^)<P(0^. 

o 

Solving certain singular integral equations with the aid 
of the Mel tin transformation. 

Let a function /(/) be defined for positive t and let it 
satisfy the conditions 

1 00 

J I i"'" 1 dt< + oo, l\f(t)\t°*-i<U < + <* (1) 

1 

for a proper choice of the numbers a 1 and o s . The function 

00 

/ r (s) = S/(0^-'d/ (s = a + it, 0l <a<a t ) (2) 
o 

is the Mellin transform of the function f(t). The inversion 
formula of the Mellin transformation is 

O+i 00 

W> = 5Hi j ^(s)r*ds (<>0, a 1 <o<o i ) (3) 

a-i oo 

where the integral is taken along the straight line /: Re 
s = a parallel to the imaginary axis of the s plane and is 
understood to be the principal value. When the behaviour 
of the function f(t) as / — ► and t — * oo is known, say 
from physical reasoning, then the boundaries of the strip 
(a lt o 2 ) may be established from the conditions of the 
absolute convergence of the integral (2). But if the beha- 
viour of /(/) is only known at one end of the interval 
(0, + oo), say as /— >-0, then only o t is defined, the straight 
line of integration / in (3) must be chosen to the right of 
the straight line a = a x and to the left of the closest singu- 
larity of the function F (s). 

The Mellin transformation is closely associated with the 
transformations of Fourier and Laplace and many theorems 
which refer to the Mellin transformation can be obtained 
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from the corresponding theorems for the Fourier and Laplace 
transformations by means of a change of variables. 

The convolution theorem for the Mellin transformation is 
of the form 

M||/(0<p(f)-f}=/ : '(s)-<D(s) (4) 

From this we can conclude that the Mellin transformation 
is convenient in the solution of integral equations of the 
form 

00 

<p (*)=/(*)+ J K (jjvWj (5) 
o 

Indeed, let the functions <p(x), f(x) and K(x) admit the 
Mellin transformation, and let (p (x) — -* <P (s), f(x) — >F(s), 

K(x) — <• K (s); the domains of analyticity F (s) and K (s) 
have a common strip o 1 < Res = a < a 2 . Taking the Mellin 
transform of both sides of equation (5) and utilizing the 
convolution theorem (4), we get 

<t>(s) = F(s) + K(s)G>(s) (6) 

whence 

(P(s)= F( l } (K(s)=£l) (7) 

l— /C(s) w v 

This is the operator solution of the integral equation (5). 
Using the inverse formula (3), we find the solution <p(x) of 
this equation: 

a-« oo 

G>nsider the integral equation 

00 Jt 

q>(x) + ! j e_T <P(0T- («>0) (9) 

o 
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Taking the Mellin transform of both sides of (9), we get 

00 CO 

M {<?-«} - j e~' x X s - 1 dx = a" j e~ z z*' 1 dz = F (s), 



M{±.e-*}=±Y(s)^K(s) (Res>0) 

so that the domains of analyticity of F{s) and K(s) coin- 
cide. The operator equation corresponding to equation (9) 
will have the form 



0( S )=£M+i-r(s)0(s) (10) 

(D(s) : 



a* ~ 2 

whence 

r (s) 

'^[i-i-r(s)] 
By the inverse formula (8) we obtain 

a+tco 

We find the integral (11) with the aid of Cauchy's integral 
formula. 

For <xx> 1, we include in the contour of integration the 
semicircle lying in the right half-plane. In this case, the 
sole singularity of the integrand lies at the point s = 3 at 
which 

i-|r( S )=o 

Then 

2 

T v ' (a*) 3 1|> (3) * 

where i|)(3) is the logarithmic derivative of the r -function 
at the point s = 3: 

, r'(3) 3 

(y is Euler's constant). 
For ax<l, the singularities of the integrand are the 

6 -1583 
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negative roots of the function 1 — £- T (s) so that 



cp (x) = — 2 V, « \ ; ; . ax < 1 

*Ti (ox)'* * < s *> 

where \p (s ft ) are values of the logarithmic derivative T(s) 
at the points s = s k (k= 1, 2, ...). Thus, 



fW = ' 



4 

(3— 2y) (ow) 3 ' a * > !> 

CO 

— 2 S — < 1 



Let us consider an integral equation of the form 

<p(x) = f(x)+\K(xt)<f>(t)dt (1) 



(the Fox equation). Multiplying both sides of (1) by X s ' 1 
and integrating with respect to x between the limits 
and oo, we get 

oo 00 00 oo 

J y(x)x s ~ 1 dx = lf(x)x s - 1 dx+ i j l ip(t)di J K (xt) x 1 ' 1 dx 



Denoting the Mellin transform of the functions q>(x), f(x), 
K(x), by <J>(s), F(s), K(s), respectively, we get, after 
simple manipulations, 

00 

®(s) = F(s) + K(s) i j i cp(t)t- s dt (2) 



CO 

It is easy to see that \ q>(t) t~ s dt = Q) (I — s) so that (2) 



will be written in the form 

0(s) = F(s) + 0(l-s)K(s) (3) 

Replacing s by 1 — s in (3), we get 

<b(l-s) = F(l-s) + ®(s)K(\-s) (4) 



FREDHOLM INTEGRAL EQUATIONS 163 

From (3) and (4) we find 

<b(s) = F(s) + F(l-s)R(s) + <b(8)K(s)-R(l—s) 
whence 

f( S ) + f(l-s)/(( S ) 
\-K(s)-K(l-s) K ' 

This is the operator solution of equation (1). 
Using the inverse Mellin formula, we find 

i e r fw+zo-.)^ .^ 6) 

2ni J 1-/(( S )/[(1- S ) v ' 

c- Zoo 

which is a solution of the integral equation (1). 
Example. Solve the integral equation 

— * 

<p(x) = f(x) + Xy ^- f <p(t)cosxtdt (7) 

o 

Solution. We have 

/ — 00 

AT(s) = X y jx^cpsxdx (8) 

o 

To compute the integral (8) we take advantage of the 
fact that 

CO 

]e- x x z ~ 1 dx = T{z) (9) 



If in formula (9) we turn the ray of integration up to 
the imaginary axis, which by virtue of the Jordan lemma 
is possible for <z< 1, we arrive at the formula 



ITU 



\e-'' x x*~ 1 dx = e a r (z) 



Separating the real and imaginary parts, we get 

^x z ~ 1 cosxdx = cos^j-T (z), (10) 
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j sin X dx = sin j T (z) (11) 



Thus, by virtue of (8) and (10), 

~K(s) = X|/£r(s) cos -™ (12) 

Also, 

(s) • /C (1 -s) = ^ jA| r < s > cos T ' k r ( ! - s ) sin ? = 
= ^-2cosy siny r(s)r(l— s) = X 2 

since T (s)-T (\—s) = -£-. Hence, ii M If (*)} = F (s), then 

sin jxs 

by formula (5) (for |a.|=^ 1) 

0(s) = F(s) + F(l-^ (S ) 

and therefore 

CT + ioo 

+ F(1 — s)X Y 4" r ( s ) cos ^]x- s ds = 

-fieri I ^ 

a- ice 

O + ice 

+ T=Tpy 7H i rWcos2F(I-.)*-ds (13) 

a-i<» 

In the second integral on the right of (13) replace F(l— s) 

co a+(ou 

by §f(t)t~ s dt and note that gL j F(s)x- f ds = /(*). 

a-(co 

Then formula (13) can be rewritten as 
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a + (oo 



J_ 

O-iai 



+rhYir<k [ rwco8^.(x/)-dsJ/(o<tt (H) 

a-ioo 

By Mellin's inverse formula, 

J r(s)cos-^-(jc/) _,s ds = cosA;/ 
so that, finally, we have 

b 

Solve the integral equations: 

00 

1 1 f 

299. <f ) W = -rq-^ + y^"J <V {t) cos xtdt. 

o" 

CD 

300. q>(je) = /(*) + *, j/~ -| j q> (/) sin xt dt. 

o 

OD 

2 r 



301. 



cp(x) = — e x + -y- J (p (/) cos ^ d/. 



CHAPTER III 



APPROXIMATE METHODS 



24. Approximate Methods of Solving Integral Equations 

1. Replacing the kernel by a degenerate kernel. Suppose 
we have an integral equation 

b 

<t{x) = f(x) + %\K(x, t)<p(t)dt (1) 

a 

with arbitrary kernel K(x, t). The simplicity of finding a 
solution to an equation with a degenerate kernel (see Sec. 15) 
naturally leads one to think of replacing the given arbitrary 
kernel K(x, t) approximately by a degenerate kernel L(x, t) 
and taking the solution cj>(x) of the new equation 

b 

(f>(x) = f(x) + \ ^ L(x, t)(f>(t)dt (2) 

a 

as an approximation to the solution of the original equa- 
tion (1). For the degenerate kernel L(x, t) close to the given 
kernel K {x, t), we can take a partial sum of Tailor's series 
for the function K (x, t), a partial sum of the Fourier series 
for K (x, t) with respect to any complete system of functions 
{u n (x)} orthonormal in L 2 (a, b), etc. We shall indicate some 
error estimates in the solution (1) that occur when replacing 
a given kernel by a degenerate kernel. 

Let there be given two kernels L(x, t) and K(x, t) and 
let it be known that 

b 

\\K(x, t)-L(x, t)\dt<h 

a 

and that the resolvent kernel R L (x, t\ X) of the equation 
with kernel L(x, t) satisfies the inequality 

6 

l\R L (x, t; k)\dt<R 
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and also that \f(x) — <t]. Then, if the condition 

1-|M*(1+|M*)>0 
is fulfilled, it follows that the equation 

6 

<P (*) = X^K( X> t) q) (t) dt + f (x) 

a 

has a unique solution <p(x) and the difference between this 
solution and the solution $ (x) of the equation 

b 

<J) ( X ) =Mx) + X$L (*.<) <P (0 # 



does not exceed 



l<PW 9WI ^ l-|X|Ml + |X|/?) + T1 



(3) 



where JV is the upper bound of |f(x)| (see [8]). 

For the degenerate kernel L(x, t), the resolvent kernel 
R L (x, t\ %) is found simply (to within the evaluation of the 

n 

integrals); namely, if L (x, t) = 2 ( x ) (0> then, put- 
ting 



we get 
where 



\x k {x) T s (x)dx = a sh 



R L (x, t\ X) — D ( *' ( ^ X) 



(4) 





o x,(o .. 


X n (x) 








T t (t) l-la u .. 




» 


(5) 




r„(/) -Ki-- 


i -k» 




* 




1 — A.a n — A.a 12 . 


. . — Xa 


In 




D(X) = 


— A.a 2l 1— Xa 2SS 


. . — X-a 


2n 


(6) 




— Xa nl —Xa ni 


••• 1-H, n 
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The roots D (X) are the characteristic numbers of the kernel 
L(x, t). 

One more estimate (A,=>1). Let 

K(x, t) = L(x, t) + A(x, t) (7) 

where L(x, t) is a degenerate kernel and A(x, t) has a 
small norm in some metric. Also let R K (x, t), R L (x, t) be 
resolvent kernels of the kernels K(x, t), and L (x, t), res- 
pectively, and let ||A||, \\Rl\\ t> e the norms of the 
operators with corresponding kernels. Then 

||«p-9||<||A||.(l + ||/? A .||).(l + \\R L \\)- H/ll (8) 

The norm in formula (8) can be taken in any function space. 
The following estimate holds true for the norm of the 
resolvent kernel R of any kernel K(x, t): 

ii^ii< i-KiUii, (9) 

And in the space C(0, 1) of continuous functions, on the 
interval [0, 1], 

i 

max [\K(x, t)\dt 

||/||= max |/(x)| (10) 
o<*< 1 

In the space of quadratically summable functions over 

i 

b b \ — 



= (j/*(*)d*) 



(11) 



Example. Solve the equation 

i 

<p (jc) = sinx-|- } (1— x cos xt) q> (t) dt (1) 
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by replacing its kernel with a degenerate kernel. 

Solution. Expanding the kernel K(x, t) = l — xcosxt in 
a series, we get 

Let us take the first three terms of the expansion (2) for 
the degenerate kernel L(x, t), 

L(x, f)=l-x + ^ (3) 

and solve the new equation 

i 

9(;c)=sinje + j(l-jc + ^p-)$(f)d/ (4) 
o 

From (4) we have 

cp (x) = sin x + C t (1 — x)-\-C 2 x i (5) 

where 

i i 
C, = C 2 = -ij>?(/)d/ (6) 



Substituting (5) into (6), we get a system for determining 
C, and C 2 . 
We have 
i 

C,= jlsinZ+C^l-O + C,/ 3 ]* =-^Ci+tC, + 1 -cos 1, 
o 

i 

C a = y J [<» sin / + C, (<* - f ') + C 2 / 5 ] d/ = 



= ^C 1 + ^C 2 + sin 1 — 1 +yCOsl 

or 

. yCj— ^-C 2 = 1— COS 1, 

— ^ C i + T5 C 2 = sin 1 + yCOS 1 — 1 

6*- 1583 



(7) 
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Solving this system, we get 

^ = 1.0031, C 2 = 0.1674 

and then 

<p(*)= 1.0031 (1— x) + 0.1674jc 3 + sin x. 

The exact solution of the equation <p(x)=l. 
Now let us estimate ||<p — <p|| using the formula 

||cp-H<IIA||-(l+||flJ)(l + |l«JHini (8) 

In the metric of the L 2 -space we get 

i 

( l 1 IT 



\\A\\<^)]x»t°dxdt> = WTf < 238. 



{ 1 1 ) T 
\\ K \\^ I ^ (i — x cos xtydxdti = 

loo j 
= |2cosl— icos2 + Tgsin2— 2 <|-, 



i i 





1 



= |jsin^r= ^ 2 - si " 2 



We estimate the norms of the resolvent kernels R K and 
/? L using the formulas 

|| n || < _J1A!L_ II D II <- HM 

II^H 5 ^ • ll^ill^a 1_|X|.|| L || 

where |X| = 1. Hence, || 4- . ||/? t ||<4- but then 



2 ' 11 ' x tii^= 2 

238 1 



H-^ll<i(l + T)(l+T)-T<0.016 



Find the solution of the integral equation in the following 
cases by means of substituting a degenerate kernel for the 
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kernel, and estimate the error: 
1 

302. (p(x)=e x —x— < \x(e xt —\)(p(t)dt. 

o 

i 

303. <p(x) = x + cosx + } *(sin*f — l)(p(t)dt. 



1 

304. (p(x) = j(e-* + 3x— 1)+ ^(e~ xtl —l)x(f(t)dt. 



1 

305. q> (je) = — -|- — sin jc+ \ (1 — cos xt 2 )xq>(t)dt. 

o 

2. The method of successive approximations. The method 
of successive approximations (iteration method) consists in 
the following. 

We have an integral equation 

b 

<p(x) = f(x) + x\K(x, t)y(t)dt (1) 

a 

We construct a sequence of functions {<p„(*)} with the aid 
of the recursion formula 

b 

<p„ (x) = / (x) + X S K (x, t) tp,,., (t) dt (2) 

a 

The functions cp„ (x) (n = 1 , 2, . . .) are considered as appro- 
ximations to the desired solution of the equation; the zero 
approximation q> (jc) may be chosen arbitrarily. 
Under certain conditions 



f b b 

\X\<±, B = ]/ ^K*(x, t)dxdt 



(3) 



the sequence (2) converges to the solution of equation (1). 
The magnitude of the error of the (m+l)th approximation 
is given by the inequality 

I <P (*)- <tW W I < FC t B-* -^q 1 ^ + OQfl- 1 1 XB r 1 (4) 
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where 



F = ]/"$ / 2 (x) <&, O = ]/" J q>; (x) dx, 

a n 

C x =y max \K z (x, t)dt 



Solve the following equations using the method of succes- 
sive approximations: 



306. <p (x) = 1 + \ xt 2 <p (t) At. 

o 

i 

307. y(x)=-jx + -j^xt<p(t)dt. 



308. Find the third approximation q> 3 (x) to the solution 
of the integral equation 

i 

<p(x)=l + lK{x, t)y{t)dt 
o 

where 

It, x^zt, 
t ) = \ x, x^t 

and estimate the error. 

Observe that the basic difficulty in applying the method 
of successive approximations consists in computing the inte- 
grals in formulas (2). As a rule, it is performed with the 
aid of the formulas of approximate integration. Therefore, 
here too it is advisable to replace the given kernel by a 
degenerate kernel with the aid of a Taylor expansion and 
only then to introduce the iteration method. 

3. The Bubnov-Galerkin method. An approximate solution 
of the integral equation 

q> (*) = /(*) + t)y{t)dt (1) 
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by means of the Bubnov-Galerkin method is sought in the 
following manner. We choose a system of functions \u„(x)} 
complete in L 2 (a, b) and such that for any n the functions 
u, (x), u 2 (x), u n (x) are linearly independent and we 
seek the approximate solution q> n (x) in the form 

n 

Wn (x) = 2 a k u k (x) (2) 

The coefficients a k (k=l, 2, .... n) are found from the 
following linear system: 

(cp„(x), u k (x)) = (f(x)u k (x)) + x(^K(x, t)<p n (t)dt, u k (x) 
(fe=l, 2, ... , n) (3) 

b 

where (f, g) stands for J / (*) £ W and in place of q> n {x) 

a 

n 

we have to substitute 2 If the value of X in (1) 

is not characteristic, then the system (3) is uniquely solvable 
for sufficiently large n and as n — ► oo the approximate so- 
lution cp„ (x) (2) tends, in the metric L 2 (a, 6), to an exact 
solution <p(x) of equation (1). 

Example. Use the Bubnov-Galerkin method to solve the 
equation 

i 

<p(x) = x+ \xty (t)dt (4) 
- i 

Solution. For the complete system of functions on [ — 1, 1] 
we choose the system of Legendre polynomials P„(x) 
(n = 0, 1, 2, ...). We shall seek the approximate solution 
qp„ (x) of equation (4) in the form 

cp 3 (x) = a x -l+a 2 x+a 3 — — 

Substituting cp s (x) in place of q> (x) into equation (4), we get 

l 
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or 

! i 3a;*— 1 . 2 , K . 
fl^+fl^ + a, — ^- = x+x-^-a i (5) 

Multiplying both sides of (5) successively by 1, x, — ^ — 

and integrating with respect to x between the limits — 1 
and 1, we obtain 

2a x = 0, 
_2 __2_ _4_ 
3 3 9 

2 n 
-5^ = 

Whence aj = 0, a 2 = 3, a 3 = 0, and so cp 3 (a;) = 3*. It is easy 
to verify that this is the exact solution of equation (4). 

Use the Bubnov-Galerkin method to solve the following 
integral equations: 

i 

309. <p (x) = 1 + J (xt + x 2 ) <? (t) dt. 

-l 

i 

310. <f(x) = l+-^x+ j(*< 2 — x)q>(t)dt. 

-i 

311. (p(x)= 1— x(e x — e~ x )+ $ x*e xt <p (/) dt. 

-i 

Note. The Bubnov-Galerkin method yields an exact solu- 
tion for degenerate kernels; for the general case it is equi- 
valent to replacing the kernel K(x, t) by the degenerate 
kernel L(x, t). 

25. Approximate Methods for Finding Characteristic Numbers 
I. Ritz method. Suppose we have an integral equation 

b 

<p(x) = X $ K(x, t)<p(t)dt 

a 

with symmetric kernel K(x, /)==/((/, x). 

Choose a sequence of functions {ty n (x)}, (x) g L z (a, b) 
such that the system {i|)„(x:)} is complete in L 4 (a, b) and for 
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any n the functions ^ (x), ^(x), t|) n (*) are linearly 

independent on [a, b]. We assume 



q> n (x) = 2 a k y k (x) 



(1) 



and we subject the coefficients a k to the condition || q>„ || = 1; 
under this condition we seek the stationary values of the 
quadratic form 

CC<P«> <P„) 

We arrive at a homogeneous linear system in the coefficients 
a k (a is a Lagrange multiplier): 



2 {(K% %)-a (ty, %)} a h = 
(/ = !, 2 n) 



(2) 



For a nonzero solution (2) to exist the determinant of the 
system (2) must be equal to zero: 



CT M>i> ¥•)— °M>i. 

(/0|>„ o(t|3 2> (/0|>„ <T(t„ *,).. 

• • -(^2. *„)— *J 



= (3) 



The roots of equation (3) yield approximate values of the 
eigenvalues of the kernel K(x, t). The largest of the roots 
of equation (1) yields an underestimate of the largest eigen- 
value. Finding a from (3) and substituting into (2), we seek 
the nonzero solution a k (k = 1, 2, .... n) of the system (2). 
Substituting the values of a k thus found into (1) we get an 
approximate expression for the eigenfunction corresponding 
to the eigenvalue that was found. 

Example. Using the Ritz method, find the approximate 
value of the smallest characteristic number of the kernel 



K(x, t) = xt; a = 0, 6=1. 
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Solution. For the coordinate system of the functions ty n (x), 
choose the system of Legendre polynomials ip„ (x) = 
= P„(2x — 1). We confine ourselves to two terms in for- 
mula (1) so that 

<p, (x) = a 1 P (2x-l) + a 2 P, (2x — 1 ). 

Noting that 

ij, 1 = P (2a;— 1)= 1, ^ = Pj (2x— l) = 2x— 1 
we find 

i i 
OPi. TM=ldx=l; t,)=S (2x-l)dx = 0; 



1 

(%> = f(2*-l) 2 dx = y 

o 

Further 

i i ii 

*i) = j( (*)d* = J j , .^^i*ef/ = 4■ 
o \o / 
I 1 

(Kip,, = j jx/(2*— l)dxdt=± 

o 

I 1 

= j jx/(2< — 1) (2a:— l)dxdt = 



'36 



o o 



The system (3) then becomes 



T~° 12 



±_J_ 
36 3 ° 



= 



or 



(h + t)=° 



Whence ^ = 0, a 2 =-i- . The largest eigenvalue cr 2 =y , hence, 
the smallest characteristic number X = — = 3. 
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Use the Ritz method to find the smallest characteristic 
numbers of the kernels (a = 0, 6=1): 

312. K(x, t) = xH*. 

313. K(x, t)=l U * J'; 

{ \x(2-t), 

314. K(x, t) = \ \ 

\±t(2—x), jc>/. 

2. The method of traces. Let us use the term mth trace 
of the kernel K(x, t) for the number 

b 

4- t)dt 

a 

where /C m (^, /) stands for the mth iterated kernel. 

The following approximate formula holds true for the 
smallest characteristic number K lt given sufficiently large m: 

Formula (1) yields the value \ in excess. 

Traces of even order for the symmetric kernel are com- 
puted from the formula 

b b b x 

A 2m =\\ K* m (*, t) dxdt = 2 J $ K* m (x, t) dt dx (2) 

a a a a 

Example. Using the method of traces, find the first cha- 
racteristic number of the kernel 

«<*•<>-{*: *</■ "- 1 

Solution. Since the kernel /((*, /) is symmetric, it is 
sufficient to find K t {x, t) only for t <*. 
We have 
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I t 

K t (x, t) = §K(x, z)K(z, t)dz = §z*dz + 



+ \ j ztdz+\ j xtdz = xt- 

t X 





xH t 3 



Then, from formula (2), we find for m = 1 and m = 2, 
respectively, 

1 X I X 1 

A t = 2^dx^K\{x, t)dt = 2 f dx ^t 2 dt = 2^~dx = 1 , 



1 X 

,4 4 =2 Jdx $/C5(*, 0* = 

' 





1 

3 + 12 + 7-36 3 15 + 30 J \i=o' 

o 
l 



V 3 +12 + 7-36 3 15 + 30J flX_ 630 

o 

Then by formula (1) 

K ™\ ~IL =2AS 

' 630 

Use the method of traces to find the first characteristic 
number of the following kernels (a = 0, b=l): 

315. K(x, t) = xt. 

316. K(x, t) = xH\ 



317. K(x, t) = 



~x(2—t), 
jt(2-x), x^t. 
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318. K( 



\-Vxt 



Vxt\nt, x</, 
lnx, x^t. 



3. Kellogg's method. Let K(x, t) be a symmetric kernel 
which for definiteness we will consider to be positive defi- 
nite, and let to (x) be an arbitrary function in L 2 (a, b). 
We construct a sequence of functions 

b 

coj (x) = J K (x, t)a>(t)dt, 



to 2 (x) = lK(x, t)a x {t)dt, 

a 
b 

co„(%) = S/((x ( 0«»»-i(0* 



(1) 



and consider the numerical sequence 



< II <a»-i II \ 
I llcojl / 



(2) 



Let cp x (x), cp 2 (x), ... be orthonormal eigenfunctions of the 
kernel K (x, t) and let ^ ^ X 2 ^ . . . be the associated 
characteristic numbers. Further, let the function co(jc) be 
orthogonal to the functions cp x (x), cp 2 (x), . . . , cp ft _! (x) but 
not orthogonal to the eigenfunction cp A (x). Then the sequence 
(2) has as its limit the kth characteristic number % k . 

The sequence of functions | |j p | converges in this 

case to some function, which is a linear combination of 
eigenfunctions that is associated with the characteristic 
number X k . The sequence 

~ (3) 



{Kll«>-ll } 



converges to the same limit as does the sequence (2). If 
(co, <pj) 0, we get two approximate formulas for the smallest 
characteristic number: 



^~ IN. II 



(4) 



180 



INTEGRAL EQUATIONS 



X, « „ 7 (5) 

Formula (4) yields the value of X t in excess. If the kernel 
K{x, t) is not positive definite, then formulas (4) and (5) 
yield an approximate value of the smallest absolute magni- 
tude of the characteristic number of the given kernel. For 
an apt choice of u>(x) the Kellogg method is comparatively 
simple for computation. 

The drawback of the method is that we do not know 
beforehand which of the characteristic numbers has been 
determined. 

Example. Use the Kellogg method to compute the smallest 
characteristic number of the kernel K(x, t) — xH i , O^a:, 
/< 1. 

Solution. We take <x>(x) = x. Then 



to j (x) = j* 2 / 2 /d/=^, 



1 



1 



Further, 

IKWIIt- \ x>dx - 



Thus, by (4), 

1 l I 



4.5"-!V 5 



11 -= 5 

4 * 5"-» ' yi 

Using the Kellogg method, find the smallest characteristic 
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numbers of the following kernels: 

319. K(x, t) = xt; 0<jc, /< 1. 

320. K(x, t) = s'mxs'mt; — n^.x, t^.n. 



321. K(x, /) = rj X x ^ t : 0<jc, 
1*(2— 0, 

322. K(x, /)= J , 0<Cx, 1. 
4<(2— 
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9. <p(x) = — x+\(t—x)(p(t)dt 



X 

10. q>(jt)=l + $<p(/)d< 



x 

11. <p (x) = cos jc — ^(x — t)(p(t)dt 



X 

12. (}>(*) = 5— 6x + \ [5— 6(x— /)] q>(/)# 

o 

13. <p (x) =cos* — x — ^(x — t)(p(t)dt 

o 

x 

14. <p(x) = x—s'mx + e x (x—\) + j [s'mx— e*(x— t)]<p(t)dt 

X 

15. (p(x) = cosx— 2x(l +x 2 )— $(1 + x 2 )(;c— t)q>(t)dt 

o 

16. cp (x) = xe* + 1— a; (a; 2 — 1) — 

X 

-j [je+-5-(jc»-Jt)(x-/)»] «p(0* 
o 

17. (p(a;) = a:(a;+1) 2 + 5a:(x:— tf<f{t)dt 

o 

19. pi=- sinh VT(x— t)(l>0) 20. e <i+wu-» 

21. eX<*-» e * 2 -t* 22. l±^ e x«*-« 

23 ?±£2i£ e xu-« 2 4 cosh * c * <*-» 

2+cos< ' cosh* 

25. a*"'^ 1 *-" 
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26. <?*-'(*— 1 + 2) 27. L e x-t_^ e - 3 u-n 

28 2xe xi - tl 29 _i^i±_L LJ?__ A e - a (*-»] . 

2 (2/+ 1)2 [4< 2 +l J" 

one of the solutions of the corresponding differential equa- 
tion y 1 (x) = e~ 2x . 

31. ^L- sinh V~2(x—t) 32. 1 33. (*— 0«- u -° 

34. «T~ Jcosh (x—t) + y=- sinh t) 

35. 2e JC - t (l+x— 36. <p(x) = e 2x 

1 1 2 

37. <p(x) =-§e zx — -g- cos sin* 

38. <j>(x) = 3*(l— e~ x ) 

39. q>(je) = e*sinx + (2 + cos*)e*ln 5 ^ n 

40. cp(jt) = e* 2 -*— 2a: 41. <p (a:) = e* 2+2 * (1 +2*) 

42. cp(x) = e*(l -H^ 2 ) 

43. <P(*)= +* arctanA:— -j In (1 -fx 2 ) 

44. <p(x) = e~(*+l) — 1 45. <p(jc)=e-* + l) 

46. q>(x) = s'mx 47. (p(jc) = cosx 

48. (p (x) = cosh x 49. cp (a:) = 1 50. cp (a:) = a; 
51. <p (*) = <?* 52. «p(a:) = 2 53. <p(a:) = 2 
54. <p(x) = x 2 —2x 56. <p(*) = 

57. (p 2 (A-)=l+* + |-* 2 + 4* 3 +2l^ + T x5 + TS^ + ^ ? 



x 4 , x 10 



58. <p,(*) = _ * + T ~ 14 + T60 59, , PW= 1 
60. cp(x) = A:— y 61. (p(x) = y(3e 2JC — 1) 

62. (p(A:) = sinx 63. cp (x) = y (2 cos V3 x + 1 ) 

x ? 

64. (p (a:) = 2a: + 1 65. (f(x) = x + -^ 

1 1 JC* 

66. (p (a;) = y sin x + y sinh a: 67. (p(aj) = a: — -g- 
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68. y(x) = e* 69. q>(x) = ^=sinh-^-A; e a 

70. cp (a:) = 1 + 2xe x 71. cp (x) = e* ( 1 + x)* 

72. cp W = !l±££!il±iilf 

73. cp, (x) = sin x, cp 2 (x) = 

74. ((> 1 (x) = 3e x — 2, <p 2 (x) = 3e x — 2e* x 

1 — e 2x 

75. <p x (x) = e 2x , <p 2 W = — — 



76. 



q>i (x) = (* + 2) sin x -f (2x + 1 ) cos x, 
(*) = — g- cos x y~ sin •* 



77. cp t (x) = 2 sin x, cp 2 (x) = 2 cos x— 1 , cp 3 (x) = x 

78. q>! (a;) = cos x, cp 2 (x) = sin x, <p 3 (a:) = sin x + cos # 

cp t (x) = ( 1 + cos x + y sin x, 

79. ■! cp 2 (x)==l— x + y sinx— (l+|-)cosx ( 
cp 3 (x) = cos x — 1 — ^l-f-y^sinx 

80. (f(x) = e x —l 81. q)(A:) = — e* 

82. cp (x) = sin x 83. cp(x)= 1— e-* — xe~ x 

84. tp (jc) = 1 — cos* 85. cp (.*:) = 1 — x + 

+ 2 (sinx— cos x) 
86. cp (*) = (]— x)e~* 87. cp (x) = c + 2e~ x 

a e (a-l)x 

88. <p(x) = cosx— sinx 89. cp (x) — ^—^ ^— j — 

90. cp(x) = cosx — sinx 91. cp(x)=l— xln3 
92. q> (*)=/'(*)— f(x)\na 93. q>(x) = xe x ' 

94. (f(x) = xe i 95. cp(x) = e 2 (x 2 + 2) — 1 
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107. <p ( x)=±f-i£iL<tt 

n J V x—t 



* 1 

108. (pW =i^ +e -j , ^r»d/) 

109. <p(*) = | 

110. q) (,,y) =ag l i (p_g), where 

(• f(u,v)dudv 



and Dj is a right isosceles triangle with vertex at the point 
(x, y) and hypotenuse on the Ou-axis of the UOV-p\ane. 



111. cp(x) = - 



113. (p(x): 



2x'> 



112. (pW = -^= 

y x 



(I 



2* 4 



r -T- * 



114. cp(x) = 3 

115. <p(jc) = sin;c 

118. cp(x) = ^A; 

120. (p(x) = 2x — x 2 

122. <f>(x) = 3 !(*<>"*— x 2 e~*) 123. y(x) = J (x) 

124. <p(*) r l— £ 125. <p(jc)=l+2x + f + f 

126. We have x 2 — * a = x 2 — 2x/ + * 2 + 2x/ — 2/ a = 
= (x— /) 2 + 2*(x— /). Therefore 



116. q>(x)= 1 117. <p(x) = e -* 
119. q>(x) = cosx — 2sinx 
121. <p(x) = 2sinx 



Taking transforms and applying the product theorem and 



7-1583 
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the theorem of differentiation of a transform, by virtue of 
which t<p (t)=— O' (p), we get 

or 

®' (/>) = -]> 

Solving this differential equation, we find 

0(p) = Cp + l 

Since O(p) is an image function, <!>(/?) should approach 
zero as p oo, so that C = and, hence, (p) = ^- , whence 

q>(*)= T . 

127. <p(x) = C— x 128. <p(x) = C + /,(2j/jt) 

129. q>(x) = C-f-x 130. q)'(Ar) = 2 + 6 (at) — 6' (*) 

131. q> (x) = 8 (jc) — sin* 

132. «p (x) = 8 (x) + 3 133. (p (x) = 1 +x + 8 (x) + 8' (x) 
134. <p(*)= 1 135. <p(x) = J 1 (x) 

136. <p(x) = — h( x )> h( x ) i s a modified Bessel function 
of the first kind. In Problems 141 and 142 the indicated 
functions are not solutions of the corresponding integral 
equations, but in Problems 137 to 140 and 143 to 145, the 
indicated functions are solutions. 

2x— t+(x+t— 2xt — ^] X 
146. R (x, t; X) = ^ — xw 



147. R (x, t; %) = V * b > 

^240 

148. R (x, t; X) = sin x cos t 

D , , ^ . _ sin x— sin t — n (1 +2 sin x sin t)X 

149. K (X, t; K) — l+2n*Jl* 
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x+t + l+2(xt + ±)x 

150. R (x, t; X) = ±j 6 -L- 

1— 2X--^X* 



l+3xt+(3^—3xt—\) 
151. R(x,t;X) = * — — '- 



1— 2X+± X 2 

Axt—x*-(2xH — ^rX i + x-^r-xt) 
152. R(x,t;k) = ^ ^-js — i- 



e*- 1 



153. R(x, t; X) = yzzi 

154. R(x,t; a) = »»"(* + 0+^ »(*-') 

,„ r> / y . 1A x— sin h 2 (g-H-x sinh fl A, 
1J>{). f, A) = [_|_4 e -ix,2 

156. (p(je)= 1 

, . I [ , ($x—2)X—X*x-\ 

157. q>(*) = -g-|* + * X »-k+6 J 

158. <p (x) = cos 2x 

159. cp(x) = yC JC 

.en / \ 3x(2X— 3X* + 6) 

161. /C^.JX, 0=(-y) n_1 (^-0. 

K 2 „(x, = 2(-l)-(y)"" 1 (x< + y)(n=l, 2, 3, ...) 

162. f) = 8ln( 2 + — ^cos(x-0, *.(*. /) = 

= ~16~ Sln (JC_ ^ 

163. /C,(x, = |-(^ + 3 + 2^ 2 + 4^+4 

/c 3 (x, o=g(*«+/«)+-§-*v»-f **+^ 

164. /C 2n -i(x, *) = (2n)«»- , (x + sin0 

/C,.(x, t) = (2nr>-i(l+xsmt)(n = l, 2, ...) 

7* 
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165. K n (x, t) = xe f 

'(>*+'+ 



166. K n (x, /) = (-l)"- 1 (^tiy-^cost 



167. K 2 (x, = \ e x+t i e i-x-t 

1 - — ±| \-(x—t—l)e"- t ,t^x^l 

e2 + l e t ~ x , — 1 <x<0 



168. t) = \ e *\ x 



e i+x , 0<x< 1 



2 

169. /; ^) = 2 _^_i )X : IM<^±T 

170. /?(*; /; X) = 2si 2 n ^ os ' ; |M<2 

171. R(x, t; X) = ^; |M<{ 

172. ^) = 3(1 t^-° ; IM<| 

173. /?(*, /; ^) = ^J; N<| 

174. R(x, t; ^) = 3^x; IM<i 

175. R(x, t; k) = sin xcos / + cos 2x sin 2t 

,76. *(,, /; X) = r l r + 3(2g -3 1 lf- 1) ; |X|<1 

180. <p(A;) = ^ r sin 2 A: + 2A;— ji 

181. (p(jc) = tanx 

182. <p(x) = -j- -X + cotjc 

183. cp W = r i±^ r 

184 . „ ( , H __^ + __L_,^1 

185- «PW= i-xr(p + l) 

2^ + C^. + ^ > ) In x 
186. <pW = ^9— ^ +i(l_4x) 



7 ..„. _ ,., . 15 .. , 5 
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187. cp(A:) = j-|^- sin x; 

188. cp (x) = Xn 3 sin x + x 

lon , . „ 2 cos x + jiX sin * 

189. = 2 4+ ^ 2 

190. q> (x) = Xn s'm x -\- cos x 

191. «p(*) = l§(x+l)* + ^ 

192. (pjW^ 0; (f> 2 , 3 M = ± 

193. (PAX)^0; (p 2 (*) = yX 2 ; (p 3 ,< (*) = ± jy=r x + ^ x' 

194. (p, (x) = 0; q> 2i 3 (x) = ± 3x 2 195. <p (x) = 
196. Has no solutions. 

198. X t = n _2 < <Pj (x) = sin 2 x 

199. No. 200. cp^x^sinx 

201. X,= — -|, ^ 2 = ^-; (p 1 W = sinx, 
cp 2 (x) = cos X 

202. There are no real characteristic numbers and eigen- 
functions. 

203. Xj = X 2 = — 3; (p(x) = x — 2x 2 

204 K = y~' <Pi(*) = |* + y* a 

1 3 

205. (pj (a;) = y a; + x 2 

206. X x = — -|- *, <}>! (a:) = sinh * 207. None 

208. There are no real characteristic numbers and eigen- 
functions. 

209. a.„ = — n 2 n 2 ; <p„ (x) = s in nnx (n=l, 2, ...) 

210. \ =1; <p (x) = e*; X n = — n 2 n 2 ; 

cp„ (x) = sin nnx -J- nn cos nnx (n = l, 2, ...) 

211. = — i^; q)„(x) = sinM, n x + n„ cos(i„x, where 
fi„ is a root of thvi equation \i — - = 2 cot pi. 
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212. X„ = 4n 2 — 1; <p„(x) = sin2nx (n=\, 2, .-..) 

213. ?,„=(n + y) 2 — 1; q>„ (x) = sin (n + y) * 



214. X„ 



sinl ' <P»(*) = sln M» + *) ( B=1 - 2 ' •••)' 
where (i„ are roots of the equation tan 2n\i — — fitanl. 
215. X„=l — (p„(x) = sinn„x + n„cos^ n x, where \i n 

are roots of the equation 2cotnn = n — -. 

4^-"; cp n (x) = sinn n x + n„cosu.„x, where u.„ 



216. K 



are roots of the equation 2 cot (j, = . 

217. k„ = — 1 — \x\; <p„ (x) = sin fi.„ x, where n„ are roots 
of the equation tan \a = |x (n > 0). 

221. (a)g; (b) £-J- : (c) i+fl* 

222. <p 1 (x)=l; <p a (x) = 2x— 1; 

223. cp^x^x; <jp.,(x)=x 2 

224. <pj (x) = 1 ; tp 2 (x) = x; cp 3 (x) = 3x 2 — 1 

225. ^0 = -^-. ?»W = 1; ^ = 1". <J>{ l, (x) = cos 2x, 
<p<*> (x) = sin 2x 

227. ^ = 2^, cp (x) = l; X„ = (-l)"^ 

tpn 1 ' (x) = cos <Pn 2) (x) = sin nx (n = 1 , 2, . . . ) 

228. (a)y, <p(x) = ±K3x; (b) | ; cp (x) = ± |/f x; 

(c) 1; q>(x) = ± ]/ -^-p 229. K = 3 

230. There are no bifurcation points. 



231. cp(x) = 



C-sinx, X = 

C • cos x, K = — 



o, 
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232. <p(x)-- 

234. <f(x) -. 

236. <p(je) = 

238. y(x) = 



233. <p(x) = C 



j C • arccos x, X—l 
"\0, 

■■C\x\ 235. <p(x) = C{x— x*) 

--sm^- 237. <p(x) = x— 2 + 2e* 



sin n*-)-sin |x (* — 1) — jtcosfijc . - n 
j r- , A. U, 

^ cos T( cos T + -2- sin !-) 

sinh (ix + sinh (x (x-\- 1) — |x cosh [ix ^ q 
2p, cosh ^cosh — -^-sinh-^-j 

where n = j/2X 



nan / \ r> i ^ V « sin 2/lX 

239. <p (*) = cos 2x + 4 2* (4n »_i )(4n »_ 



3) 



240. <p(x) = < 



A,cosV r A,+ l (Jt— *) + cos V%+ 1 re 

(X+l)cosn>^X+T 
A, cosh — X— 1 (n — x)+cosh |^ — X — 1 it 
(A,+ l)-cosh it V — X— 1 

JW+ 1, 



, ^> -1 

,*<-l 
X = -1 



241. 



<p(x) = 



3 (sinh |x + ft cosh n*)-|-sinh [x (* — 1) — 2(x cosh |x (j;— 1) 
(1 +2n 2 ) sinh n+3(i cosh|x 

fc>0 (|i = 2/ X) 
3 (sin cos [x Jc) + sin ft 1) — 2, cos \i (x — 1) 
(1 — 2n 2 ) sin n +3|x cos fx 

Jt < ([1 = 2 K^Jt) 

242. <p (x) = - 1 243. <p(x) = e-jinhVjx 

sinh V^2 +V? cosh]^2 

-sinh^cos,, ^, ((A= ^— ) 
244. cp(*) = <| sinh ■■cosh,, ^ < , ^ = 
no solutions if X = 1. 
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245. <p(x) = > 



cosh u y x —-^ 



cosh^ — ^sinh 



if H = V2X, %>0, 



COS|l(*— 

* if n = 2k, l<0 



llJl , UJl . UJl 

C0S T+T sin 2 



(f(x)=l if X = 0; u. is not a root of the equations 

, UJl UJl • t UJl n 

cosh !y — ^-sinh^- = 0, 
cosf +f sinf = 0. 
246. <p (x) = 1 + 2~^X COs8 *' ^^TT - ^° so ^ u t' ons * or 

ji 



247. <f( ;l(: ) == ^i^X x ' ^^T* ^° so l u ^' ons f° r — -j 

248. cp (x) = x + ^yjr^^ ' » ^^js- No solutions for A = ^ 

249. 3,(2^- 2 ^-5X--6) + (X + 3)^ ? ,^_ 3 Nq 

solutions for K = — 3. 

l 3(4X+5)-* . f . 3 . 3 

250. <J>(%) = < ii 3 

{ xS-llx + Cx* if k = -j 

3 

For X = — -j there are no solutions. 

25i ( x ) — / sin x if A=?M 

W\ ) — < c i cosx + C 2 sin2x + sinx if X=l. 

252. q>(je)=— 4" + T — tanh 1 if 1; 

( ) ( (u'-l)coshu* , jM 

f w \coshu— usinhutanhl ^'j u a 

if X = n 2 — 1, where u. is not a root of the equation 

cosh (i = |isinh |i- tanh 1; <P (*) - jr« [ J^ u'sin^tanh 1 ~ 1 ] 
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if X= — (n 2 + l) where \i is not a root of the equation 
cos|j, + (i sin n-tanh 1 =0. In the remaining cases there are 
no solutions. 

I— 2)x, 0<x<£ 

(i-i)x-i, s<x< i 



253. G( 



x, 6) = { 



254. It is obvious that the equation y"(x) = has an 
infinity of solutions y (x) = C under the conditions y(0) = y(l), 
y' (0) = y' (1). Therefore, Green's function does not exist for 
this boundary- value problem. 

255. Green's function does not exist. 



256. G(x, l) = 



257. G(x, I) 



258. G(x, t) = 




(31— x), 0<x'<6 
(3x—l), E<x< 1 

{x—xl + 2%), 0<x<£ 
[x(2-x)(t-2)+l], £<*<1 



259. Green's function does not exist. 260. Green's func- 
tion does not exist. 



261. G (x, l) = - 



sinh k (I — 1) sinh kx 

k sinh k 
sinh k% sinh k (x — 1) 

k sinh k 



262. G (x, |) = 



cos ( x— £ + 



2 sin 



cos ( 1—x + y 
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263. Green's function does not exist. 



264. G(x, l) = 



265. G(x, %) = ■ 



266. G(x, £) = 



267. G(x, !)=■ 



268. G(x, i) = 



269. G(x, 1) = 



270. G(x, £) = 



271. G(x, i) = 



272. G(jc, 1) = 



h + H+hH ' u 



(AE+l)[//(x-l)-ll 
h + H + hH 

a+1— j, 0<x<£ 



1 



In |_1 _£!§_!)?, 0<x<| 
4-(*-l).E<x<l 

(4)" -(f)"] , 0<x<? 



l 

2n£ 



[«>•-(!•)" 

xln| 



(lng— l)^ 
In * 
"6(lnE-l)» 



, 1 



In j, 
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273. G(x, !) = <; 0<*<£ (1^1=^1) 



2(1 +X) 2 

For A,= l, G(x, £) = — -j e_l does n °t depend on /. 
For X = — 1 Green's function does not exist. 

274. y = x — J sin* 275. y = ^(x* — 4x + 6) 

276. y = l[(l— e *)\nx + x 2 — 1] 

277. t/ = -^(2.x:— l)sinjuc. 

278. t/ = 2[sinhx— sinh(jt:— 1) — sinh 1]. 

279. y = smhx + (l — x)e x 



280. </ = 2cosx + ^2— sinx + je 2 — 2 

n 

2 

281. yW = X jG(x, £)i,(£)d£ + 



f (f-lU 0<*<£ 
where G(x, = { K % x \ 

i 

282. */(x) = a.$G(x, E).jr(E)d| + «*— + 1 

where G(*. E) = j (*_!)£, 
i 

283. y(x) = X$G(x, I) y(l) d% + ± sin^ + 1 cos ^ 

-sin -2- — — 1<*<| 



where G (*, |) = 



n 2 
isin-2.(*-S), 
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284. y(x)=-k§G(x, l)y(l)dt + ±(2x° + 3x*~\7x~5) 

where G(x, |)=<j (E _ 1)jc _i f 
i 

285. y(x) = X§G(x, I) y (£) d\ + g (x 2 - 4x + 6) 

o 

f |(3|-x), 0<x<| 
where G(x, £) = < ?2 

I V(3*-S). 

i 

y(x) = — k$G{x, t)y(l)dl + Y 2 x(x-l)(x' + x-\) 
o 

where G (x, |) = 



286. 



7x(Jt-E)(E-l). 0<x<£ 

£<X<1 



1 

287. y(jc) = c* — X Jg(jc, 6)0(£)d£ 

. r , ° J (l+x)g, 0<x<| 
where G(x, £)=<j (1+ | )jc> l 

292. cp(x) = cosx 293. (p (x) = x 2 — (cosh x — cosx) 

294. (p (x) = — ^p- x 8 + cosh 2x 

295. q,(x) = r i I5 [e«-X(e«-l)] 

296. q> (x) = jzrjj ( cos x + A. sin x) 

297. <p(x)= ! _i^ B [cosx + ^(x— sinx)] 

298. cp(x) = -f^ 299. <p (x) = V ne~* 

300. (p{x ) = -J^r + T ^ r y ^f(t)smxtdt 
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301. <p(x) = e- 



V n(\+x*) 

302. <p (x) = e x — x — . 50 1 Ox 2 — 0. 1 67 1 x 3 — 0. 0422x 4 ; 
| (p — q>|<0.18; the exact solution q>(;c)=l 

303. <p (x) =cos x + ^ [78 — 78 • sin 1 —24 • cos 1 + x (84sin 1 + 
+ 108cosl— 84)]; |<p— <}>| < 0.040; the exact solution <p (x) = 1. 

304. <p (x) = ^(e- x +3x— 1) — 0.252x 2 +0.084x 3 ; |q>— <p|< 

< 0.016; the exact solution <p (x) = x. 

one -/ \ * , 1 . , /58 16 . , 52 ,\ , 

305. <p(x) = Y + -2 smx + [ 9" — "3 sin 1 — T5 C0S J 

| <p — (p | < 0.0057; the exact solution (p(x) — x. 

306. 9W=H|x; <P.(*)=1 

307. (p„ (x) = ^ 1 — ^ <PoW = 0; the exact solution 
<p (x) = x. 

22 1 2 1 1 1 

308. q> 8 W=l + -x- Y x*--x3 + ^x*__x'__x°; 

<p (x)=l; the exact solution (p (x) = cos x -{- tan 1 -sinx. 

309. cp„ (x) = 6x 2 -f- 1 is the exact solution. 

310. <p 3 (x)=l is the exact solution. 

311. q> 3 (x)=l is the exact solution. 

312. X x = 5y (exact value ^ = 5). 

313. ^ = 2.486; ^ = 32.181. 314. X x = 4.59 
315. X 1 = 3 316. X t = 5 317. A, x = 4.19 
318. 3^ = 5.78 319. X 1= 3 320. K t = 4 
321. X 1 = 2.475 822. X l = 4.998 
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SURVEY OF BASIC METHODS FOR SOLVING 
INTEGRAL EQUATIONS 



I. Volterra Integral Equations 
Volterra integral equations of the second kind: 

X 

<p (*) = /(*) + A. S K(x, t)<p(t)dt (1) 

a 

1. Method of resolvent kernels. The solution of equation 
(1) is given by the formula 

X 

<(,(x) = f(x) + X^R(x, t; X)f(t)dt (2) 

a 

The function R (x, t; k) is the resolvent kernel of the integral 
equation (1) and is defined as the sum of the series 

R(x, i; l)= J k n K n+1 (x, t) (3) 

n = 

where the iterated kernels K n+1 (x, t) are found from the 
recursion formula 

X 

K n+1 (x, t) = ^K(x, z)K n {z, t)dz 

KAx, t) = K(x, t) (n=l, 2,...) 

2. Method of successive approximations. The solution of 
equation (1) is defined as the limit of a sequence {y„(x)}, 

n = 0, 1, 2 the general term of which is found from 

the recursion formula 

X 



It is often convenient to take the function f(x) as the zero 
approximation <f (x). 
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3. Volterra integral equations of the second kind of the 
convolution type 

<P (x) = /(*) + J *(*— t}<p(t)dt (6) 
o 

are- solved with the aid of the Laplace transformation. 
Let f(x), K(x) and q>(x) be original functions and let 

f(x) = F(p), K(x) = R(p), cp(x) = ®(p) 

Taking the Laplace transforms of both sides of equation (6) 
and utilizing the product theorem, we find 

(p(p)= F i p) , K(p)^=l (7) 

The original function q>(x) for will be the solution of 
equation (6). 

4. Volterra integral equations of the first kind 

\K(x, t)<p(t)dt = f{x) (8) 



where 

K(x, x)=£0 

are reduced, by differentiation, to Volterra integral equations 
of the second kind of the form 

o 

5. Volterra integral equations of the first kind of the 
convolution type 

X 

S*(*—/)q><0 # = (10) 
o 

are solved with the aid of the Laplace transformation. If 
f(x), K(x) and y(x) are original functions and 

f(x) = F(p), /C(x) = /C(p), <p(x) = 0(p), 
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then, by taking Laplace transforms of both sides of (10) 
and applying the convolution theorem, we obtain 



H! Kip) 



(ID 



The original function <p(x) for the function O (p) will be 
the solution of equation (10). 

II. Fredholm Integral Equations 
Fredholm integral equations of the second kind: 



(fW- t)y(t)dt = f(x) 



(12) 



1. Method of Fredholm determinants. The solution of 
equation (12) is given by the formula 

b 

<p(x) = f(x) + x\R(x, t; k)f(t)dt (13) 



where the function 



K(x, t, K) dJX) ' 



(14) 



is called the Fredholm resolvent kernel of equation (12). 
Here, 

n=l 

The coefficients B„(x, t), C„ are denned by the relations 

B„ (x, I) = 

K(x, t) K(x, I,) ... K(x, t n ) 
K{t lt t) K(t lt t,) ... K{t lt t n ) 
K(t 2 , t) K(l,, t,) ... K(t„ t n ) 



D(x, t;X) = K (x, + L B » (*• X " ( 15 > 

n=l 

CO 

°W= 1 +E L ?V (16) 



(17) 
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B (x, t) = K(x, t), 

b b 

c„=s...s x 



K(t lt U) K(t lt t t ) 
K(t t , t t ) Kit,, t t ) 
K(t 3 , /,) K(t 3 , t t ) 



x 



Kit,, t n ) 
Kit,, t„) 
Kit,, t n ) 



Kit n , tj Kit n , t. 2 ) ... Kit n , t„) 
The recursion relations are 

b 



dt t ... dt n 



(18) 



B n ix, t) = C n Kix, t)—n[K(x, s)B„_ 1 (s, t)ds, (19) 

a 

b 

C„=Sfl„-i(s, s)ds (n=l, 2,...), 



C,= l, t) = Kix, t) 



(20) 



2. The method of successive approximations. The integral 
equation (12) may be solved by the method of successive 
approximations. To do this, assume 

cp (*)=/(*) +2 (21) 

n= 1 

where ip„ix) are determined from the formulas: 

b 

^ (*) = $*(*, t)fit)dt, 
a 

b b 

ix) = [K ix, t) ^ (/) dt = S K % ix, t) f it) dt, 

a a 
b b 

= t)* t (t)dt=lK,(x, t)f(t)dt and so forth. 
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KAx, () = lK(x, z)K(z, t)dz, 



K 3 (x, t)=\K{x, z) KAz, t)dz 



and, generally, 



K n (x, t)=^K(x, z)K n - 1 (z, t)dz, n = 2, 3,. 



(22) 



where K^x, t) = K(x, t). 

The functions K„(x, t) defined by formulas (12) are called 
iterated kernels. 



3. Integral equations with degenerate kernel: 



<j>(*)-^S 



2 a* (x) b k (t) 



k=l 



<p (t) dt=f (x) 



The solution of equation (23) is 

<p(x) = f(x)+X S C k a k (x) 



(23) 



(24) 



where the constants C k are found from the linear system 

(l-HJC.-^A- • • • -u ln c n = h, 
—la 21 C 1 + ( 1 — Xa 22 ) C 2 — . . . —Xa tn C n =/„, 



•k.Ci-k ! c--+(i-k»)C„ = / J 



(25) 



Here, 



« to =S a * /.= S*«(0/(0<tt (26) 

a a 

(ft, m= 1, 2, . . ., n). 

If the determinant of the system (25) is not zero, then 
equation (23) has a unique solution (24). 

4. Characteristic numbers and eigenfunctions. The value 
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of the parameter X ^ for which the homogeneous Fri'dliolin 
integral equation of the second kind 



q> (*) — *, j/C"(JC, t)(f(t)dt = 



(27) 



has a nonzero solution <p (x) ^ is called the characteristic 
number of equation (21) or of the kernel K(x, t), while the 
nonzero solution of this equation is called the eigenfunction 
associated with the characteristic number X. 

If the kernel K,{x, t) is continuous or quadratically sum- 
mable in the square Q{a<!jt, t^.b}, then to every chara- 
cteristic number X there correspond a finite number of line- 
arly independent eigenfunctions. 

In the case of an equation with a degenerate kernel 

9 (x)— X J J 2 a k (x) b k (/)j <j> (/) dt = (28) 
the characteristic numbers are roots of the algebraic equation 



A(X) = 



1 — Xa x 

— X-a, 



— Xa l 
1 — Xa, 



—Xa 1 
— Xa, 



-Xa ni —Xa n2 \—Xa nn 



(29) 



where A (X) is the determinant of the system (25); the degree 
of this equation is p^n. 

If equation (29) has p distinct roots (1 ^Cp^n), then the 
integral equation (27) has p characteristic numbers X l ,X i , . . . , X p 
to which correspond the eigenfunctions 

<P„ (x) = 2 Cra k (x) (m = 1 , 2 p) (30) 

Here, C' ft m) (fe=l, 2, ...,n) is the solution of the system 
(25) that corresponds to the characteristic number 
X m (m = 1 , 2, .... p). For an arbitrary (nondegenerate) ker- 
nel, the characteristic numbers are zeros of the Fredholm 
determinant D(X), i.e., they are poles of the resolvent 
kernel R(x, t; X). 
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If the kernel K(x, t) is Green's function of some homo- 
geneous Sturm-Liouville problem, then finding the characte- 
ristic numbers and eigenfunctions reduces to solving the 
indicated Sturm-Liouville problem. 

5. Nonhomogeneous symmetric Fredholm integral equa- 
tions of the second kind: 

b 

(p(x)-X\K(x, 0<P (/)# = /(*), (31) 

K(x,t) = K(t, x). 

Let X„ (n— 1, 2, . . .) be the characteristic numbers and 
q>„ (x) (n = 1 , 2, ...) the associated eigenfunctions of the 
kernel K(x, t). 

(a) If the parameter X=£k n (n = 1, 2, . . .), then the inte- 
gral equation (31) has a unique solution continuous on [a, b]: 

where 

b 

a n =\f(x)y n {x)dx (33) 

a 

The series on the right of (32) converges absolutely and 
uniformly on [a, b]. 

(b) If the parameter X coincides with one of the characte- 
ristic numbers, say X = X k , of index q (the multiplicity of 
the characteristic number X k ), then equation (31) has an 
infinite number of solutions when and only when f(x) is 
orthogonal to all eigenfunctions of the characteristic number 
k k , i.e., when the q conditions are fulfilled: 

b 

\f(x)<f m {x)dx = Q (m=l, 2 q) (34) 

a 

All these solutions are given by the formula 
y(x) = f(x)-l £ 5[ ^i_ ( p nW + 

+ Cm (x) + Cm (*)+...+ C ? <P ? (x) (35) 
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where C u ...,C q are arbitrary constants, <f ■ , (a) '!',,( ») 

are eigenfunctions of the kernel associated with the chara- 
cteristic number % k . 

If even one of the q conditions (34) is not fulfilled, (hen 
equation (31) has no solutions. 

If the function f(x) is orthogonal to all eigenfunctions 
<p„(x) of the kernel K(x, t), then this function will itself 
be a solution of equation (31): <p (x) = / (x). 

6. Fredholm integral equations of the second kind whose 
eigenfunctions are classical orthogonal functions: 

(a) y (x) — % \ K (x, t) <p (0 dt = (33) 



where 

*('-/) 



, x^.t, 



The characteristic numbers are 



I 



The eigenfunctions are <p„ (x) = sin (n=l, 2, ...) 
i i 
(b) (p (x)-X j K (x, t) <p (t) dt + y f «p (/) d/ = 

where 



- 1 



K(X, t): 



1 , \+X . 
1 1 1+' 4 ^ 



The characteristic numbers are A.„ = n (n-j- 1). The eigenfun- 
ctions are <p„(.x:) = P„(;c), where P n (x) are Legendre poly- 
nomials defined by the formula 
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Since P (x)=l and P 1 (x) = x, it follows, using the recur- 
sion formula 

(n + 1) P„ +1 (x) = x (2n + 1) P n (x) + nP„, t (x), 

that it is possible to find Legendre polynomials of any 
degree n = 2, 3, ... 



(c) <p(jc) — %\K{x, t)<p(t)dt = 



where 



K(x, t) = 



2v V * / ' 



The characteristic numbers are k n — a\, where a„ are roots 
of the transcendental equation /, (a) = 0. The eigenfunctions 
are <f n (x) = J v (a n x), where J v (x) are Bessel functions of the 
first kind of order v. Bessel functions of the first kind are 
defined by the formula 



'.<*>= (4 YE 



(-')* 



* = 



k\ T(v + A:+1) 



(d) <p{x)—X S K(x, t)sp(t)dt = Q 



where 



K(x, t) = 



x+t +• 

» 2 



t 

I- co 



The characteristic numbers are 

K = n + l 

The eigenfunctions are 



(p n (x) = e" 2 L n (x) 
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where L n (x) are Chebyshev-Laguerre polynomials defined l>y 
the formula 

M*) = «* si (*"«-*) 

Using the recursion formula 

L n+l (x) = (2n+l-x)L n (x)-n*L a - 1 (x) (n> 1), 

it is possible to obtain Chebyshev-Laguerre polynomials of 
any degree n, knowing that L (x)=\ and ^(x) — — #-f- 1. 



(e) q>(x)— X J K(x,t)<p{t)dt = 



where 

X* + t' X 



K(X, t) 



y=e 2 j e-^di j e-^dx, x^t, 

-00 t 
X'+f t + 00 

e 2 J e~ z * di J e -t * dx, / ^ jc 



The characteristic numbers are 

% tt = 2(n+l) (n = 0, 1, 2, ...) 
The eigenfunctions are 

where //„(*) are Chebyshev-Hermite polynomials defined 
by the formula 

= (2jc) "_^=-L) (2^-« + «(/»-D(«-2)(«-3) (2jc)n _ 4 + _ 

Knowing //„ (x) = 1 and # x (x) = 2x, it is possible to obtain 
Chebyshev-Hermite polynomials of any degree by using the 
recursion formula 

H n+1 (x) = 2xH n (x)-2nH n - 1 (x) 
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